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FOREWORD 


The precision of a gun system clearly involves the 
dynamics of the gun carrier, ground characteristics, and 
interior and exterior ballistics. It is a problem of 
enormous complexity and is often divided into different 
phases for investigative purposes. While the division of 
the task is convenient aad often necessary, one should 
always keep in ш! 14 that the different phases interact and 
the dynamic forces are usually coupled. This fact 
neiessitetes an interactive process or, better yet, а 
complete system approach, if at all possible, to the 
precision problen. 

During recent years, one has witnessed great strides 
in various branches of continuum mechanics, kinematic 
designs, and numerical and computer techniques for solving 
problems of great complexity as well as in the areas of 
experimental mechanics and instrumentation. It appears 
feasible now more than ever to gain understanding and to 
improve the design of gun systems for greater accuracy by 
exploiting the new technological advances. The present 


Symposium represents the continuing interest of the U.S. 


Army in this direction. 





These proceedings contain nearly thirty-five papers 
presented at the Symposium held at the Hilton Inn of the 
Palm Beaches, Riviera Beach, FL, during 7-9 May 1985. The 
papers represent the current research efforts on gun 
dynamics and its effect on precision and design by 
industrial, university, and Department of Defense 
Laboratories in the United States and two allied nations - 
the United Kingdom and the Federal Republic of West 
Germany. 

The editors gratefully acknowiedge the work of Ellen Fogarty 


in preparing volumes I and II of Gun Dynamics, and her 


assistance in the collection of the papers and the required 


clearances. 








e 





LIST OF ATTENDEES 


Gary L. Anderson 

Army Research Office 

Р.0. Вох 1221! 

Research Triangle Park, NC 27709-2211 


Robert S. Becker 

S&D Dynamics 

755 New York Avenue 
Huntington, NY 11743 


Gerald A. Benedetti 
Sandia National Laboratories 


East Avenue 
Livermore, CA 94550-0096 


Willard R. Benson 

Aerojet Ordnance Company 
Rockaway Townaquare Office Park 
33 Ht. Hope Avenue, Suite 103 
Rockaway, NJ 07866 


Philip Benzkofer 
Armament Research & Development Center 
SMCAR-SCA 


һа er te 


Dover, NJ 07801-5001 


John Bostonian III 

Armament Research & Development Center 
SMCAR-LCU-SI 

Dover, NJ 07801-5001 


Robert X. Brennan 

Armament Research & Development Center 
SMCAR-LCN 

Dover, WJ 07801-5001 


Jack Brooks 

Armament Research & Development Center 
SMCAR~LCS 

Dover, NJ 07801-5001 


David N. Bulman 

Royal Military College of Science 
Land Systems Group 

Shrivenhan 

Swindon Wilts SN6 BLA England 


Bruce P. Burns 

Ballistic Research Laboratory 

AMXBR~-IBD 

Aberdeen Froving Ground, MD 21005-5066 


S. H. Chu 

Armament Research & Development Center 
SMCAR~LCA 

Dover, NJ 07801-5001 


Kok Chuug 

Armament Research & Development Center 
SMCAR-LC 

Dover, NJ 07801-5001 


Susan A. Coates 

Ballistic Research Laboratory 

AMKRR~IRD 

Aberdeen Proving Ground, MD 21005-5066 


Herbert E. Cohen 

Aray Materiel Systems Analysis Activity 
АМ SY¥-MP 

Aberdeen Proving Ground, MD 21005-5071 


P. Ae Cox 

Southwest Kesearch Iastitute 
6220 Culebra Road 

P.O. Drawer 28510 

Sau Antonio, TX 78284 


Thomas E. Devidson 

Armament Research & Development Center 
SMCAR-LC 

Dover, NJ 07801-5001 





Тһолав Dolce 

Aray Materiel Systems Analysis Activity 
AMXSY-GB 

Aberdeen Proving Ground, MD 21005-5071 


Patricia M. Dutko 

Armament Research & Development Center 
SMCAR-LCA~G 

Dover, NJ 07801-5001 


Thomas F. Erline 

Ballistic Research Lavoratory 

AMXBR- IBD 

Aberdeen Proving Ground, MD 21005-5066 


Roger К. РапсеЕ( 

Royal Armament Research & Development 
Establistir:ent 

Fort talstead 

Seven Oaks, Kent TN14 7BP 

Suz Land 


Stephen С. Floroff 
Armament Research & Development ‘ater 
СМСА DLC UF 


“> แจ ด อ ะ ม อ ฑ์ 


Dover, № 07801-5001 


Jeffrey M. Fornoff 

Armawent Research & Development Center 
SMCAR-TSB 

Dover, NJ 07801-5001 


Paul L. Fritch 

Armament Research & Development Center 
SMCAR-LC 

Dover, NJ 07801-5001 


Robert F. Gartner 
Honeywell, Inc. 

3901 South County Roa4 18 
Edina, MN 55436-5000 





Ronald Gast 

Armament Research & Development Center 
Benet Weapons Laboratory 

Watervliet, NY 12139-5000 


Daniel Gleason 

Air Force Institute of Technology 
Wright Patterson Air Force Базе 
Deytcn, OH 45433 


John Groff 

Arry Materiel Systems Analysis Activity 
ANXSY-GB 

Aberdeen Proving Ground, MD 21005-5071 


Richard G. Hasenbein 

Arnawent “3search & Development Center 
В.гес Weapons Laboratory 

Watervliet, NY 12189-5000 


Railey T. Haug 

Ballistic Research Laboratory 
AMXBR-IBD 

Aberdeen Proving Ground, MD 21005-5066 


David High 

Hercules, Inc. 

Allegany Ballistic Laboratory 
P.0. Box 210 

Cumberland, MD 21502 


Tonney Hung 

Armament Research & Development Center 
SMCAR-SCS-E 

Dover, NJ 07801-5001 


Sidney S. Jacobson 

Armament Regearch & Development Center 
SMCAR-TDT 

Dover, NJ 07801-5001 








We P. C. King 

Royal Military College of Science 

School of Mechanical Materials & Civil 
Engineering 

Shrivenham 

Swindon Wilts SN6 8LA England 


Herbert B. Kingsbury 

University of Delaware 

Department of Mechanical & Aerospace 
Engineering 

Newark, DE 19711 


Bruce Knutelsky 

Armament Research & Development Center 
SMCAR-LCA-M 

Dover, NJ 07801-5021 


George Kobori 

FMC Corporation ~ Ordnance Division 
1105 Coleman Avenue 

Box 1201 

San Jose, CA 95108 


Bent. Mat 

гъа амд зад 

Armament Research & Davelopment Center 
SMCAR-SCF-RE 


Dover, NJ 07801-5901 


Scott Langlie 

FMC Corporation 

Northern Ordnance Division 
4800 E. River Road 
Minneapolis, MN 55421 


Neil A. Lapetina 
Sandia National Laboratories 


East Avenue 
Livermore, СА 94550-0096 


Norman T. Lionetti 

Armament Research & Development Center 
SMCAR-LCW-E 

Dover, NJ 07801-5001 


Rurik K. Loder 

Ballistic Research Laboratory 

AMXBR-IBD 

Aberdeen Proving Ground, М2 21005-5066 


John M. Miller 

Harry Diamond Laboratories 
DELHD-DE-FT 

2800 Powder Mill Road 
Adelphia, MD 20783-1197 


Albert H. Mitchell 

Armament Research & Development Center 
5МСАВ-5СА 

Dever, № 07801-5001 


Lionel Pasiuk 

Naval Sea Systems Command 
SEA-62R41 
Washinzton, D.C. 20362-5101 
Edward M. Patton 

Battelle, Pacific Northwest Labs 
Battelle Road 

Richland, WA 99352 


P. H. G. Penny 

Royal Armament Research & Development 
Establishment 

Chobham Lane 

Chertsey, Surrey KT16 OEE 

England 


Robert Peterson 

Armament Research & Development Center 
SMCAR-SCS-EW 

Rock Island, IL 61299-7300 


George Pflegl 

Armament Research & Development Center 
Benet Weapons Laboratory 

Watervliet, NY 12189-5000 





8. E. Powell 

Royal Mílitary College of Science 
Land Systems Group 

Shriveuhaa 

Swindon Wilts 586 SLA England 


Robert J. Radkiewicz 

Armament Research & Development Center 
SMCAR-SCS-EW 

Rock Island, IL 61299-7300 


Frederick E. Saxe 

Armament Research & Developme Е Center 
SMCAR-LCÀ-M 

Dover, NJ 07801-5001 


Edward M. Schmidt 

Ballistic Research Laboratory 

AMXBR-IBD 

Aberdeen Proving Ground, MD 21005-5066 


Jimmy Q. Schnidt 
Ballistic Research Laboratory 
AMXBR-IBD 


E мек TS 
Aberdeen Provi 


12 Ground, Mh 21005-50656 
Mark J. Seymour 

Hunting Engineering Linited 

Reddings Wood, Ampthill 

Bedfordshire MK45 2HD England 


Thomas E. Simkins 

Armament Research & Development Center 
Benet Weapons Laboratory 

Watervliet, NY 12189-5000 


Martin T. Soifer 

S&D Dynamics 

755 New York Avenue 
Huntington, NY 11743 





G. D. Stilley 
Honeywell, Inc. 

5901 S. County Road 18 
Edina, MN 55436-5000 


Martin Thomas 

Royal Armament Research & Development 
Establishment 

Chobham Lane 

Chertsey, Surrey KT16 OEE 

Eng land 


Philip E. Townsend 

Armament Research & Development Center 
SMCAR-SCS 

Dover, NJ 07801-5001 


Tiea-Yu Tsui 

Army Materials & Mechanics Research 
Center 

Watertown, MA 02172 


David C. Uhrig 

Air Force Armament Laboratory 
Guns and Projectiles Branch 
AFALT/DLIG 

Eglin Air Force Base, FL 32542 


John D. Vasilakis 

Armament Research & Development Certer 
Benet Weapons Laboratory 

Watervliet, NY 12189-5000 


Cris Watson 

Ballistic Research Laboratory 
AMXBR-IBD 

Aberdeen Proving Ground, MD 21005-5066 


Emma M. Winebolt 

Ballistic Research Laboratory 
AMXBR-IBD 

Aberdeen Proving Ground, MD 21005-5066 








E 
` 
Ze A 


Julian J. № 


\\ US Army Research, Development, & 


“. 


Standardization Group (UK) 
Box 65 
FPO, NY 09510 


Adam Zak 

Aeronautical & Astronautical Engr Dept. 
University of Illinois 

104 S. Mathews Avenue 

Urbana, IL 61801-2997 


William T. Zepp 

Armament Research & Development Center 
SMCAR-LCW-S 

Dover, NJ 07801-5001 


xiii 


John В. Zimmerman 

Armament Research & Develc ment Center 
SMCAR-TSE-L 

Dover, NJ 07801-5001 


Kathleen L. Zimmerman 

Ballistic Research Laboratory 

AMXBR-1BD 

Aberdeen Proving Ground, MD 21005-5066 








sr —————— —— — — ұла зы қыл ФИНИ даладан ы қ лл — uM EUM ЧАЛ ЧЕН ЧЫ аллы қалалы A Ш aS P VL eee ed REN. M SN s. ИМА, a. L... ` 


TITLE: Technology Review on Projectile Gun Disengagement 
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Dr. Roger K. Fancett 
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ABSTRACT: 


This paper outlines the importance of an accurate projectile launch 
definition to gun accuracy, projectile-guu dynamics models, and gun system 
diagnostics; reviews the technology which is available for the investigation 
of the projectile disengagement mechanism from the gun barrel; discusses the 
relevant instrumentation and measurement methods; and describes the ongoing 
research efforts, thus providing the introduction and basis to the other 
papers. 


This paper was unavailable at the time the Proceedings were published, For 
further iaformation, contact the authors. 
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Jimmy Q. Schmidt 
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ABSTRACT: 


Acquistion of data pertaining to the launch parameters of a projectile ís 
very important for understanding the underlying causes of projectile-gun 
system inaccuracy. These experimental data are also necessary in the 
validation of the theoretical modeiling of projectile-gun systems, During the 
past several years, the Ballistic Research Laboratory has been developing a 
new technique to measure projectile and gun tube motion during the launch of 
the projectile. As part of the ongoing program, the Ballistic Research 
Laboratory, US and the Royal Armament Research and Development Establishment, 
UK, recently conducted a joint collaborative firing experiment with a 40-mm 
gun. The primary purpose of the experiment was to assess the BRL system by 
performing paraliel measurements using standard electrooptical methods. This 
paper presents a brief description of the measurement techniques used, the 
firing exneriment, and a brief discussion of the preliminary results. 
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DESCRIPTION OF THE JOINT BRL-RARDE 40-ММ EXPERIMENT TO DEFINE PROJECTILE 
LAUN Н 


MR. JIMMY Q. SCHMIDT 
U.S. ARMY BALLISTIC RESEARCH LABORATORY 
ABERDEEN PROVING GROUND, М 21005-5006 


MR. THOMAS O. ANDREWS 
ROYAL ARMAMENT RESEARCH AND DEVELOPMENT ESTABLISHMENT 
FORT HALSTEAD, UK 


1. INTRODUCTION 


In May 1982 a new measurement technique [1] which provides measurements 
of projectile transverse displacement during muzzle exit was presented by the 
Ballistic Research Laboratory (BRL), at the Third U.S. Army Symposium on Gun 
Dynamics. This new technique is a four sensor coil arrangement derived from 
the Radio Frequency Oscillator (RFO) velocimeter commonly referred to as the 
"Muzzleschmidt". This method aroused considerable interest among the visiting 
delegation of the Royal Armament Research and Development Establishment 
(RARDE), from the United Kingdom (UK), because of its high potential of 
providing the experimental means for defining the projectile launch. 


The importance of accurately defining projectile launch parameters and 
some of the previcus experimental methods used to obtain such data are 
discussed in the opening paper of the Session," iechnology Review on 
Proj3ectile Disengagement” by Dr. R.K. Loder (BRL) and Dr. Roger K, Fancett 


(RARDE) . 


During subsequent discussions between BRL and RARDE personnel it was 
evident that both the US and the UK would benefit from a joint BRL/RARDE 
Program to assess the ability and accuracy of the Schmidt displacement 
transducers to measure both the projectile motion with respect to the gun 
muzzle and the gun muzzle motion with respect to the ground during the 
projectile-barrel disengagement. BRL would provide measurements using the 
Muzzleschmidt technique while the UK made measurements using their optical 
methods. By performing the test in the UK, RARDE personnel would be able to 
obtain a direct comparison of the measurement techniques and to assess the 
relative merits of each. The US while also obtaining the direct comparison 
would in effect obtain a dynamic calibration of the BRL method. This would be 
highly desirable since the sensitivity of the RFO sensor to mechanical 
displacement is at least an order of magnitude greater than what is easily 
obtained with a mechanical test calibrator. 
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In May 1983, Mr. J.Q. Schmidt, BPL, visited RARDE to coordinate with Mr. (59 
T.O. Andrews and consider all the technical details and the time table for 
preparing and conducting the joint firing experiment. This joint experiment 
was carried out at RARDE in June 1984. This paper covers the experimental 
aspect of the joint collaborative firing experiment. A paper on the data 
analysis will be given later in this Session. 


The objective of this trial, fired in RARDE/GR2's range at Fort Halstead 
using a 40-mm gun, was to compare the results obtained from Radio Frequency 
inductive loop devices, developed at BRL, with the results of observations 
using optical devices developed in the UK by Hunting Engineering Ltd. (HEL) 
and by RARDE. 


2. GENERAL OUTLINE OF THE TEST 


A total of twenty-one 40-mm 160 ogival nosed projectiles were fired at 
velocities of approximately 500 m/sec, using the RARDE Ord. Q.F. 40/70 gun. 
One half of the projectiles were essentially unmodified and the other half 
deliberately unbalanced to provide a yaw in excess of the maximum yaw normally 
encountered with this projeccile-gun system. Twelve of the rounds had a 0.254 
mm undercut from 5,00 mm behind the bourrelet to 5.00 mm in front of the 
rotating band. This undercut is used to accurately relate longitudinal points 
on the recorded waveform corresponding to longitudal points on tbe projectile. 


BRL provided the instrumentation for the measurements of the displacement 
of the projectile with respect to the bore axis during shot exit using che o 
Muzzleschmidt technique. Measurements were made at two locations, one at the 
muzzle, and one 2.54 cm forward of the first sensor. The use of dual sensors 
allows the rate of yaw determination. BRL also provided a gun tube motion 
measurement with respect to ground during shot exit, vsing a modified 
Muzzleschmidt detector. Shot exit time was obtained by the BRL projectile 
displacement sensor. Since the gun tube motion after shot exit may exceed the 
clearance between the gun tube and the tube motion sensor, RARDE built a 
mechanical sliding platform to pull the gun tube motion sensor clear of the 
gun after shot exit. 


RARDE provided measurements of gun tube motion at two stations some 
distance to the rear of the BRL sensor mounting collar, using the HEL 
electrooptical device incorporating a collimated light beam and knife edges. 
During the latter stages of the test RARDE measured the projectile in-bore 
yaw, using Peter Fuller's iaser and optical grating method. In addition to 
the specific test requirements, RARDE also measured breech pressure, 
projectile velocity, and target strike. RARDE provided the data recording 
facility. Projectile displacements were recorded on Nicolet digital 
oscilloscopes with backup recording by a high speed digital data logger. Tube 
motion measurements, both BRL and КАКОЕ, plus skot exit time were recorded on 
a second data logger with backup recording on an analog tape. Data recorded 
on the data loggers were transferred to the Nicolet oscillcscope for storage 
on floppy disc. This enabled BRL to hand carry the recorded data back for 
processing. The RARDE in-bore yaw data were recorded en film. 


One of the prime areas of concern was the mechanical compatibility which © 
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is quite often taken too lightly. When fabricating parts to fit a gun tube 
based on mechanical drawings, quite often problems are encountered such as 
improper thread mating and unexpected variation in tolerances, To eliminate 
any last minute problems, RARDE fabricated the mounting collar in the UK based 
on BRL requirements, assured proper fitting of the collar to the gun tube and 
sent it to BRL for the sensor and electronics to be tastalled, 


3. MEASUREMENT TECHNIQUES 


As previously mentioned, BRL's measurements are based on the RFO 
(Muzzleschmidt) techníque while RARDE's measurements on optical methods, Each 
method has its relative merits and, by conducting parallel measurements, a 
better assessment of each method is obtained. [п addition, in view of the 
information exchange program between the US and the UK, the expertise gained 
by the parallel test will he valuable in future information exchanges. 


A. BRL, Dual Projectile Displacement Measuring System (Muzzleschmidt) 


Several years ago the RFO technique was developed by BRL to provide a 
method for measuring projectile velocity directly at the mvzzle, primarily on 
tapid fire guns or moving gun systems. The basic concept is quite simple, an 
inductive loop excived by a radio frequency oscillator and clamped to the 
muzzle face detects the metallic parts of the projectile as it passes through 
it [2]. The detected signal is an electrical pulse representative of the 
geometric configuration of the projectile. Knowing the projectile length, the 
measurement of the tíme duration of the detected pulse provides a measurement 
of muzzle velocity. Experience gained in the development of the RFO 
velocimeter lead to the development of the RFO system tc measure projectile 
transverse displacement during muzzle exit [3]. The basic operation of the 
system is explained in detail іп the reference so only a very brief 
explanation vill be given here. 


In the RFO displacement measuring system, the single circular sensor used 
in the velocity measuring system Is replaced by a four segment sensor which is 
configured to provide two semi-circular loops in the horizontal plane and two 
semi-circular loops in the vertical plane as shown in Figure 1. Each sensor 
Joop is responsive to the proximity of the metallic parts of the projectile to 
it. By electronically differencing the signal from the left and right loops a 
measurement of the projectile displacement from the center axis is obtained as 
a function of time. Vertical projectile displacement is obtained in the 
identical way. The displacements from the center axis in the horizontal and 
vertical plane as a function of time yield the angle and the orientation of 
the projectile yaw. By electronically summing the detected signals, a pulse 
is obtained from which velocity can be computed or used simply as a very 
precise muzzle exit time. 


Shown in Figure 2 is the basic block diagram of the sensors and tle 
summing and differencing circuits. The detected waveforms shown correspond to 
the detected signal of a cylindrical projectile perfectly aligned in the 
horizontal plane. Should the projectile be displaced from center but have no 
yaw, a positive or negative pulse would be obtained, the polarity being 
dependent on the direction of the displacement from the axis and the magnitude 
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Figure 1. The Basic Sensor Oil Configuration for Measuring Projectile 
Displacement 


corresponding to the amount of displacement. The signal shown in the vertical 

channel represents the projectile ' ith a fixed yaw in the vertical plane as 

indicated by the linear slope of the pulse. If a rate of yaw were present at 

the time the projectile was passing the sensor, the slope would be non- 

linear. The summed signal is used to provide a profile of the pulse to assess 

the acquired data, obtain muzzle velocity, and for use as a muzzle exit 

trigger. The summed signal is used to trigger the recording device. This is 

necessary because the differenced signals can be either positive or negative- 

Boing or even essentiallv zero, depending on the orientation of the projectile 

as it passes through rhe sensor. © 


By using two sensors, spaced a known distance apart, the rate of yaw can 
be obtained from the simultaneous displacement measurements of the front and 
rear section of the projectile. This of course is an oversiaplification of 
the complexity of the analysis, when one considers the complex motion caused 
by the longitudinal displacement, the spin of the projectile, the rate of yaw 
during disengagement, gun tube motion, aud minute sensor to gun displacements 
caused by stress waves being propagated down the gun tube. A complete 
mathematical analysis of the projectile motion will be given in the data 
analysis paper to follow this paper. 


Shown in Figure 3 is the block diagram of the dual projectile 
displacement measuring system. The sensor and amplifier assemblies are 
contained in two electronic housing boxes mounted on the collar at the gun 
muzzle. The horizontal and vertical difference signals plus the horizontal 
summed signal from each sensor are coupled to the recording facility via co- 
axial cables. Unity gain isolation amplifiers in the recording facility are 
used to minimize 50 Hz ground loop currents. The horizontal and vertical 
difference signals from each sensor were recorded on Nicolet digital s: 
oscilloscopes, using the horizontal summed signal as a trigger. A. 


B. BRL, Gun Tube Motion Measurement System 
Since the measurements of projectile displaceme = are with respect to the 


gun tube on which the sensors are mounted, gun tube motion measurements ey 
relative to ground must be made in order to accurately determine projectile b 
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Figure 2. Block Diagram of the Signal Processing Circuit 
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Figure 3. Block Diagram of the Dual Projectile Displacement Measuring System 


motion in reference to ground. The RFO technique can also be used to measure 
gun tube motion. By fabricating identical sensors and circuitry as used to 
measure projectile displacement and mounting them in a metal loop surrounding 
the gun tube, accurate measurements of tube motion can be obtained in the 
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horizontal and vertical planes. Ideally, the clearance between the gun tube 
and the sensor ring should be large enough that the gun tube would not strike 
the sensor when counter -recoiling to the rest position. If the spacing 
between tne sensor ring and the muzzle face, or in this case the mounting 
collar, is increased, the radiated field depth is greater, permitting the 
sensor tc be further from the gun. This causes a loss of definition of the 
signal which is undesirabie for projectile masurements; however, for tube 
motion measurements an accurate definition of the longitudinal displacement of 
the gun is actually not necessary. This permits increasing the spacing 
between sensors for gun tube notion measurements. By placing two additional 
rings of printed circuit (PC) substrate between the horizontal and vertical 
sensors and one between each sensor and the metal mounting collar, the 
radiated field depth was increased to provide measurements with a sensor ring 
which is 10.0 mm larger in diaweter than the gun, 


Electrically, the sensor and associated circuitry were identical with 
that used to measure projectile displacement with the exception of a lesser 
signal bandwidth requirement. Mechanically, the sensor consists of a 
vertically oríented sensor loop mounting attached to a small box housing the 
electronics. The box was mounted to an adjustable base providirg horizontal 
and vertical adjustments to position the sensor loops equidistant from the 
outer surface of the gun. Shown in Figure 4 is the basic configuration used 
for gun tube motion measurements. 


— e —n s 





ELECTRONICS 
BOX 


ADJUSTABLE BASE‏ سح 


Figure 4, The Basic BRL Configuration Used for Gun Tube Motion Measurements 


Since this was the first test of the RFO technique to measure the 
transverse motion of the gun tube at projectile exit and the required 
sensitivity was unknown, it was decíded that 5.0 mm clearance would be 
sufficient for tube displacement before and immediately after muzzle exit. 
However, provisions had to be made to pull the sensor forward several 
centimeters after projectile exit so that the tube , after recoil, did not EY 
strike the sensor, It may be possible to increase the diameter of the sensor 
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rings sufficiently to prevent interaction between the tube and the sensor in 
future tests. x 


C. RARDE, Gun Tube Motion Measurement Technique 


One optical method to measure gun tube metion is the use of a collimated 
light beam passing over a knife edge and partially illuminating a large area 
light sensitive diode through a narrow slit [4]. To measure transverse gun 
tube motíon, kuife edges (razor blades) are mounted to the gun tube parallel 
to the tube axis and perpendicular to the light beam, The razor blade is 
positioned to partially Llock the light beam passing through the slit masking 
the photo-diode. If the razor is mounted in a vertical posit on, partially 
obscuring the light beam, than any vertical tube motion will move the razor 
blade and cause a lesser or greater amount of light to strike the sensor. The 
detected signal is amplified and recorded. The signal is proportional to the 
gun tube motion in the vertical plane. By using two light sources and sensors 
positioned orthogonally, both the horizontal and vertical components of tube 
motion can be measured. The basic test set up used to measure tube motion is 
shown in Figure 5. 


LIGHT SOURCE FA 
SENSORS (VERTICAL) 
) Q BARREL 


BRL SYSTEM 


BREECH 














p EDGES 
SENSORS (HORIZONTAL ) LIGHT SOURCE 
Figure 5. The Basic Electrooptical System To Measure Gun Tube Motion 


D. RARDE, Projectile In-bore Yaw Measurement System 


This system determines the in-bore yaw of a projectile by the deflection EN 
of a parallel beam of laser light reflected from its nose [5]. The light from s 
a l6mW Неће laser (wavelength 0.000632 mm) is focused by a lens fitted to the 
laser to give a point source of light in the focal plane of a concave ті. сог 
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300 am in diameter (focal length 1220 шп, 48 inches). The lens is chosen so QJ 
that the mirror produces a beam of parallel light about 25 mm in diameter. 


This beam of light, Figure 6, is reflected by plane mirrors to lie alorg 
the tore axis. It is then reflected back by & reflective optical grating 
fixed to the nose of the projectile and set to be accurstely perpendicular to 
the projectíle's axis. The grating has equal widths of reflecting and non- 
reflecting surface, the reflecting strips being spaced at intervals of 
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= 
ORL PROXCTLE MOTION LOOPS, © ы 
INSIDE MUZZLE ATTACHMENT 


BRL ТИШЕ MOTION’ LOOP 


Figure 6. The Complete Instrumentation Layout 


0.4233 mm. The area of the grating is divided into two, with the grating 
lines оп one half being perpencicular to the lines on the other half. The 
reflected light then consists of a central beam and four beams representing 
the two first-order diffraction beams from each half of the grating. When 
these beams, which return via the plane mirrors to the concave mirror, are 
brought to a focus they produce five spots in the form of а "quincunx" (5-spot 
die pattern). 


A half-silvered mirror is placed between the laser and the concave mirror 
to act as a beam splitter, This reflects the returning beams away from the 
laser to where the spots can be recorded on a cíne film, at about 6000 frames 
per second. The ends of fibre optic iight guides are positioued in the focal 
plane of the camera to provide reference points for film reading cnd film 
orientation. 


The angular divergence of the beams producing the spot pattern depends 
only on the grating spacing and on the laser wavelength and so provides an "e 
automatic yaw angle calibration. With the laser and gratings used for this La 


111-59 











SCHMIDT, ANDREWS 


trial, each outer spot is at an angle of 0.0855 degrees from the central spot; 
but, as a yaw angle of 0.1 degrees will deflect the spot pattern as a whole by 
0.2 degrees, the yaw scale is obtained by taking the diagonal pairs of outer 
spots as being 0.0855 degrees apart. 


If the light beam reflected down the bore is truly parallel, then 
deflection due to yaw, as viewed in the focal plane of the concave mirror, is 
independent of the projectile's travel. Any change in scale which may occur 
will be shown by a corresponding change in the spot pattern size. 


E. Data Recording 


The prime consideration in choosing the method of recording the trials 
data was the need to provide copies of the data for BRL in a form suitable for 
analysis in the United States, and preferably available to be taken to BRL 
immediately at the end of the trial. This requirement involved more effort 
than all other aspects of RARDE's preparation for the trial. 


The recording system assembled was largely digital, with an analog tape 
recorder as back-up for the low-frequency records, i.e., barrel motion and 
chamber pressure, with one projectile profile as a synchronizing signal. 

Three Nicolet Explorer III two-channel digital oscilloscopes were used to make 
the primary high-frequency recordiags of the outputs of the Schmidt 
displacement transducers, with a Datalabs DL2800 transient recorder as back- 
up. The primary low-frequency records were made using an eight channel 
Datalabs DL1200 transient recorder. 


A Digital Equipment Corp. MINC 11/03 computer was used to control the 
transfer of data from the Datalabs recorders into one of the Nicolet 
oscilloscopes and the recording of all the data on Nicolet disks. This could 
readily be done in a few minutes between rounds. As opportunity offered the 
data were copied from the Nicolet disks onto MINC disks for later analysis at 
RARDE so that the Nicolet disks were available to be taken to BRL for 
analysis. 


4. DETAILS OF THE BRL SENSORS 


As agreed, preliminary and final drawings of the mounting for the 
projectile displacement sensor were exchanged between BRL and RARDE. The 
mounting cullar was then fabricated by RARDE. During the time the collar was 
being fabricated, the electronics were constructed as far as possible without 
the sensor assembly. Upon receipt of the collar, sensor rings were fabricated 
by BRL. Each sensor consist of three rings constructed from PC substrate, 
Approximately 0.6 mm of copper was removed from the ID of the rings and the 
adjacent ring undercut by the thickness of the copper clad as shown in Figure 
7. This placed the ID of the sensor rings behind a protective layer of PC 
substrate and provided a seal when epoxied together. This protective layer 
was to prevent any metallic parts of the rotating band from striking the 
sensor elements. After the sensor loops were formed from the copper clad and 
leads attached, the sensor rings мете epoxied into the sensor housing апа held 
in place by rings which screwed into the housing. 
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Figure 7. Projectile Displacement Sensor Rings Before Sensor Loop Formation 


The complete sensor mounting collar was then assembled and connected to 
the oscillator circuitry. The electronics for the muzzle sensor were housed 
in a box on one side of the collar, while the electronics for the front sensor 
were housed in a box on the opposite side. 


The muzzle attachment consisted of a collar screwed onto the muzzie with 
a ring bolted to the forward surface containing two sets of BRL loops, one 
close to the muzzle face and the cther 25.4 mm further forward, for the 
observation of the position of the projectile. A further ring bolted to the 
front of the attachment provided a measuring surface for the “tube motion” 
loop. 


Shown in Figure 8 is a drawing of the collar assembly and the relative 
positions of the sensor elements and the gun tube. 


The complete projectile sensor assembly was electrically checked out and 
mounted on a dummy gun tube. Sensitivity checks were made using a mechanical 
calibrator consisting of a steel cylinder with the diameter of the projectile, 
which was inserted into the sensor area. The steel cylinder was displaced 
from the center line of the bore axis іп 0.127 mm increments and measurements 
were made in the radial direction every thirty degrees by inserting the face 
plate of the calibrator into guide piu holes. The following sensitivities 
were measured. 


Muzzle Sensor Front Sensor 
Left - Right 8.6 V/mm Left - Right 8.32 V/mm 
Up ~ Down 8.58 V/mm Up - Down 7.76 V/um 


With sensitivities in the order of 8.0 V/mm it is easy to see the need for 
an extremely precise mechanical calibrator which at this time is not 
available. Therefore, at the present time, the calibration constants are 
obtained via statistical data analysis. More precise calibration relations 
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SENSOR RINGS: 2PAIR , 3EA- 1.58 тт THICK 


Figure 8. Relative Positions of Projectile Sensors and Gun Tube 


based on the statistical analysis are used in the actual data analysis of the 
recorded data. 


Shown in Figure 9 is the complete projectile sensor assembly before 
mounting to the gun tube. The assembly replaces the flash hider which is 
normally screwed onto the gun. 


The tube motion sensor performs electrically identical to the projectile 
sensor; however, in this case the sensor surrounds the muzzle collar rather 
than the emerging projectile. Since longitudinal definition is not necessary 
for tube motion measurements, and a greater spacing is required between the 
sensor and the gun, the spacing between the sensor elements was increased. 
Instead of three PC rings with a one PC ring spacing between elements as used 
for the projectile sensor, six PC rings were used with the sensing elements 
spaced two PC substrate thicknesses apart. This provided a usable depth of 
the radiated electromagnetic field of at least 6.00 m. 


The tube motion sensor was mounted to the dummy gun tube fixture with the 
projectile sensor collar. Calibration was performed by recerding the 
horizontal and vertical output voltages while repositioning the sensor 
assembly with its horizontal and vertical adjustments on the base. Again this 
is a rather crude method to calibrate, but at this time no other calibration 
rig is available. The horizontal sensitivity was found to be 5.2 V/mm and the 
verttcal sensitivity was 4.8 V/mm, These sensitivities were with reduced gain 
in the amplifier. Additional gain can be provided by the implifier but this 
would limit the dynamic range of the measurement because of the amplifier 
output Limiting due to the voltage swing. This sensitivity provides a dynamic 
range of tube motion of approximately + 2.5mm. Figure 10 is a photograph of 
the cospleted tube motion sensor. 


111-62 








SCHMIDT, ANDREWS 


BRL, Projectile 
BRL, Projectile 
BRL, Projectile 
BRL, Projectile 
Sensor (Front 
BRL, Projectile 
(Front U-D) 
BRL, Projectile 
Muzzle Sensor 


BRL, Horizontal 
BRL, Projectile 
RARDE, Vertical 


Horizontal) 


Horizontal) 








Figure 9. Complete Projectile Sensor Assembly 
5. DETAILS OF THE TEST 


The actual firing tests were conducted at the RARDE cest range during the 
two week period of 29 May through 12 June 1984. Interfacing the BRL equipment 
to the RARDE facility proceeded extremely well with no problems being 
encountered either mechanicaily or electrically. The following data were 
recorded, the longitudinal position given here being only approximate, 


Horizontal Displacement at the Muzzle (Muzzle L-R) 
Vertical Displacement at the Muzzle (Muzzle U-D) 

Summed Horizontal Displacement at the Muzzle (Muzzle L+R) 
Horizontal Displacement 2.54 cm forward of the Muzzle 

L-R) 

Vertical Displacement 2.54 cm forward of the Muzzle Sensor 


Summed Horizontal Displacement 2.54 cm forward of the 
(Front L*R) 


BRL, Vertical Tube Motion at the Muzzle (BRL U-D? 


Tube Motion at the Muzzle (BRL L-R) 
Muzzle Exit 
Tube Motion, 32.5 cm from the Muzzle (HEL Front Vertical) 


RARDE, Horizontal Tube Motion, 32.5 cm from the Muzzle (HEL Front 


RARDE, Vertical Tube Motion, 44.5 cm from the Mizzle (HEL Rear Vertical) 
RARDE, Horizontal Tube Motion 44.5 cm from the Muzzle (HEL Rear 


RARDE, Breech Pressure 

RARDE, Projectile Velocity 

RARDE, Projectile In-bore Yaw 

RARDE, Projectile Strike Coordinates at the Target 
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@ In addition to these primary recordings, all projectile and tube motion 
measurements were put on backup recording devices to prevent loss of data. 


Shown in Figure 11 is the BRL and RARDE instrumentation attached to the 
40-mn gun. The laser set up to measure in-bore yaw is not shown in this 
photograph. 


Nineteen rounds were fired during the two week test period with two 
additional rounds being fired by RARDE, after the BRL personnel left, to 
obtain additional projectile in-bore data with Fuller's optical lever/grating 
system. The projectiles were all fired at a velocity of approximately 525 
m/sec. As previously mentioned some rounds fired were standard projectile 
configurations and some rounds were deliberately unbalanced, Details of the 
gun and projectile configuration fired are given below. 


The gun and ammunition used were chosen because a considerable volume of 
data had been obtained for these firing conditions in previous trials. 


Gun: 40-mm 170, Barrel No. E3169, with 1 in 18 twist of rifling, 
muzzle face machined to remove chamfer, BRL sensor assembly P 
fitted to muzzle in place of "Flash Hider." Adapted into 6 Pdr. га 
7 cwt breech rings and recoil mechanism mounted on range stand. E 


Charge: All, 95 grams SC/Z, 0.02 mm scroll propellant. 
@ Ignition: No. 12 M4 Percussion. 
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Figure li. BRL and RARDE Instrumentation Attached to the 40-mm Gun 


Projectiles: 


All 


Type A 


A shell body, 40 ша/Т Ж 4, Part №. SX252, fitted with a steel 
base plug of local manufacture in place of a tracer. 


Fitted with an ogival aluminum alloy nose plug and filled with 
88 grams of H.E. Substitute wax giving a total weight of 
approximately 865 grams. The center of gravity is 64.6 ши from 
the base of the shell. 


ALUMINIUM AUCH 
NOSE MUG 


As A, but with the outer surface of the body turned down to form 

an "undercut" 0.25 om deep from 5 mm behind the shoulder to 5 mm 

in front of the driving band, a length varying between 34.62 mm 

and 35.05 ши. This reduces the shell weight by 8 grams and 

moves the center of gravity 0.17 mm towards the base. » 
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As A, but deliberately unbalanced, with a 9.53 mm (3/8") 
diameter steel rod, 77.4 mm long, placed against one side of the 
shell cavity, the cavity then being filled with Bees Wax, of 
density 0.95 gram/ml, rather than H.E. Substitute wax of 
density 1.674 gram/m;. The rod weighs 43 grams, and the Bees 
Wax filling 45 grams, giving the same weight of filling (88 
grams) and the same total shell weight. With a cavity diameter 
of 27.53 mm and a rod diameter of 9.53 mm the rod axis is 9 mm 
from the shell axis. The center of gravity of the complete 
shell is 0.39 mn off axis and possibly 0.1 mm nearer the base 
than in type A. 


QuiSt! маз NOD #5) = f. 71 4mm LONG 
PLACED AGAINST МОЕ C SHELL Cavity 


ROO ARIS a fen OFFSET PROM Bi\ 
& Of G ARIS = 0. Nom OFFSET FROM vu Vi 





WI 


With a machined outer surface as B, and unbalanced as С. 


As A, but with the aluminum alloy nose plug replaced by a flat 
steel nose plug fitted with a reflective optical grating for in- x 
bore yaw measurement. These shells also had a second groove j 
machlned in the body, behind the crimping groove, to take a 

rubber obturating ring, to reduce the obscuration by gas leakage 

of the light reflected from the optical grating. The total 

shell weight is still 865 grams but the center of gravity is 

possible 0.6 mm nearer the base than in type A. 


ым ни TO SMOOT: PRONI? 
ERECTIE OPTICAL 
GRAING 


With a machined outer surface as B, and with an optical grating 
and obturating ring as E. 


With a machined outer surface as B, and unbalanced as C, with an 
optical grating and obturating ring as E. 
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6. PRELIMINARY RESULTS 


Since the analysis of the recorded signals for the projectile and tube 
motion is quite complex and covered in detail in the following papers in this 
sessicn, no attempt is made to present precise data analysis or assessment of 
the comparative results. However, preliminary assessment of several rounds of 
the recorded data gives some insight on the quality of the data recorded as 
well as the general trends in the behavior of various configurations of the 
projectile fired. 


Comparison of the recorded data of tube motion seems to indicate that the 
tube motion of similar rounds is fairly consistant. The unbalancing of the 
projectile mass distribution produces a distinctly different motion near 
Muzzle exit yet is very similar within its particular grouping. Unbalancing 
the projectile in the opposite direction produces an opposite effect but still 
similar within its particular grouping. Since only six rounds were fired with 
the center of gravity (00) offset in the vertical plane (three with the G; up 
and three with the GG down), there is not enough data to come to any firm 
conclusion. Even with a greater number of data rounds statistical analysis is 
necessary before any real degree of confidence can be placed on the results. 


This is also true of the analysis of the projectile motion. Quick look 
comparison of the projectile motion of the unbalanced rounds shows some 
similarity within their pariiculat grouping (three rounds with the G offset 
up and three rounds with the G; offset down). However, examination of the 
records of the balanced rounds shows a wide variation in the epparent amount 
of projectile yaw. Since these measurements are relative to the gun tube and 
the gun tube motion during the launch phase of the projectile is distinctly 
different for balanced rounds versus unbalanced rounds, complete analysis is 
again necessary. 


Data from the single BRL tube motion sensor provides only transverse tube 
displacement and not the angle of the tube axís relative to the line of 
sight. In general, however, the transverse tube motion during and immediately 
after projectile exit using balanced rounds is small as compared to when 
unbalanced rounds are fired. 


Shown in Figure 12 are the recorded difference signals of projectile 
displacement at the muzzle sensor for a balanced round with a slight undercut 
from 5.00 mn behind the bourrelet to 5.00 mm in front of the rotating baud. 
The recorded vertical signal indicates that the projectile is well centered in 
the vertical plane with virtually no indication of projectile yaw relative to 
the gun tube. The recording of the horizontal displacement indicates the 
projectile is displaced from the center axis but, as in the case of the 
vertical signal, with very little indication of projectile yaw. It should be 
noted that these measureuents are made with respect to the center axis of the 
sensors and there can be a slight displacement between the center axis of the 
sensor and the center axis of the gun tube. This will be taken into account 
in the actual data reduction. The signal is clearly defined from the front of 
the bourrelet to the rotating band. Figure 13 shows the horizoatal summed 
signal of the same round. The undercut, the rotating band, and the groove for 
crimping the case to the projectile are clearly defined. 
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Figure 12 Figure 13 


The recorded signal of the projectile displacement at the front sensor 
for the same round, Figure 14, shows a substantial amount of yaw in both the 
horizontal and vertical nlanen. Since this уаш interpretation 18 relative to 
the gun tube and the front of the projectile has disengaged from the gun, it 
must be resolved as to whet’.er the amount of yaw is completely projectile 
related or whether the muzzle of the gun is moving and affecting the 
trajectory of the projectile. Both sets of information are important to gun 
tube accuracy, since, even if the projectile was perfectly aligned with the 
gun tube axis but the tube axis has actually moved, the trajectory will be 
affected. 


Shown in Figure 15 is the front-summed signal. Here again the 
longitudinal reference points such as undercut, rotating band and groove are 
very clearly defined. The difference in signal amplitude in the rotating band 
area can be attributed to the relative response of the sensor to the copper 
band as compared to the steel projectile case and any asymmetrical engraving 
of the band. These are only typicai records of the test and direct 
correlation between the projectile yaw and the location of the projectile 
center of gravity or gun tube motion will have to be resolved in the complete 
analysis of data. 


A rough indication of the angle and the angle of orientation of yaw can 
be obtained from this quick look data. The data were recorded on the + 4 volt 
range or 8 volts total range of a Nicolet oscilloscope. The 12 bit resolution 
of the Nicolet oscilloscope provides approximately 2 mV per point, therefore 
the measured voltage can be read from the plot. A more accurate quick look 
method is to read the voltage differential of the slope frcm the Nicolet 
recordings. This provides a voltage differential, from the beginning of the 
undercut to the end of it, of 668 mV for the horizontal and 678 mV for the 
vertical signal. This agrees roughly with what can be read from the plot in 
Figure 14 when plotted to this scale. The length of the undercut on this 
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prolectile was 35.05 mm. Taking into account the substrate thickness of the 
sensor, the recorded values were measured after the raised shoulder of the 
projectile had passed the sensor by the thickness of the substrate in front 
and before the rear raised shoulder passed ty sensor by one substrate 
thickness. Therefore the baseline is actually 35.05 mm minus the thickness of 
two substrates (3.175 mm) or 31.875 mm. Using a horizontal sensitivity of Ө 
8.32 V/mm and a vertical sensitivity of 7.76 V/mm. The horizontal and 
vertical component calculates to be 0.0813 mm and 0.0889 mm. By transposing 
these values to zero and solving for the vector the displacement calculates to 
be 0.120 mm down and to the left, looking downrange, at an angle of 42.4 
degrees or 222.4 degrees from the tweive O'Clock position. The actual angle 
of yaw is then simply calculated to be 0.216 degrees using .he arctangent 
function. This of course is an oversimplification of the analysis but still 
gives a rough indication of the orientation of the projectile during the 
launch phase. 


The tube motion data acquired appears to be very good, The BRL loop and 
RARDE sliding mount arrangement worked very well. After each round was fired 
it was a simple matter to slide the sensor loop back into position. As a 
common reference, prior to each shot, the sensor was slid back so the rear of 
the metal loop was 2.0 mm in front of the joint between the «end сар and the 
sensor housing. 


Figures 16 through 20 are examples of tube motion data at the muzzle as 
detected by the BRL sensor. Only the BRL data are presented, because the 
RARDE measurements were made well to the rear of the muzzle and require 
extrapolation of the data to obtain motion at the muzzle. Shown in the 
following examples are gun tube motion with a balanced projectile and with 
projectiles unbalanced and loaded in the gun with the center of gravity (05) 
up, down, left, and right. These records are the actual recorded signals not 
yet converted to actual displacements. 
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>, 
d Shown in Figure 16A are the horizontal and vertical displacements of the | 
tube when firing a balanced round. Also shown is the precise muzzle exit time ห 

of the projectile. Figure 16B is the same measurement with the X and Y scale 
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Figure 16A Figure 16B 


expanded to show more detail at muzzle exit time. Аз can be seen, even though 
the gun has moved before projectiie exit, especialiy in the vertical plane, 
a there is very little movement during the actual time of exit. There are 
noticeable difference when examining the records of an unbalanced round. In 
Figures 17A and 17B it can be seen that the gun tube is moving down at the 
time of projectile exit and abruptly swings up during projectile exit. The 
same effect is noticeable ir. the horizontal plane during exit, the muzzle 
abruptly moves left. Several other motions are noticeable. The projectile 
during its in-bore travel rotates three full turns plus one hundred and forty 
degrees. With the unbalanced projectile there is a modulation of the tube 
displacement for approximately three and one half cycles before exit 
indicating that the unbaianced projectile causes the muzzle to move in unison 
with the spin as the projectile traverses the gun tube. This modulation 
starts at approximately the start of projectile motion although it is hard to 
ascertain from the visual interpretation of these records. Also by comparing 
records of balanced versus unbalanced projectiles, it appears that this same 
modulation is detectable on the summed projectile signal (muzzle exit 
puise). The higher frequency modulation is probably caused by the stress 
waves generated during engraving of the band. Figure 18A and 18B show the 
displacement caused by an unbalanced round but loaded with the 05 to the 
right. Tne same effects are noticeable although in the opposite direction. 
In the vertical plane the gun tube motion going down increases during 
projectile exit while in the horizontal plane the tube abruptly moves right. 
[t appears the effect of the projectile unbalance is forcing the tube into the 
direction of the G; offset of the projectile and then as the projectile 
disengages, the tube abruptly swings the opposite direction. For example, 
with the projectile making approximately three and one half turns during in- 
Gy bore travel and the G; offset before firing to the left, the 05 of the 
projectile would be to the right at muzzie exit. If this forced the gun tube 
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e muzzle right then upon releese, it would swing left. This same effect is 
evident in the tube motion record when the projectile was loaded with the 03 
to the right. Also noticeable in these two examples is how long after 
projectile exit the tube motion continues in that same direction. The same sort 
of effect is evideut in the records shown in Figures 19 and 20, although 
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Тіриге 20A Figure 208 


the major difference in motion is in the vertical plane which seems to confirm "s 
the cause of motion. In these cases the projectile OG was initially EE 
positioned up and down. Ог e again this is aa oversimplification of the 

required analysis but it should provide some insight on how the tube motion is 

related to other factors such as projectile in-bore moticn. Also, with 

further analysis it should be possible to calculate the tube motion at the 

points in time the projectile displacement was measured. This should 

ascertain how much of the measured yaw is in reference to ground and how much 

apparent yaw is due to tube motion. 


Shown in Figures 21, and 22, are composite plots of the horizontal and 
vertical tube motion for a balanced round, an unbalanced round with the 05 to 
the left, aud an unbalanced round with the G; to the right. Approximate 
projectile exit times are shown with "tick" marks since a plot of the muzzle 
exit signal would have made the displacement records hard to read. However, 
the records clearly show the same effect producing tube motion in opposite 
directions. Hopefully, further analysis wili either confirm this effect or 
determine what the actual cause is. Since the tube records are very clean, 
furti. >r details should be obtainable. The records presented here are shown on 
an arbitrary scale, which is sufficient for providing an example of the 
motion, 


A typical record of the observed in-bore yaw, from Round 21, is shown in 
Figure 23. This does not show the behavior of the projectile during the last 
Gx. full turn of the rifling, but it appears by then to have settled down to a 
Yay constant yaw of the maximum allowed by its fit within the bore, turning at the 
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rifling rate. The clearance betweer the projectile and the bore allows a 
maximum in-bcre yaw of the order of 0.2 degrees to 0.25 degrees, The initial 
angular motion, until this yaw angle is reached, occurs before the first 
linear motion (7.5 milliseconds before time zero). The hesitations and 
reversals of the precessiun of the yaw which follow ny be due to variations 
іл the contact force and the frictional force between the shoulder of the 
projectile and the bore; the last occurs shortly after the time of peak 
prersure, 


The result shown in Figure 23 exhibits a considereble difference in rcale 
between the vertical and horizontai directions. This may arise from curvature 
of the supposedly plune mirrors. One of these, Sitvated necessarily on the 
bore axis, is destroyed by the shot, it cannot be supported anywhere but at 
its edges. 
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(ARROW LENGTH REFAC. XNTS 0.0858 DEG.) 
Figure 23. 
7. SUMMARY 


At this time all preliminary avalysis results indicate that the joint 
test was very successful; but, how successful really is never knowr until the 
data analysis is completed. The quality of the dala appsars to be very good 
and, with few excepticas, 511 the data desired are recorded. No problems were 
encountered in interfacing the BRL and RARDE systems cr periorming the test. 
in fact, the only problen in the whole program was a delay caused by the 
exchange of drawings which were temporarily sidetracked in the official US 
comunicatior chaunel, 


This was the first time спас tube motion measurements were made using the 
Sch, dr displacement transducers and a great deal cf informatlon wae obtained 
which chould be useful in future test. Based on th apparent high quality uf 
the tube mofion daca, future tube motion measurements ucing the Schmidt 
aisplaccrent transducer will employ at least (wc sensors tc obtain the angle 
of tube dis)Jacement at muzzle exit. The recorded signals of the projectile 
displacements are ล โฮ ง of good quality and it is hoped that compariscn with 
rhe in-Lore laser meagureneats made by КАКОЕ will further confirm the BRL 
weusurereatu. Future tests will be conducted using a symmetricel mounting 
collar with the projectile measurement electronics ізсзбей oft of tre gun. 
This will provide a symmetrical measuring surface for dial tube mot lor 
measurements and protect the electronics. 
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TITLE: Resulte of the Joint BRL-RARDE 40-mm Firing Experiment to Define 
Projectile Launch 
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USA Ballistic Research Laboratory 
Aberdeen Proving Ground, MD 21005-5006 


Dr. Roger K. Fancett 
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ABSTRACT: 


Using the BRL and the two HEL displacement data sets individually as 
well as combined, the muzzle motion of the instrumented 40-mm L70 cannon is 
determined and evaluated for the totai interior ballistic cycle from the 
initiation of the firing until shortl: after the projectile has left the 
muzzle. The validity of the indirect y calculated muzzle recoil is examined 
by comparing it with recoil даса obtained by an earlier experiment with the 
ваше gua and ammunition. The time window corresponding to tbe projectile-gun 
disengagement is expanded to the same temporal resolution as used in the 
analysis of the projectile motion, and the linear displacement vector and the 
pointing direction of the coordinate system which moves with the muzzle are 
computed as functions of time. 


The projectiie motion - the vectors tor the linear displacement of 
the center of mass and the angular orientation of the axis and their 
derivatives with respect to time — during the projectile gun disengagement 
time is computed by employing the previously described data analysis program 
to the two recorded data sets, graphically presented with respect to both the 
muzzle and the ground frame of reference, and discussed. The axial component 
of the projectile velocity obtained from the analysis of the Schmidt 
displacement transducers is compared with the velocity recorded downrange and 
the velocity gain in the muzzle blast regime is estimated. The projectile 
yaw, i.e., the orientation of the vector for the projectile axis with respect 
to the velocity vector of the center of the mass, is displayed with the data 
from Fuller's optical lever arrangement overlaid, whenever applicable. 


The variations in the initial condition - projectile unbalance and 
its orientation in the gun chamber — are correlated to the linear and angular 
motion of the muzzle and the projectile immediately before and after the 
projectile separation from the gun barrel and related to the recorded impact 
location on the target. 
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This paper was unavailable at the time the Proceedings were published. For 
further information, contact the authors. 
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TITLE: Data Analysis Procedure for the Schmidt Displacement Transducer 
to Extract Projectile Launch and Muzzle Motion 


Dr. Rurik K. Loder and Emma M. Wineholt 
USA Ballistic Research Laboratory 
Aberdeen Proving Ground, MD 21005-5006 


ABSTRACT: 


This paper delineates the analysis concept and procedure used in the 
development of the data analysis program. The projectile-related analysis 
program contalns: (í) the data base describing the frame of reference, the 
bore geometry of the muzzle in reference to the locations of the two sensors 
used in the experimeut, the contour geometry of the projectile, and the 
recorded data; (1i) the analysis of a subset of data to establish points of 
reference for the temporal alignment of the projectile contour with the data 
sequence and to extract the axial component of the projectile motion; (111) 
the analysis of the remaining data to spatially align the coordinate systems 
of the two displacement sensors, extract the transverse displacement of the 
center of the projectile contour from the sensor axis from the data in the 
presence of large amplitude colored noise, and convert these displacements 
into the vectorial description of the motion of the projectile's center of 
mass and the angular motion of projectile's axis about the direction given by 
the velocity vector of the center of mass; and (iv) the utilization of the 
results of the auteceding analysis to determine the pivotal location on the 
projectile axis and tne orientation of the tube axis with respect to the frame 
of reference, represent the projectile motion with respect to the tube axis, 
and estimate the asymmetry of rotating band engraving. The frane of reference 
is the coordinate system given by one of the two sensor systems and, hence, 
moves with the muzzle. The analysis program allows the superposition of the 
muzzle motion and the description of the projectile motion with the respect to 
the frame of reference in which the muzzle motion is defined. 


The analysis program for the muzzle displacement data is straightforward. 
The data are converted into displacements and superimposed on the displace- 
ment equaticn for a moving, flexural beam in the least squares sense. 
Essentially, the transverse displacements of the muzzle section of the gun 
barrel are set forth as polynomials in the axial position, with the 
coefficients varying in time and satisfying the end conditions for a free beam 
at the muzzle face. Since the sensor position with respect to the muzzle is a 
function of gun recoil, the correlations between the data and the model 
parameters are expressed such that five of the six degrees of freedom are 
obtained. Only the rotation about the bore axis is not extractable. The 
program permits the analyis of data sets from one or more measurement 
stations. If the data of only one station are used, the muzzle section is 
considered a rigid base which moves in space. 
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TITLE: Results of the Joint BRL-RARDE 40-mm Firing Experiment to Define 
X Projectile Launch 
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ABSTRACT: 


Results from the 37-mm and 40-mm caliber firings are used to 
evaluate the performance of the Schmidt displacement transducer for the 
measurement of the muzzle motion as well as the projectile dynamics during 
the final phase of the interior ballistic cycle. Limitations and shortcomings 
in the current instrumentation design and measurement technique which includes 
the data analysis procedure are díscussed, together with possible avenues of 
its improvement and research ín progress. 





This paper was unavailable at the time the Proceedings were published. For 
further information, contact the authors. 
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ABSTRACT: 


This work presents new measurements of transverse motion of the 30 

ma (GAU-8) gun tube first reported at the Third Gun Dynamics Symposium in 
1932. The measurements are unique and fully corroborated through the use of ` 
two Independent measuring devices. In particular, the work discusses three n 
items of interest. First, there definitely exists a 'base-line' transverse | 
tube movement, the cause of which has yet to be determined. The magnitude of 
this motion is of the same order as that due to other causes intentionally 
introduced for study. Second, intentionally int:roducing an eccentric breech 
mage produces muzzle digniacements in good agreement with theoretical models, 
provided the base-line component is accounted for. Finally, the muzzle 

| rotation created by the moving projectile - though insignificant when 

& operating alone ~ is strongly coupled to, and capable of greatly modifying the м 
rotation due to other causes. This coupling does not appear to strongly 4 
affect muzzle displacement. Іс is concluded that predictions from gun » 
dynamics models which agree well with displacement measurements, may err 
greatly when used to predict muzzle rotation, tbe quantity of which is of 
central interest. 
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1. INTRODUCTION 


The work herein reports further measureme: :8 of the transverse 
muzzle motion of ап elastically suspended 30 mm gun tube. When last reported 
[1], there were significant but unexplained muzzle movements prior to shot 
ejection. As a result, the credibility of the measurements of these motions 
was brought into question and it was decided that further analysis could not 
be justified until fully corroborated measurements were in hand. It was 
decided that two completely independent methods of measuring muzzle 
displacement would be employed simultaneously and that close agreement of 
these measurements would be demanded for acceptability. This quest for 
corroborated - and hence believable - measurements was highly successful аз 
the results herein will demonstrate. 


Having the proper instrumentation and measurement techniques in 
hand, however, did not guarantee immediate success ia totally explaining the 
motions of the 30 mm tube during firing - the central problem of gun dynamics 
to which we have finally been able to return. Latest measurements, although 
greatly improved from those reported in 1982 [1], still reveal an unexplained 
traasverse movement of the tube prior to shot ejection. Furthermore, the 
magnitude of this motion is roughly the same as that which is to be studied 
from intentional causes such as tube curvature, tube and/or projectíle 
ecceutricity, etc. There can be little doubt that the cause of this 
underlying motion will be found, but as yet this is not the case. 


2. EXPERIMENTAL DESIGN AND INSTRUMENTATION 


Referring to Figure 1, the 30 mm/GAU~8 gun tube was suspended as 
described in Reference 1. Briefly, the tube is suspended by two pairs of 
wires from a virtually rigid overhead structure. The wires are 6.026 inch in 
diameter and each ts approximately 50 inches in length. The recoiling tube 
thus behaves as an elastic bar-pendulum, stretching each wire approximately G 
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0.010 inch at shot ejection. A NASTRAN model of this means of support e 
indicates that the wire support loads create negligible transverse tube motion 

prior to shot ejection, i.e., that the suspension as modeled caa be regarded 

as 'soft.' Following shot ejection, the recoil motion is arrested by contact 

between a circular buffer plate attached to the breech and a block of 

open-cell silica foam. (This foam has been found far superior to the 

'styrafoam! used previously.) 


The propellant charge used in these experiments is half that 
normally used ín a standard GAU-8 round and yields a maximum pressure of about 
10,000 psi. The round із separated from its case and manually started into 
the origin of rifling with a light tap. The cartridge case containing the w 
propellant is then inserted and forced home as the threaded breech cap is : 
tightened. The breech cap contains a small central recess into which ап ; 
electrically actuated, propellant-driven Firing pin is inserted.* 


Instrumentation used to pertorm the measurements presented herein 
consists of two late model optical trackers (trade name - 'Optron'), and one 
so-called ‘eddy probe' manufactured by Dymac Division of Scientific-Atlanta. 
One Optron is used with a light source which is interrupted as the shot leaves 
the muzzle. This, plus a breakwire mounted across the muzzle provides two 
independent measurements of the time of shot ejection. Agreement between the 
two is within 0.0002 seconds. It is hoped that in the future an inductance 
device will improve this measurement. The remaining Optron is used to record 
the vertical motion of the tube approximately six inches rearward of the | 
muzzle. From this point the muzzle protrudes into a tube which exits through © | 
а wall of the Coon, The tube protects the instrumentation from muzzle flash 
and smoke. The end of this tube also serves as a mounting place for the eddy 
probe which is positioned 0.050 inch directly over the point of the tube being 
followed by the tracker. Thus the tracker and the eddy probe both follow the 
motion of the same point on the upper surface of the tube at ali times. Owing 
to recoil, of course, different material surface points are being tracked in 
time, l.e., the measurement frame is Eulerian, not Lagrangian. This 
difference Ls of little consequence for the purpose at hand, however. 


Measurements are monitored and recorded digitally on а four-channel 
Nicolet system. Approximately twenty milliseconds of data are recorded from 
each round using the muzzle breakwire and a trigger. 


3. MEASUREMENTS OF GUN TUBE MOTION 


Figure 2 shows the moticn of the muzzle as measured by an optical 
tracker and by the eddy probe. Shot ejection is taken at t = 0.0. It is 
estimated that most of the time agreement із within 0.0002 inch (0.2 mils). 
Assuming sufficient care Ls taken to duplicate all of the preliminaries, a 
second shot will produce measurements in close agreement with those of the 
first shot. This is exemplified in Figure 3. 


*These 'mini-actuators' are manufactured by ICI Americas, Inc., Atlas 
Aerospace Division, P.O. Box 819, Valley Forge, PA. 9 
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Presently there is no explanation as to the cause of the motions 
depicted in Figures 2 and 3, i.e., no eccentric masses have been applied to 
the tube, and the curvature of the tube "ав manufactured’ is not sufficient 
(according to our two computer models) to create motion of this magnitude. 
Likewise, the support reactions applied by the wires are not supposed to 
change radically during recoil. (A thorough check of these reaction forces is 
presently underway.) 


Figure 4 shows the displacement of the muzzle when an eccentric mass 
is added to the breech eud cf the tubo. The eccentric mass is created by 
simply replacing the 7.75 lb circular buffer plate with an identical one whose 
center is located below the central axis of the tube. The weight of the 
entire recoiling mags is 90.8 lbs. The eccentric location of the buffer plate 
is 1.3 inches directly below the bore axis of the tube. In Figure 5 the 
response shown in Figures 2 and 3 is taken as a 'base-line' and subtracted 
from that of figure 4. The result is in reasonable agreement with that 
predicted by a NASTRAN model of the system. This is shown in Figure 6. 


4. MODEL VALIDATION 


The agreement shown in Figure 6 Ls at least as close as that in 
Figure 2, if one disregards the higher frequency components of the motion 
(probably radial vibrations); i.e., the experimental measurements and the 


a ал ft^ бил nama จ 4> ¿x A โด a ль 
NASTRAN=sredtctcd müzzte displacements agree to the same extent as our 


measurement devices. It would appear, therefore, that we have a useful and 
validated computer model. Assuming this to be the case, one of the most 
expected uses of a computer model of a gun system is to predict muzzle 
rotation. It 1з of much more interest to kaow muzzle rotation than it ts to 
know muzzle displacement, target error being associated most strongly with the 
former. Can our NASTRAN model be crusted to predict muzzle rotatíon? The 
answer is very negative. One reason for this will be shown ~ the fact that 
our NASTRAN model does not account for moving projectile mass. Unfortunately, 
there may be others. 


In Figures 7a through 7d the predicted muzzle response due to an 
eccentric breech mass with and without the presence: of the moving projectile 
mass are compared. (These predictions were produced by a second computer code 
іп which moving mass effects are easier to include.) It is observed that 
although there is close agreement between muzzle displacement predictions with 
and without woving mass effects, the prediction of the (all-important) muzzle 
rotations aie totally different. Thus the effect of moving projectile mass 
cannot be ignored. 


While lt has been shown that moving projectile mass must be included 
in any mathematical model from which muzzle rotation predictions are expected, 
this is by no means a sufficient condition. Other effects normally overlooked 
and expected to be negligible, may have to be included also. Most important, 
one may not know what they are, or in what detail to describe them. For 
example, the muzzle rotation predicted when th» moving projectile is included 
may well depend on the detail with which it is modeled. (In Figures 7c and 
7d, the projecttle was simply modeled as a moving point mass.) It may be 
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concludad, therefore, that the only way to be certain that а по1г1 can be 
trusted to predict muzzle rotations is by experimental verification ~ the very 
tesk one hoped to avoid! In conclusion, complete motal validation requires 
naasure aent of displacement and rotacions. 


That the inelvston of movin, n"rojectile mass in any msthewatica? 
modei of a gun system wonld have such а draratic effec. on predicting muzzle 
rotations ณะ ท como as a suvpris: to sew investigators. After all, 
predictions of maczie motion due to a moving projectiles wass іл the азсепсе of 
ether load functioso (such as eccentric breech mass, support reactions, 
curvature, etc.) hav: been virtually vegii, ible. One mar not forget, 
howeve.', thet the tube respouse due tco combinations of lead functions which 
are mozion-Cependent, ie not simply the. sum ot the respuuzes due to eacl. load 
acting indepentently. Even though one із dealing in generel with z Last 
model based on soue particular Linear partial differentia? cquation (pec c.), 
the differential operatur is altered every tim. a motion-depindent 1034 із 
incluósd or excluded. The operator must remain the same in order for Linear 
Superposition tc hold. For exemple, consider the partial cifiereutsial 
aquations correspundiag to Figur з Ti. through 74. If tue movin, proldectiie 
mass and an eccentric breech mast are both ізсінісі, the р.д.с. is 


๑ ๓ ๑ 


(ETy")" + py = «Мщ(с)б'(х) = nig (vty! bei y" eg) fine C) (1) 


where aft} is the recoil ace:leraticn of a breech mes M located distance € 
from the bar? axis (see Figure 7a), m, ie the masc of the projectiie located 
ai a Cintauce (ty from the breech end ог Lie cube, v i: che trzangverse 
disp.acemeni of tue tuba, 6 із the Dicuc fenction, and g is tha grevitaticnal 
constat. Note that the projectile. mass creates а norion-dependent load 
funclio: = it depands оп y and it derivatives, Tie solution to Eq. (1) at 
the muzzle is whawn ya Віри: es 7a and 7b. 


Эа the other hand, Lf che moving projectile is neglectcá 


(EIy")" + ру = :Ma(t) (x) (2) 


aad the solution is shes . ia rigures 7c «па d. A quick glance shere that the 
displacement is hardly changes from Figure га, but the rotation (J'Y is 
coupictely different. 


Now соатійег the cese where the moving ร ย 43 ๑60316 is the only 
«Сай acting). The o.d.e. js 


m DI 


CEL,")" + ny = ap(y+2y'+E2y"+g) 6(%=6) (3) 


Yue muzzle displacement and rotation corresponding to this equation are shown 
in Figure 8. Ав can be seen, their magnitudes are very small. ТЕ these are 
added to those from Eq. (2), hardly any change to Figures 7c and 7d would 
resuit, demonstrating that solutions to Eqs. (2) and (3) cannot be summed to 
yield the solution to Eq. (1). [his is shown in Figure 9. Іп effect, the 
motion-dependent term affects the operator of the p.d.e. and one cannot add 
solutions to equations with diiferent operators. 
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5. CONCLUSION e 


l. The measurement of transverse tube displacements during firing 
can now be accomplished with a high degree of confidence. Eddy probes and 
optical trackers may be expected to disagree as much as 0.0005 inch (0.0.27 
mm). Neither device can be assumed superior to the other, but eddy probes are 
relatively inexpensive - almost expendable, and very easy to use by 
comparison. 


2. There Ls no way to know if all motion-dependent load functions 
have been accounted for in a given mathematical model or if those which have 
been accounted for have been described in sufficient detail. It has been 
shown herein that motion-dependent load functions, which are unimportant when 
acting alone, become very important when acting in concert with other loads. 
The importance of such motion-dependent loads cannot be assessed frow their 
effects on tube displacement which may be minimal. Tube rotation, however, 
may depend significantly on the inclusion and proper description of these 
loads. Thus, a model should not be considered validated until displacement 
and rotation have been verified experimentally. Weedless to say, all models 
must be validated to be of any use. In view of this, it is important to the 
future of analytical gun dynamics that means be found to obtain measurements 
of tube rotations - particularly at the muzzle. In addition, such 


measuremcats should be corroborated. 
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INTRODUCTION 


Whenever a gun is fired a projectile is propelled violently along the 
barrel, whilst at the same time the barrel recoils, The magnitude of the 
resultant accelerations can cause large forces and moments to be applied 
to the barrel in the lateral direction and these then produce flexural 
waves that travel along the barrel. As a result, by the time the projectile 
is launched, the muzzle may have considerable transverse velocity and 
displacement and thus the projectile departs in a slightly different 
direction to that which was intended. 


A second effect is that whilst the shot travels up the barrel it 
may undergo considerable vibration relative to the barrel due to the finite 
stiffness of its driving bands, gyroscopic motions, shot/barrel clearances 
etc. Therefore even in a perfectly rigid barrel the shot may depart with 
considerable lateral disturbance due to its internal motion within the 
barrel. 


In practice both of these effects are present at the same time and 
they are probably coupled together as well. The calculation of these 
effects poses some considerable computational difficulties. The 
vibrational frequencies of the shot whilst in the barrel are quite high, 
usually in excess of 1KHz, but the flexural motion of the barrel 
usually has a fundamental frequency of only a few Hz. As a result using 
standard numerical integration techniques, the equations of both the 
shot and the barrel would have to be solved using very small time steps, 
sufficiently small to ensure that the high frequencies of the projectile 
dynamics were handied in aa accurate and stable manner. Typically step 
sizes of the order of 10u secs would be required to compute the shot motion 
and on a computer such as the VAX the solution of the motion of the shot 
alone would take quite a considerable time. Add to this the model of 
the gun barrel with 100 or mere degrees of freedom and the solution time 
would escalate considerably. 
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This report addresses the problem of how this excessive solution 
time can be reduced. In addition it attempts to use the already established 
SHOCK-AID solution techniques Ref. (2) in such a manner that advantage is 
taken of the ability of SHOCK-AID to model barrels simply and efficiently. 
The major problem however with the latter approach is that SHOCK-AID is 
based upon e modal su,erposition system that takes advantage of the linear 


nature of the flexural equations of the barrel. Motion of the shot is 
non linear and hence does not lend itself to modal style solution, 


and as a result some method has to be found to couple the modal solution 

of the barrel flexure with the numerical integraiton of the shot motion. 
Details of the shot model were given in Ref. (3) along with very brief 
details of the interaction technique. This paper considers the computational 
details cf the shot/barrel interaction in detail. In order to demonstrate 
the nature of the problem and its solution, a simplified moda! of shot/barre! 
interaction is considered so as to remove non essential detail from the 
problem. 


2.9 Ап initial simpiified modei 


The shot is a relatively small mass connected to the barrel via 
stiff springs, whe, sas the barre] is a very large mass having a relatively э, 
low stiffness. As a result the simplest model of the two elements 
combined together is a two mass spring system as shown below 


Ig к, 
1 " Ne ipe 


Xg —— xç‏ مو 
ms and ks represent the mass and stiffness of the shot whilst Mg and‏ 
kg represent the mass and stiffness of the gun. External forces applied‏ 
to the shot (they may include dynamic effects such as gyroscopic forces‏ 
etc.) are represented by F.‏ 


The equations of motion of the system are 
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assuming х. (0) = х (0) = ks(0) = Ха(0) = 


and using Laplace Transforms we can write, 





(ко + k )/nç | + k.)/m 
eee ๕ Б -1 + ыл COSWit 
n м: (wî - wi) 


(k. + k 
" * ARIES соѕм, t (1) 
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and w and w, are given by 


k. (ky 





k 
- 4-2 si (3) 
9 


Now in order that we have representative values for the constants 
in the above we assume the following: 


К; = 3.95 E10 N/m 
л = 10 kg 

kç = 6.32 E8 N/m 
Mm, = 100 kg 

F = 6.32 Eb М 


These figures mean that the natural frequency of the shot mass on its 
own is about 10 KHz whilst that of the barrel element on its own is 


about 400 Hz. The load F produces a static defiection of the barrel s 
element of Im. Kai 
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Using the above data in equation (3) gives 


% w 


1 


k 
; z 
> wš п (4) 


nt 
“е” 
ә 


and since wŠ >> wi we сап write 


cs RS 
o? ға cosw,t) (5) 
g 
Ë k 
Xç š l = coswit - d cosw,t) (6) 
and thus 
x -x sË coswet (7) 
s g к 2 


In order to represent the shot barrel model using SHOCK-AID the 
obvious approach would be to solve the motion of the gun mass using 
the SHOCK-AID modal technique and solve the shot motion via a numerical 
solution such as Runge Kutta or a predictor/corrector technique. In 
order to link the two solutions some form of iterative scheme would 


@ be required. 


If we denote the force in the shot spring by the letter f,i.e. 
f = k (x; - ха) (8) 
then the equations of motion become, 


T ыы (9) 


m. น = Е - f (10) 
If we consider the barrel mass and spring alone we can write, 


k 
f f 

Kose 9 х ET + — (11) 
т т 9 m 

and in order to solve (11) we need to know the variation of "f" with 
time,however, as equation (8) shows "f" is itself dependant. upon the 


solution of (11). Do 
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Ignoring for the moment the implicit variation of "f" with "x 


we can write the solution of (11) as, 


t | 
f(t) sin w(t ~ 1)d 
xj 0) + Í Cou Sin M corde 


xX = 
g w 
0 9 
t 
A = v (0) + f f(t) cos w(t - т)дт 
g و‎ м, 
o 


Let b = f/m, then expanding the integrals thus, 


t t t 
( b(1)sin w(t - t)dt _ ОША [ blx)cas (ลา mE í 


w 
0 0 9 


and the integrals can be represented by 


c2 


t+At t+At 
| Ь(т)ѕіп wrdt = k(t + At) = k(t) | b(r)sin мтат 
0 t 


A linear approximation to "b" is used by ietting b = mt + c over the 


step At, then, 


= b(tg + At) - b(tg) 
ШЕ 1 1 


o 
H 


= b(t) - nt, 


We can thus write (15) and (16) as 


ü 


J(t + At) = 0060) + JINC(to) 


k(t) + At) 


" 


k(t ) + KINC(t,) 
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б мһеге 
ds mtotAt 
to*At А ; 0 
JINC(t, ) =f (тт + c)cos wrdt = mcos WT + MSIN WT , <510 мт (21) 
t м? м w 
o ty 
0+4 | dio 
KINC(t ) -f (m + c)sin widr = [min wr _ mcos wr _ csin wr 
0 2 (22) 
w w w 
to to 
and thus the final sclution for Ke and Ve is, 
2 sin wt Š 
ม อ (ไอ + At) = xg (t) + TE (to) + JINC(t} 
cos wt г x i 
- SOSIE (to) + KINC(to)} (23) E 
f ent 
v (to + 44) = vg (to) * COS wt qo) + KINC( ty D N 
@ + sin wt $c " лю) (24) m 
The above equations represent a SHOCK-AID moda! type solution to the 
motion of the gun mass. 
Having produced г SHOCK-AID style solution for the gun element 
and a standard numerical solution for the snot, the next question is 
how to integrate the two calculations together. A possible scheme 
is as follows: 
(i) Estimate the acceleration of the gun element uver the interval 
ty > to + AT 
% = (коху - f)/m, (26) 
(ii) Compute the shot motion using the predicted motion cf x 
т у 2 
ха (1) = xg (to? + v (to)at + $v (tp)at (27) 
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to calculate the variation of "f" over the step to > to * AT. 
At the same time compute "avsf" the average value of f over 
the interval. 


(iii) Compute the motion of the gun element using f - avsf as the 
applied spring force over the interval. 


(iv) Compare xg(to + AT) from step 3 above with xg(t + АТ). If they 
agree (within a suitable tolerance) then advance the solution 
otherwise re-estimate Vg at step 1 and repeat the calculations. 


A program to implement the above was produced, It used the average 
force over the interval not a linear variation. Ол?“ representations 
of the force variation were investigated, but all versions of th program 
appeared to have a common dràw back in that it was not in general possible 


to drive Xg and xó into agreement. No matter how long the program 
iterated it always Seemed to converge оп a solution Which had а smali aan 
but steady error between Xg and хо. 3 usn 


The problem was that the modelwes too simple in that 
it only uses two mass spring elements. A more representative model 
is a three element model. The reasons for this are explained in the 
next section. EE 
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3.0 Three Mass Spring System 


As mentioned above a twc mass/spring model for shot barrel interactions 
is too simple. The reason for this is that there are in fact THREE 
types of elements in such a problem. They are: 


(i) The shot element or elements 
(ii) The contact elements, i.e. those elements that the shot applies 
forces to directly. 
(iii) Internal gun elements that have no direct contact with the shot. 


The reason for treating each of the above elements differently 
is that they experience forces of very different frequency content, A 
shot element is subjected to very high frequencies due to its low mass 
and its high contact stiffnesses, an internal element on the other 


hand is only exposed to relatively low frequencies associated with 
the barrel dynamics. The contact element however is exposed to both 
the high and the low frequencies. As a result our three mass/spring 
model of barrel/shot interaction is as shown below. 


displacements 





= D —g — f spring forces 


The two solutions (modal and numerical) now overlap in that the 
modal solution computes the motior of all the internal elements and 
all the contact elements, whilst the numerical solution computes the 
moticn of all the contact elements and all the shot elements. In order 
to do this some assumptions are required about the behaviour of elements 
excluded from each calculation. 


A possible solution scheme is as follows: 


(i) Predict the motion of mg over the next time step AT 
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* ๑ (1) = xglty) + Kelty at + Zen 25 (8) (7) 


to < t < % + AT 


XGay is a guess of the average acceleration over the interval AT, initially 
set this to Xe(t.). 


(ii) Solve the equations of motion of masses m, and m, using the 
assumed xg(t) above so that spring force g can be calculated. This 
will use small time steps 5T where in general 6T << AT. During this 
calculation monitor the variation of the spring force f and calculate 
f у» an "equivalent" spring force that can be applied to the modal 


eq 
solution. 


(iii) Apply ци to the modal solution of the contact elements 
and the internal elements. 


(iv) Compare seit + AT) from (iii) with the predicted value from 
(i), if they agree within a reasonable tolerance then advance the solution 
another AT and repeat trom (1); otherwise re-estimate Ход, used in а=, 
(i) and repeat the calculation, C 


In order to carry out the above, two items are requirec, these 


a + an algorithm for calculating “Т " and another for updating Gay: 


qv 
Considering the estimation of XG first. The initial estimate 

of Xeay is Xg(t,) since no data about the future variation of XG is 

available. On the next iteration however we have available Xg( to * AT) 

that was calculated by the modal solution and thu a reasonable estimate 

for Xg,, would be that value that would produce a predicted value of 

Xa tg + AT) equal to that observed. Thus using (28) we can say, 


or 
XGAV = a(x. ét + AT) - xg (ty) - Ke (tat) /at? (30) 


The calculation of "Ten is a little more complicated. The standard 
SHOCK-AID modal solution assumes that the variation of applied force EI 
over the time step is linear. Thus the calculation of "ไอ ย ุ ห must еҙ 
produce a linear force variation over the time step AT that is equivalent 
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to the (in general non-linear) variation of f over AT (f is available 
at intervals of ST). 


The criteria of "equivalent" in this cas. is that the displacement 
and velocity of the large mass should be the same at the end of the 
long time step irrespective of whether the linear or non linear force 
is applied to it. This also means that both the potential and kinetic 
energies are the same for the "real" and the "equivalent" forces. In 
order that the two forces produce the same defiections and velocities 
at the end of the AT time step, examination of the SHOCK-AID equations 
shows that the values of the following two integrals must not change 
if either the "equivalent" or real forces are used. 


Tu A tat 
f(t)cos wtdt ; | f(t)sin wtdt (31) 
% % 


for all w'5 (w's are the gun system eigenvalues) 


Now let T, = En qun system period(s) 


If т > AT then we can write 


.tg*AT ЗАТ 
2 Ant 
#(t )cos~— dt = f(t)dt (32) 
t е | 
о o 
tatl tas AT 
° ; ¿Ont ui un 9 2nt 
f(t)sin?.- dt z f(t) 7 dt 
S w t d 
o e 
tytAT 
= F £(t)tdt (33) 
“1 
0 


and the integral іп (33) сап be written, 
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t +éT 


f f(t)tdt 


% 


t +aT š t +AT 
[ td |: #(т)ат Mt 
to to 


. +AT T+AT -t 


uf (ов - f f cona 
to % 


% 


" 


Thus the "real" and the "equivalent" forces must have the same values 
for the first and second inteqrals, 
t rat $ +AT t 


f F(t)dt and С [сова 


% % % 


If we call these two integrals for the "real" force 1+ and Ilẹ respectively, 
thes we require for the linear equivalent force, 


t +AT 
° maT? 
(mt + c)dt = + с АТ = If (36) 
t 
0 
t tat 
2 42 c 
[^ eÉ cta = м, Ат „ү, (37) 
t 
0 


where cy (mt. + c) 


These car, be written 


; 5 | 
Ir] _ }АТ°/2 М т 
PIER ว | | (38) 
í ie А1572) Lead 


the determinant of (38) is Bit? and the inverse is, 


12 | 12/2 Ar 
15 [273/6 дт2/2 
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d thus 
auch ы TE qp (35) 
AT? AT? 
2 6 
Желле, (40) 
CT E cb 


Direct implementation of equations (39) and (40) would mean that 
the equivalent force could be discontinuous, e.g. 


fw 





What is more in keeping with the original form of SHOCK-AID is 
| for "со" to be determined by the previous time step. Then in that 
a case, if tne first integral is correct 


m = (If - CoAT) (41) 


AT? 
and the value of the force at the end of the step is given by 


_ 2 
maT + C. = т (iç - CAT) + Cy 


= 514 - Co (42) 


Alternatively if the second integral is correct 


6 CAAT? 
m= — (Ilf - äs (43) 
AT? 
and 
6 à 
ПАТ + c, = — II, - 2c (44) 
0 AT 0 


HG 
n 
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and thus it is not possible to satisfy both equations (42) and (44) 

at the same time if "co" is given by the previous step. 9 
fhe question therefore arises, if we can only satisfy one of the integrals 

which one is best? 


If we satisfy the first integral then the cosine terms in the 
solution are correctly handled (see equation (32)) but there is an error 
in the sine expressions. From equations (33) and (34) the error in just 
using the first integral gives an error in the sine integral of, 


Фо+АТ t 
maT? AT2 
є = т, | J f(r) та Е =e = Com 
o o 


where m is given by equation (42) for a c prescribed by the previous step. 


Thus 
an f АТ oiT? 
, с 
є = T. ff - KR Ig br ) (45) 
It is not easy to tell from the above how the alternatives behave @ 


and a more convenient test is to assume that 
f = sin T (46) 


and observe how the different techniques approximate the function. 


First of ail using equations (46) we can write for a time step AT, 
Ig = cos(T) - cos(T + АТ) (47) 


Ilẹ = sin(T) + ATcos(T) ~ sin(T + AT) (48) 


Using equations (39) and (40) in conjunction with (47) and (48) 
gives the results shown in Fig. 1. The curves show that 
for both large and small values of AT the approximation to f is well 
behaved, and the discontinuities are required in order to match the 
original function. 


Using equations (41) and (42) with (47) and (48) removes the 
discontinuities but only the ist integral is satisfied. Curves of the ว 9 


6 
Е - 
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results are shown in Fig. 2 for a range of AT values and 
the approximation is well behaved for all values of aT, there being larger 
discrepancies at large AT as would be expected. 


If equations (43) and (44) are then used so as to satisfy the 
second integral only, the results shown in Fig. 3 are obtained. 
In this case the discrepancies are very large and somewhat surprisingly 
increase as AT decreases, there being some form of unstable behaviour 
apparent. 


That this is indeed unstable behaviour can be illustrated as 
follows. | 


For the Ist integral case we can write equation (42) as 


_ eI 
014+1 - E <. Ci (49) 


and if AT << 1 then Ig = aTsinT = aTsin(iaT) 


91+1 = @sin(iAT) - cj 


and taking z transforms gives 


z(c(z) z co) Е 2zsinaT E c(z) 
2? = 2zcosaT + 1 
or | 
Blah se ๐ 02 , 1 2251іпАТ (50) 


(2+1) {z+ 1) (22 - 2zcosaT + 1) 
the poles of (50) are given by 
2 = -1 and z = COSAT + i sinaT 


i.e. on the unit circle and thus the solution is stable. 


For the second integral case equation (42) can be written, 


= SIF - әс; (51) 


AT? 


Cist 


and if AT << 1 then 
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Ig _ sin(iat) 


T2 2 en) 
4 - 
and thus 99 
C; 4 = 3sin(ist) - 261 (53) 
and taking z transforms 
z(c(z) Е Cy) = _— Ə3zsinT  — - 2c(z) 
22 - 2zcosaT + 1 
or | 
๐ (2) = EH + 1 __ 3zsinaT (54) 


(2 +2) (z+ 2) (22 ~ 2zcosAT + 1) 


z=-2 and 2 = cosaT + i sinaT 


i.e. outside the unit circle and therefore unstable. 


Thus the conclusions are that ideally both integrals should be 

satisfied and the results will be valid for all step sizes, or for small 

steps the ist iniegral alone сап be satisfied. Satisfying the second 

integral alone results in unacceptable behaviour for ail step sizes. © 
The preceding sections provide a basis for the solution of the three 

mass problem. The development is split into three phases as for the two 

mass problem. Firstly a modal solution is developed for the internal 

mass and the contact mass, making suitable assumptions regarding the motion 

of the shot mass. Secondly a numerical solution is derived for the motion 

of the shot mass and the contact mass, making suitable assumptions 

regarding the motion of the internal mass. Finally, the two are integrated 

together using the ideas of the previous paragraphs for Хбдү and fa 


qv" 
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3.1 Modal Solution of Gun Masses 


The equations of motion of the three mass model are, 


mç X6 = -6 + 9 
= = f 
mç xç 9 + 
ША X. = =f + Е 
where, 
g = ко xg) 
f= к.(х- - xg) 


The equations of motion of the gun masses alone are, 


0 m x kg › “ky D f 
or x = M kx + kf 
and 
м2 = [Um 0 
о 1% 
Мк = |-(kG + kg) kg 
пс mG 
kg . Kg 
Mg mg 
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(55) 


(56) 


(57) 


(58) 


(59) 


(60) 


(61) 


(62) 


(63) 


(64) 
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the eignvalues of (64) are g ven by 


momo 


E PIDE 


The standard SHOCK-AID notation is 


Lk 
wend кб + kg $ d т “akg = 0 


+ Ax = b(t) 
thus 
А = Mk 
b= Mf 


The left hand eigenvectors of A are 


г = 
he 
г 
๐ 
eo 
[=] 
- 
— 


where 


1 I 1 zi 
Y= ซอ ธี |. та 
‘еб: Жеб, 
where 
Det = еб: - 4:6, 


Let а; = Yy-W, and а; = И-М, 


The final solution is 
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(65) 


(66) 


(67) 


(68) 


(69) 


(70) 


(71) 


(72) 
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t 
x(t) = U, = Í b(:)sina,(t - t)dr> V; 


t 
td | Ь (т) 59 паг (& - +)dzy V, (73) 
= 2 

0 
Equation (73) can be manipulated in a similar manner to equations (12) 
and (13) so as to produce marching integrals 


3.2 Numerical Solution of Shot Motion 


The equations of motion for the shot and the contact mass, assuming 
that mg is motionless are, 


g = k (x, - xc) (74) 
f = k(x, - xg) (75) 
จู =i (-f EI T 
s 
@ A (ว 7) 
เง — (-9 + f) (78) 
q 
55 = Vg (79) 


A numerical solution to themotion of these two masses is used. 


3.3  Coupled Solution 


The two programs above were then integrated into one program using 
the interaction scheme outlined previously to produce a single program, 


Basic data used in these programs is given below. 
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Me = 100, т = 100, m. = 10 5) 
ke = 6.32 ER, ky = 6.32 E8, k; = 3.95 EIC 
Е = 6.32 E5, t%?0;F=0,t<C 


The program has two forms, 


(1) that uses first integral only 
representation of feav' 

(ii) that uses both first and second 
integrals to represent feqv' 
The ห ล ร น ไข ร for the Ist integral case are given in Fig. (4) 

for modal solution time store of 10 and 500 y secs. The shot 

numerical integration time step was | и sec in each case. For small 

time step ratios 10:1, Fig-(49), tne results are satisfactory but at 

large ratios 500:1, Fig. (4b), the solution is different. This can be 

need by comparing the data calculation points (the corners on the graph 

every 500 u sec) on Fig. (4b) with the same points on Fig. (4a). This ๐ m 

is in line with what would be expected from the preceding results. Thus MO 

the Ist integral only algorithm is suitable for modest step size ratios : 

but begins to break down when very high ratios are used. 


The results for 154 and 2nd integrals are illustrated in Fig. (5) 
for step sizes ratios of 10:1 and 1000:1. Excellent aqreement is obtained 
between the two cases, the solution points of the 1000џ sec step size 
calculation falling on top of those of the 10 u sec step size calculation. 
Again this is what would be expected from the previous analysis. 


e 
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4.0 Conclusions 


An algorithm suitable for use ss а SHOCK-AiD interactive contre] 
routine has been developed that allows ง numerical integration routine 
to be coupled to a modal solution tecanique such that each can have á 
different vime step (in intoger ratio). This means that the numerical 
integration can solve the frst dynamics of à few non linear degrees от 
freedom (the shot) with û smal? cime step, wuüilst the modal SHOCK-AID 
solution can solve the slower dynamics of 4 much longer number of 
linear degrees of freedom (the barrel), using à much larger time step. 


The algorithm, has been programmes into SHOCK-AID along with the fuli 
projectile model given in Ref. (21 and is at present undergoing validation 
tests. 
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This paper addresses the use of the SHOCK-AID computer code together with a 
recently developed subroutine for modelling the effect of barrel curvature 
profiles on the dynamic response of the gun system. It is pointed sut that 
attempting to incorporate such details at the Finite Element stage suffers 
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INTRODUCTION 


1. Gun systems are generally analysed by treating the barrel аз а large 
elastic structure with many degreees of freedom, whilst the remainder of 
the system is viewed as a set of driving mechanisms. 


2. Generally these mechanisms are interactive with mutual loads depending 
on the current displacment аға state of motion of the whole system. 
SHOCK-AID (Reference 1) is a suite of computer routines for simulating the 
response of gun systems under general firing conditions. It has been 
designed in a modular fashion to accomodate evolutionary enhancements and, 
in particular, contains the facility for a user io introduce new mechanisms 
via an Interactive Contro] Routine (ICR). 


3. A number of mechanisms have been ccnsidered for ICR applications. For 
exampie: 
a. A Shot/Barrel ICR to simulate realistic loads on the barrel as 


the shot progresses toward the muzzle. This ICR has been outlined 
at previous symposia (Reference 2). 


b. A Barrel/Cradle ICR to simulate the contact/impact mechanism as 
the barrel recoils through the cradle bearings (Reference 3). 


A. This report describes the development of an ICR for predicting the 
effect of barrel curvature on gun response. Barrel curvature arises in 
many ways: Gravity droop; differential cooling; manufacturing set; etc. 
A problem arises in simulating barrel curvature numerically (eg: using 
finite elements or differences), since the effect of curvature will be 
swamped by numerical approximations unless the model is prohibitively 
detailed. Furthermore, the performance of parameter sweeps involves the 
generation of a new model for each distinct curvature profile. The Barrel 
Curvature ICR described here provides for the imposition of smal) 
deviations in the barrel profile as an input file and interactively 
corrects for the resultant changes in mass distribution as the event 
evolves. In this respect it differs from cther ICR's in proividing a 
correction for geometrical effects rather than adding further loading 
mechanisms. 


THEORY 


5. Until now gun barrels have been represented іл SHOCK-AID by beam 
elements with forces being applied to an undisplaced configuration. There 
are two problem types that can be dealt with by the Barrel Curvature ICR: 


a. Any initial shape of the structure must be included in the model 
explicitly at the finite element modelling stage. for example, to 
model the differing curvatures оп а set of barrels due to 
manufacturing set a corresponding set of models, one model for each 
curvature would have to be created. 
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b. Calculations performed in the main SHOCK-AID response routine 
are based upon the physical geometry of the undisplaced nodes. This 
may cause an approximation when the response of the beam is 
calculated because it does not account for the modified mass 


distribution. 
ASSUMPTIONS 
6. In order ว overcome these difficulties the following assumptions are 


made in the formulation of the ICR: 


a. The effect of small changes in mass distribution сап be modelled 
by the imposition of appropriate correcting couples to the 
undisplaced configuration. The corrections concern inertial and 
external loads only. 


b. For the purpose of calculating the correction the beam element 
may be treated as being rigid. Although tnis is not strictly 
necessary it is an extremely convenient simplification. 


с. The correcting couple can be applied аз а self-equilibrated pair 
of forces applied to the ends of a beam element. This, of course, 
relies on the assumption of rigidity. 


d. Тһе correcting couple is based upon the state of motion of the 
beam elements at the previous time step and any externally applied 
loads. 

e. The mass matrix is assumed to be banded. Again this is not a 


necessary but s convenient simplification. 
MATHEMATICAL FORMULATION 


7. Fiqure 1 shows the displaced апа undisplaced beams and their 
associated loading geometries. 


8. Consider the vector from A to the point of action 
re-of, = 1 - в, 


Then the length of the beam can be written as 
25; 0 


or in matrix form 


-а -1 1 
| 4) = Ir, KR Fol 1 
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From the point of action to the beam is denoted by 


d =—r+ 
ж 


~ 


AN 


=a A + x 


~ 


hence couple c = dA R za ( a + x) КЕ, + Е.) 


E (Ei AES АКЕ T Eo) 2 
9. Similarly when an element is rotated Бу аз shown in Figure 1, so the 
couple acting is: 


a 
where (89 


We knaw that 


:1 
ГЕ ЕЛЫ 


X6 = ( cos 8 sin dÉ 
-sin 9 cos 8 

t 91 
Using equation 1 and expanding putting К -— and Го = 2 
f2 92 

q. l4 => q. i. 

i.e. a ее 

927129 


Similarly for the rotating beam it can be shown that 


10. It follows that by subtracting the couple calculated for the 
undisplaced beam from that for displaced beam we can evaluate the couple 
discrepancy 


Š, = Cp - с 


Combining equations 2 and 3 the expression for the ccuple discrepancy 
becomes 


£ 


x g "ai (L A Fo) + Gg 7" x) A (F, + Fo) 


IMPLEMENTATION 


ll. The implementation of the ICR is illustrated in Figure 2. The 
routine is called from the main SHOCK-AID loading module. The global model 
data is passed on initial entry to the ICR. The routine then loops through 
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all the elements along the barrel successively extracting local mass апа 
geometry data, calculating the required correction couple and applying the 
appropriate statically equivalent self-equilibrated force pairs. 


12. The following is a list of requirements regarding tasks to be 
performed by the Barrel] Curvature ICR. 


a. The ICR should allow initial displacements for one series of 
connected nodes. 


b. Displacements anc accelerations shculd be used to calculate 
couples, Equation 4. 


c. The response routine (in SHOCK-AID) force vector should be 
modified according to the couple calculated for the beam between pairs 
of nodes. 


d. Absolute position data must be produced for output using plotting 
routines. 


e. Geometric, connectivity, acceleration, displacement and mass data 
must all be available. 


15. The last requirement, е, implies certain modifications to the 
SHOCK-AID data structure. A new data structure has therefore been 
developed to provide facilities for all forseeable applications. 


CONCLUSION 


14, The Barrel Cuvature ICR provides a simple and effective solution to 
the problems associated with parameter sweeping on model geometry, and 
also corrects for the approximations concerning mass redistribution 
inherent in linear finite element techniques. It constitutes an example 
of the way in which, with the enbanced data structure, SHOCK-AID can cope 
with complex geometries and and boundary conditions in a realistic manner. 
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The prediction of barrel behaviour during the firing phase has in 
the past been approached iu (ме distinct ways. The first, that of 
finite element analysis to determine mode shapes and tregquencies, 
followed by the use of these to give the transient respouse(c!,2). Тһе 
second, the solution of the basic Euler-Berüoulli, or Timoshenko, beam 
equations by a finite difference or other technique (3,4). Іп both of 
these methods ditticulties arise when non-lineoarities aro introduced 
pacticulariy nonlinear bearing stirfiessoso, clearances, or damping. 

lı order to overcome these dirticulties a dedicated росхаце has 
beea produced which treats the barrel as а тоз сі-Чергео of freedom 
spring Mass system. The resultant set or second order cifferential 
equations are then integrated directly. This technique enables non- 
linearities to be introduced easily апи previous investigations have 


included the ettect ot bearing elasticity, bearing clearance, bearing 


КА 


and barrel damping (5) ead more recently shot barrel interaction 
effects by Powell (o). 

This paper describes the principles ot the programme written for 
che SIMulation of BArre! Dynamics, SIMBAD, and in particular the 
inclusion and investigation ot the effeet which barrel expansion has 
upon the barrel response. The barrel is assumed to expand within the 
bearings until the clearance is teken up. This cau be introduced te 
the barrel modei as either а forced displacement’ at а bearing, or mere 
realistically, a generalised torse derived trom the bearing stillness. 
Batres sima au su 

In the Simulation package the barrei is divided into a munher of 


Euler-bernoulli team elements with three de;rees of freedom ai eoch 











PI 


ола. These relate to the transverse displacement, the longitudinal 
displacement, and the angular displacement. This provides a typical 
single plane analysis which is considered adequate for symmetrical 
barrels and shot barrel interaction etfect3 ot non-spiuniup 
projectiles. 

Foi 1 ว ง พ โท ม entry of the barrel external and internal profiles, the 
barrel is automatically divided inco a number of elements as specified 
by the user. This also takes account of bearing position. The 
resultant equations are assembled into а stiffness matrix in the same 


wav às a finite element stiffness matrix would be composed. The fina 


basic equation is of rhe form: 


where [X] is the stitíness matrix | 


[5] represents all of the degrees of freedom and 





[F] repre: ents the external forces at each node which can 
ง : - : & Ж | 
include inertia forces, damping forces, bearing forces, etc. ~ 
The present system allows a maximum of 25 elements (78 degrees of 
freedom) but experience has shown that a solution penerally converges 
at between 20 and 25 elements. 
It the inertia effects are included the equations become 


(M) [8] + [К] [6] 


-- 
=; 
et 

il 


Where [M] is the mass matrix of the system which in the package is 
calculated automatically together with the stiftness matrix. 
Additional masses such as breech mass and moment of inertia is added to 
the mass matrix and the resultant set of simultaneous ditferential 
equatious solved by a second order Кипде-Ки са technique. During this 


process it is relatively simple to suppress a displacement, or impose a 


Dier, displacement. 
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birrel Expansion Ettecty 

Previous bearing r ntations have generally either assumed a 
pin joint, o linea.lwv elo ic bearing or an elastic bearing with a set 
clearances The latter ot these ік probably the more realistic, but tt 
still does not take account of the fact that during firing this 
clearance is changing as the barrel expands. Also, the expansion 
itself will tend to make the centre line of the barrel move relative to 
the bearing. Tais Latter case on its own, can be represented simply as 
a prescribed displacement at a bearing node, such that the displacement 
ol the node equals tho barrel expansion. This represents au 
unreaiistic, idealistic case of expansion on a solid bearing. However, | 


when it is realised that the change in radius for free expansion of a ! 





typical 120 mu gua is of the order of 0.33 mm, it is of interest to 





examine tue case as a base line, 

The complete case is more complex because conditions change as the t 
clearance is taken up or the barrel loses contact with either siue of 
the pearing. kftectively the model of the barrel represents the motion 
ot the centre line of the bore. The expansion, or ілсгеаве in barrel 
radius at the bearings, can then be added or suotracted to give the 
etfect which the expansion hes upon the bearings. it the barrel is 
initially sitting on the bottom oi the bearing rben four separate 
conditions need to be considered, The 1121151 starting position is 
shown in Dicgram 1, and from this point the displecement of the barrel 
centre line at the bearing is Xe. The displacement of the lowest 


Surtace is therefore Xe - Ex, and ot the top surface Xc + Ex. 
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The tour separate conditions «sre now: 


Ms When Хе - Ex > C and Хе + Ex © СІ, the barre! vil! not be in 
contact with the bearing, and tle bearing force, RY = 0, 
de When Nc — Ex c U and Ne + Ex < СІ, the Lower surface of the 


barrel will be pushing into the bearing, and the upper surrace 
wiil be tree. The resultant torce acting on the barrel will 
ho: 
Bro = -(No - Ex).Kb 
3. When Xe -Ех £ 0 and Nc + E! > Cl, both the upper and lower 
surface will be in contact with rhe bearings. The force on 
the barrel will hive two components such that: 
ВЕ = “(XC = Ex).Kb - (Xe + rx - С1).КЬ 
ev 
Br = -(2Xe - CI).Kb 
It will be noted tnat once the clearance has been taken up, 
the force component due to any more expansion disappears. 


When Xe - Ex > 0 and Xe + Ex ә CL, the lower surface has risen 


+ 
. 


off the bearing, but thu clearance has been taken up on the 
other surface. Therefore, 
Bt = -(Xe + Ex -C1).Kb 
The above equations within the Limits indicated, give the force to 
be applied to © barrel at the bearing node if the barrel initially 
rests on the bottom of the bearing. A similar set of equations can be 
written when the barrel initially rests at the top of the bearing. 
The barrel expansion is actually a dynamic situatioun, but 
unpublished studios by P King at RMCS suggest that for all practical 


purposes the expansion сап be considered directly proportional to the 


internal pressure. For this study, a simple finite element model was 











3 


used to obtain an expansion ceetticient, Ce, and the expansion given 
bv: 
Ex = Ce.Instantaneoug pressure БЕ 

The value used tor the bearing stiffness was derived trom earlier 
work carried out by King used in studies of a 120 mm barrel, 
Computer Simulation 

The simultion of the barrel expansion was applied to a typical 120 
mm pun, similar to the Chiertain barrel. The pressure time curve was 
also typical of an ARDS round ín such а дай. Keeping the same bearing 
positions, the following cases were considered: 

1. Ne expanston with pin jointed bearings. 


No expansion with elastic bearings. 


3. No expansion with elastic bearings and clearance of 0.25 mm. 


As 3 abeve but with very large clearance ot .0 mm.so that 


г 


impact with the other side of the bearing does not occur before 
shot exit. 
5. Barrel expansion on solid bearings. 
о. An 3 above but including barrel expansion. 
7. Ав 6 above but with cloarauce of 0.15 mm. | CH 
8. As б above but with a clearance of 0.5 mm, 
9. As + above but including barrei expansion. 
10. A case with no transverse bearing cestraint,(Free-Free). 
The simulations all included the recoil acceleration and the 
resultant moments produced by the offset breech aud muzzle reference 
sight. 


Discussion of Results 








Figures l and 2 demonstrate the efrect of various bearing 


{dealisations witbout barrel expansion included. Figure | shows the 








š 


e 


vertical displacement of the muzzle tor pin jointed hearings. elastic 
beariuss, aud elastic bearings with cleavance. Figures 2 shows the 
displocement with typical clearance, very сауыс clear:uce and no 
Sadrin; Tsiraini at all. 

The Тагие clearance rase aud Free-rFree case produce very similar 
results. This 15 because the centre of gravity of the barrel is ahead 
oi the front bearine, and therefore it normally sits on the bottom of 
the trout bearing, aad against the top of the rear bearing. When the 
breech moment is applied, the terlency is to Litt away trom each 
hearin, and thererare have пэ restraint. However, there is an initial 
compression of the tront beariag tor a very short period which produces 
the slight: ditterence between the two causes, 

Figure 3 is shown tor completeness and dispiays the probably 
unrealistic case of expansion on solid bearings compared with elastic 
bearings and 10 mm ciearance. 

Figute 4 shows Lhe different responses for expansion with elast.c 
bearings and various amounts of clearance. What is particularly 
interesting is the flattening of the responses around shot exit as the 
erearance is reduced, This has been observed in actual firings when 
Measurement oi muz-leo тоуешеп was ботад taken. Lt occurred during à 
series ul firings іп relatively rapid suceesston. It could therefore 
have been caused by the barrel expandiug as ifs temperature increased 
and Pence Che clearance reduced, Direct comparisons between expansion 
and no expansion are shown in Figures 5 and b foc 0.25 mm clearance aud 
10 เพ ต clearance. [t can be seen coat the e pansion reduced the peak 
response ot the bariel іп both cases. 


Contidence in the results is obtained when the barrel detlected 
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shapes are considered. These are shown іп Fipures 7,4, 9 and 10, and 
Indicate the barrel shapes takon up at 0,5 ms Intervals. Тһе case ot 
bearing clearance with no expansied, Figure Z, shows (һа: at the ^ cont 
bearing the barre! moves across toe clearance and, attor ccnpressing 
into the other side of the bearing, їз det lected hack across the 
Clearance, When expansion is іпеншіеі, Fipure 8, the barrel still 
moves across the clearance, but because the barrei has expended and the 
clearance therefore reduced, the centre ot line ot the barrel does not 
yet detlected back so easily. 

The effect of the rear bearing is more obvious. The motion ОГ the 
barrel dowuwards is assisted b the expansion of the barrel against the 
top or the rear bearing. The downwards motion of the breech is 
therefore increased by nearly 1.4 mm, 

The case ot a simple eiastie beiriny with no clearance ав shown in 
Figure 9. The motion at each bearing has obviously been suppressed, as 
has the overall motion of the compiete berrel. Figure 19 shows the --- 
Егес-Егес case tor comparisons 

Finaliy Figure Ll shows the jemp characteristics for a simulation 
with and witaout expansion. The figure ot jump is calculated as the 
“usalo angle plus the muszle transverse velocity divided by the shot 
exit velocisy. This is shown tor the complete time period although the 
tain area o- interest is at shot exit, It enables an appreciation to 
be made of tie possibie change in jump if the shot exit point was not 
correct. [t сап be seen thet for the simulation with barrel expansion, 


the exit time is not so eriticat. However, 0.5 ms difference for the 


no expansion case would produce 2 very difrerent jump figure. 





Gone Pus tons 

The effect of barrel expansion within the bearings and its effect 
on transverse response can be simulated relatively easily. Lt can be 
defined as additional forcing terms to the barrel, dependent upon the 
displacement ot the barrel withiu the bearing. 

Results show that inclusion of the barrel expansion prodvces a 
significant change in the barrel response and should therefore be 


considered in studies ot gun jump characteristics. 
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A SIMPLE THEORETICAL MODEL OF ZHO;/BARREL INTERACTION WITHIN A 5МООУН 
BORE GUN 


ABS PRACT: 


A model of a smooth bore barrel in one plane is developed. 
This involves the use of stiffness equations of single peam elements, 
similar to those used in finite element modelling. The equations 
of motion describing the pitch and bounce modes ol the shot are 
stated and the two sets of equations are linked to predict the 
effects on the shot of shot-barrel interaction. 

Data from a compressed air powered gun is used to validate 
the model. 
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Le INTRODUCTION. 

When considecing modern tank guns one oí the factors contributing 
to round to round dispersion is the yawing motion or the shot. Indeed с 
Andrewa(!) suggests that this motion is of major importance and 
Walker(2) makes the point that shot balloting (rattling in the bore) 
can cause Severe sheli engraving and wearing ot the gun barrel. To 
understand the cause ot гам for зи aerodyninically stable shot, it is 
necessary to studs barrel movements, projectile tn-bore motion and the 
interaction of the bare) and tue projectile. 

А considerable amount of experimental and theoretical worx has 
token place in (ін Cicld at dirrel motion. Mathematical modelling 


techiiques usiny computers are now being used. These can be divided 


iuto Caose using білісе ditverence equations, for instance Pagan(3), 





and Ming and T^oaas(a) and those using the finite element method,such 


as Taclor,Thonaasen and Vanee(5) in Great Britain and Soifer and 

Becker(6) ia the United states. çv 
А Tew theoretical borcei-projectile modeis have been produced. 

Notably by Soo Но» and Anderson(7) whe used equations of motion for the 

ргә зер! donamiess otogiued using Lagrange’s method, coupled to а 

finit. difference navel model, and Chu(8) who used Euler’s approach to 

Odtala equations of wenion for tie shot. The potentially more powerful 

finite Төлеп barrel models produced in Great Britain and the United 

States are being modiried to include the cffects of shot/barrel 

inberacib:on, 
Another approach, Suggested by Bulman(9) is to produce the 

stitinmess equatiosus similac to the finite element technique, include 

the inertia forces, and then озун the cesulting set or simultaneous 


second order difterential squat ious using one of the standard 


аїдог шә aviriable. 








ln this paper Bulman’s approach is used to produce a model of a 
smooth bore barrel in one plane. The equations of motion describing E 
the pitch and bounce modes of a shot in one plane are stated and the 
two sets of equations are linked to predict the effects on the shot of 
shot-barrel interaction. 

À compressed air powered apparatus which simulates many of the 
characteristics cf a recoiling tank gun has been developed at RMCS(10). 

Data from this gun is used to test the theoretical results. 


DERIVATION OF THE MATHEMATICAL MODELS 





Barrel modelling 


A two dimensional mathematical model of a barrel with constant 
bore and outer diameter was produced by considering a series of 
Euler-Bernoulli beam elements joined at the nodes. At this stage 
longitudinal vibrations were not included. This gives two degrees of 
freedom at each node, These are linear motion perpendicular to the 
barrel axis and angular motion about the horizontal axis. Using the 


sign convention shown below for the displacements, 
ز۷‎ 


Ө) 


\ 


node Uj = 


Fig 1. Displacement of node J 
and letting the applied external forces on che barrel element at nodes 


j aud jtl be shear forces Fy and F 1 and moments М. and Hui: with the 


jt J 
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There is no requirement that the barre! elements need to be the 


same length. Indeed the program chooses element lengths s9 that the 3 
bearings аге at barrel nodes. 
Shot modelling 

We first make the following assumptions. 5 


1. There is no friction between barrel and shot. 


2. The driving bands have no play iu the barrel. 


small such that sin 8 = 9 


and 3. Yaw deflections are 





Fig 3. Schematic diagram of a shot in a smooth bore barrel, 


Qe Using the schema: іс diagram (figure 3) and considering shot 


‘bounce’ and ‘pitch’ we can write the equations of motion for the shot. 


My = -k.(y = GL) = kply + ai = С. (У = OL) ~ Cy + Bal ............ (4) 
OL = T+ L|k,(y - OL) ~ Loky(y + 61.2) + C.L (у ~ 0L/) - BEE i ves (5) 


These equations are solved using a 2nd order Runge-Kutta routine 


and the shot and barrel solutions ar: connected using the following 


method. 
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The position of the shot in the barrel is calculated from tie 
pressure time curve and the position of each of the driving bands 
relative to elemeuc nodes is deduced. 

a) Forces from the barrel on the shot 

To include barrel effects in the equations of motion for the shot 

it is necessary to interpolate between the barrel nodes to give values 


for barrel motion at the driving bands. 


1+3 





Fig ^. Displacement of the barrel at che onot driving bands. 
Using figure 4, let the rear driving hind be between nedes x; and хі 


and the front driving Баса be Болысси anodes xj aud WI then, 
L 


= : * Ke m Cs } ft 
Ха хі ы Ui] Kd al 


t = he + А - ха, /t 
b j (хы; US 
The vertical distance moved by the barrel is then incorporated into the 


shot bounes equation (4) to give the new equation, 


Hj = -k (y-x Bh = К (у-х * OL.) = C (eA š 
bo? POSS EOL ри mX —©], = С T + Or с 
а а 1 b L 2 На | €, =, 90.) |... (6; 
The eqnation for shet pitch becomes, 
`. T ` ` ` 
Jj = P+ J k (y-x “Al, = L. K Ven ке LEE GSS E Y, uum แว MATS A WER 
ha Ee BN OS GG EE ล. 





is 





D korce From the Slot son Che barrer 
Vertical forces from a driving band are added to the barrel nodes 
either side of that driving band in proportion to the distance of the 
driving band from the nodes. Using figure 4 and letting Fib be the 
shear force from a driving band a distance d, from node i the shear 
force on node i is calculated using the formula, 
М, = ธิ น ส (1 = d./L) 


and the shear force from the same driving band on node ізі is. 

Nig, = Fap(d4/L) 
So that if the shot is within one barrel element, all of the forces 
from the shot are regarded to act at the nodes either side of chat 
element. But if the driving bands are in separate barrel elements the 
forces are applied at three nodes. 
EXPERIMENTAL PROGRAMME 

A 30mm, muzzle loaded, compressed air powered gun has been 
developed at RMCS to simulate many of the characteristics of а 
recoiling tank gun. Many parameters can be varied easily, these 
include breech mass, breech inertia, buffer ;tiríness, bearing spacing, 
shot mass, bearing clearance and breech pressure/time profile. Barrels 
can be interchanged within a few minutes and these include a smooth 
bore mild steel barrel with a wall thickness of 10 cm, a thin walled 
aluminium smooth bore barrel and a rifled barrel. 

The gun body is symmetrically balanced about both the horizontal 
and vertical axes but masses can be added above or below the axis to 
create an out of balance ‘breech’. 

Instrumentation has been designed which records the change іп 
angle of the shot, about both the transverse horizontal axis end the 
vertical axis. To facilitate this the round was required to carry a 12 


mn diameter mirror on its front face. 
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A round has been «сзіспей(11) to give а set of characteristics 
that include, 

l, a pitch frequency of between 150 Hz and 300 Hz, to vive 
sufficient number of cycles as the shot passes down the 
barrel, 

2, а maximum pitch deflection of +4 mrad, 

3, a damping ratio of .l, or less, 

and 4, symmetry about the central axis with provision to 
attach oft axis masses to create an out of balance 
round, 

А cross section of the design is shown in figure 5. It consists 
of an outer cup (АА) with machined driving bands (BB) апа (CC), а l2umm 
diameter inner beam (D), a cantilevered beam (E) and an end inass (FF). 
The position of the mirror to be used when measuring the angle of the 
front face is shown at (G). In the experiments described here the 
cantilevered beam (E) was specified in silver steel, with the rest of “- 
the shoc in aluminium alloy. However in future shots beams "D and ‘Е’ 
are to be integrally machined from one material. The ‘free’ end of the 
beam is shaped to give a larger contact area with the mass. This is 
attached by a press fit and adhesive. 

The driving bands were specified with a tolerance of .01 mm, this 
fit was designed to form a pressure seal and was considered tight 
enough for the balloting motion (rattling in the bore) of the shot 
outer cup to be neglected in theoretical predictions. 

Except for two holes in the end mass (FF) the design is 
symmetrical about the central axis. Various masses can be screwed into 


these holes, to make the shot symmetrically balanced about the bore 


axis, or to give it an off-axis centre of gravity. 








Le 


To give an initial test of the mathematical mode! two sets of 
firings of shots with a 2.1 gram mass, offset from che shot axis by 
10.5 mm and positioned below the bore axis prior to firiug, were 
chosen.These used, 

1) The mild steel barrel with a muzzle support fitted so that the 
barrel can be regarded as inftexible and the shot studied in 
isolation, aud 


2) The flexible (aluminium) barrel with a balanced “breech”. 


RESULTS AND DISCUSSION 


Figure 6 shows two theoretical traces for a shot with a 17mm x 
2.6um diameter cantilevered beam with an offset mass of 2.1 grams 
placed 10.5mn underneath the axis of the shot. One trace is for au 
inflexible barrel and the other for the flexible barrel. The 
theoretical predictions suggest that barrel motion will have a definite 
effect on the shot, In this case the effect of the flexihle barrel is 
to rotate the nose of the shot downwards about the horizontal axis. 

Im figure 7 the theoretical prediction and an experimental result 
from the ‘inflexible’ barrel are plotted. As this figure shows the 
theoretical model is now giving reasonable predictions of shot motion 
within a solid barrel. The experimental trace shows a frequency of 
approximately 222 Hz, whilst the theoretical trace has a frequency ot 
about 210 Hz. 

Theoretical and experimental traces for a shot with a 20mm x 2.6mm 
diameter cantilevered beam, fired from the flexible barrel, are plotted 
in figure 8. The experimental result confirms that a downwards 
rotation of the shot позе, about the horizontal axis, does occur when 


the inflexible barrel is replaced by the flexible barrel. 
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Опе cause of inacurracy in the flexible barrel model is that the 
gun bearings were modelled as pin joints, and it has been shown(9), 
that the bearing stiffnesses and in particular damping at che bearings 
is an important factor when predicting barrel response. 

A large scale firing programme to provide data to fully test the 


mathematical model will now be implemented. 


CONCLUSIONS 
Early comparisons with a few experimental results show encouraging 
results. 
One of the advantages of the experimental apparatus described is 
that a large number of firings can be undertaken cheaply and quickly. A 
firing program will now be undertaken to provide large amounts of data 
both to .est the simple mathematical model described and for the 
testiug of other theoretical models. 
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NOTATION 


For the barrel. 


— F; Shear force at node j 
Mj Moment at node j 
E Modulus of elasticity 
I5 Second m ment of area of the barrel cross section 
Lh Length of barrel 
L Length of a barrel element 
Uy Horizontal displacement of node j 
V; Vertical displacement of node j 
6; Angular displacement of uode j 
Mpe Mass of a barrel element 


For the shot. 
0 Angular displacement 


Lj Ly Position of the body centre of gravity relative to 
, 


d the front and rear springs. 
y Vertical displacement of the shot centre of gravity, 
T Couple produced by the effect of an offset centre of 
gravity. 


kak, Stiffness coefficients for the driving bands 

М Mass of the shot body 

I Moment of inertia of the shot about its centre of gravity 
СазСь Damping factors for the driving bands 


Shot-barrel interaction 


Nj Shear force on barrel node j due to a driving band 
Ғаһ Гогсе from a driving band 
а; Distance ot the driving band from а barrel node 


l 
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ABSTRACT: 


The sources of error in a representative modern tank are considered. 
Those processes occuring after shot ejection are treated in detail. Muzzle 
blast, sabot discard, and free flight aerodynamic perturbations are 
quantified. The influence of fire control system error is investigated. 
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ANALYSIS OF SOURCES OF ERROR IN TANK GUN FIRING 


EDWARD M. SCHMIDT, PhD# AND JOSEPH W. KOCHENDERFER, MR 
LAUNCi] AND FLIGHT DIVISION 
U. S. ARMY BALLISTIC RES ARCH LABORATORY 
ABERDEEN PROVING GRCUND, MD 21065-5066 


1. INTRODUCTION 


The capabilities of modern tank ammunition can be fully realized 
only when the rounds are delivered accurately. A primary source of error is 
the ability of the fire control system, including the gunner, to sense and 
locate the target then to lay the weapon to the proper orientation to insure 


ballistic intercept. However, even if this operation is performed to 
perfection, there are a number of perturbations which influence the launch ard 
flight of the projectile. Some involve the dynamics of the 


tank/gun tube/projectile during inbore acceleration and subsequent separation. 
Others are dominated by the prejectile aerodynamic characteristics. Also, the 
inability to sense and/or correct for the ambient conditions over the 
trajectory will cause disparity between the desired and actual points of 
impact. The present paper will address some sources of error in a 
representative modern tank, mainly considering perturbations following 
separation of the projectile from the launcher. 


An illustration of some of the problems facing the analyst of 
in-bore dynamics of modern projectiles is provided in Figure 1. The weapon is 
a masSive guide rail for the projectile; however, under the extreme loads of 
launch, the gun tube begins to respond while the shot is still in bore. Іп 
order to force a consistent sense to this response, cannon are generally 
fabricated with mass asymmetry. The tube is supported and constrained by the 
recoil system which may or may not be symmetric in design. Under gravity, the 
tube droops producing a curvature which is complicated by variations in bore 
straightness induced by manufacturing, installation, use, or the environment. 


Upon firing the charge, the pressure within the tube builds up to 
thousands of bars causing both the projectile and gun tube to accelerate. 
Fach of these moving bodies is subject to a gross rigid body motion upon which 
is superimposed flexural or vibrational modes. The linear and angular 
acceleration of the gun tube induces inertial loads which contribute to 
flexing. In addition, the motion of the projectile along the curved path of 
the tube results in transverse reaction loads which drive the gun motion. 
Historically, projectiles are assumed to undergo rigid body balloting motion 
as lateral clearances are taken up and rebound occurs. With long rod kinetic 
energy ammunition, this response is supplemented by longitudinal flexure of 
the sabot and of the projectile relative to the sabot. 


In terms of accuracy, the critical moment of the in bore cycle is 
the separation of the projectile from the gun tube (Figure 1). Anything that 
influences the magnitude or direction of the projectile velocity vector could 
resuit in significant deviation from the intended trajectory. If the launch 
process terminated at the muzzle, the parameters of interest would be the 
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transverse velocity, angle, and angular rate of the gun tube plus the 
transverse linear and angular velocity cf the projectile relative to the gun 
tube. Unfortunately, the projectile must disengage from the tube, pass 
through the reverse flow region of the muzzie blast, and then discard the 
Sabot componenis before entry into unconstrained free flight. During transit 
of this near muzzle region, the trajectory may be further perturbed. 


Once clear of the blust ani sabot discard regions, the projectile 


free flight motion is reasonably well understood. Assuming that the 
projectile nas not been damageu during launch, the trajectory сап be well 
determined if the initial state із defined. The main source of error 


associated solely with the rree flight of the round can be ascribed to the 
inabilicy of the fire control system to accurately sense and provide 
appropriate corrections for ambient conditions over the flight path. 

In order to complement the papers cescribing the contribution of gun 
and projectile in bore dynamics to launch conditicns, the present paper will 
discuss the environment seen by the projectile following separation from the 
gun tube. Muzzle blast loads will be reviewed for weapons with aad without 
muzzle brakes. Sabot discard interactions wiil be considered and related back 
to in bore properties. Finally, the free flight motion will be addressed. 


2. MUZZLE BLAST LOADS 


When problems are encountered in launching rounds from cannon, it is езі 
a common practice to ascribe the difficulty to loads encountered in the muzzle 
blast. Іп defense of this interesting gasdynamic phenomena, it must be made 
clear that care is required in diagnosing the origin of weapon launen 
problems. Muzzle blast may be the source of disturbance in some instances, 
particularly, for objects which are bluff, asymmetric, or of low density. 
However, in many cases, the extreme loads characterizing the in-bore 
environment are the real source of difficulty. 


Once free of the gun tube, the projectile must transit the muzzle 
exhaust flow which consists of an outer air blast driven by the expanding 
propellant gas plume (Figure 2). The piume has the structure of a supersonic, 
underexpanded jet and is terminated at the shock layer consisting of the plume 
Mach disc, the propellant gas/air interface, and the outer shock of the blast 
wave. While inside the Mach disc, the projectile is subject to the plume flow 
wherein the gas velocity reaches values 2-3 times that of the projectile, 
1.е,, the projectile is in reverse flow. Conventional wisdom dictates that 
for such a situation, finned projectiles would be unstable, but since the 
residence time within the plume is much lower than the inertial response times 
of the body, stability is not really of concern. Rather, it is the transverse 
impulse transmitted to the projectile that influences the subsequent 
trajectory. 


The calculated variatior. [1] of the lift force, Us with distance 
from the muzzle is illustrated in Figure 3 fer different values of the ratio 
of the projectile exit velocity to the propellant speed of sound at shot 
ejection, Vp/c, i.e., the Mach number of the propeliant gas prior to shet 
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exit. For values of the parameter equal to or greater than one, there is no 
in-bore interaction possible. For values less than one, the in-bore 
interaction occurs and may be important. Outside of the weapon, the behavior 
of the lift function is of interest. The relative flow velocity increases as 
the round moves through the plume causing the lift on the fins to increase. 
However, this increase due to the velocity, is offset by the rapid decrease in 
density of the propellant gas due to expansion. Under the influence of these 
two property variations, the lift function peaks aad then decays to negligible 
values after only two exit diameters of travel. 


Integration of the lift permits the computation of momentum 
tranferred to the projectile and the resultant deflection of the trajectory. 
The estimated trajectory deflection for a generic kinetic energy projectile is 
plotted as a function of launch velocity in Figure 4. The deflection in mils 
is ratioed by the launch angle of attack of the projectile in degrees. The 
lack of monotonicity in the plot is due to the assumed muzzle exit conditions 
which were taken from a variety of data sources. However, the result is 
striking. As the launch velocity increases, the trajectory deflection 
decreases significantly. For a typical tank round, the іл-роге yaw angle 
should be quite low. Even if the angle were as large as 0.5 deg, the jump 
induced by muzzle blast would amount to only 0.05 mils. 


3. MUZZLE DEVICE EFFECTS 


When а weapon is equipped with а nuzzle device, analysis of ine flow 
becomes more difficult. A common device, the muzzle brake, reduces gun recoil 
by venting the propellant gas rearward. Tle device has a  three-dimensional, 
confined geometry through which the propellant gases expand over the 
projectile. For sabot encapsulated rounds, the installation of a muzzle brake 
on the gun can create serious problems. First, the muzzle brake must not 
interfere with the sabot discard. Mechanical contact does extreme violence, 
both to the brake and to the round. Another, major difficulty can be the 
enhancement of muzzle gasdynamic loadings. Since muzzle brakes are under 
consideration for installation on a number of low recoil tank gun systems, it 
is of practical interest to consider the possible interactions. 


Recently, tests were conducted on a meaium caliber cannon with and 
without a triple baffle muzzle brake in place [2]. When the brake was 
installed, the dispersion of the system doubled. To determine if the flexural 
characteristics of the tube were altered, a mass simulating the brake was 
mounted at the muzzie. The mas: did not interfere with the free expansion of 
the muzzle gases. Firings делспвбгабей that the dispersion returned to the 
level obtained with a bare tube indicating that the added mass was not the 
cause of the problem. A further investigation was conducted to investigate 
the influence of brake length. Successive baffles were cut off the device. 
As the baffles were decreased from three to one, the dispersion also 
decreased. It was found that a single baffle design produced по measurable 
increase in dispersion over the bare muzzle case. Apparently, enhanced 
gasdynamic loads associated with the confinement and length of the brake were 
altering the launch dynamics of the systen. 
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Computations cf the gasdynamic loads on the projectile as it passed 
through the brake, predicted a 20% growth in dispersion. Since dispersion 
doubled, this was not the sol effect. А series of orthogonal flash 
radiographs were taken of the sabot discard process for cases with and without 
the triple baffle brake in place. These data show that the presence of the 
brake produced a change in the sabot discard trajectory which was not treated 
in the computations. Mechanical contact between the sabot components and 
projectile occurred over a longer period and asymmetry in sabot discard was 
enhanced. The resultant asymmetry in the discard increased the loading upon 
the projectile and degraded precision. The tests demonstrated that care must 
be taken in the installation of muzzle devices upon weapon systems which are 
designed for precision fire. 


4. SABOT DISCARD 


The fact that sabot discard interactions influence the trajectory of 
fin stabilized projectiles can be demonstrated by considering measurements of 
free flight yawing motion (Figure 5). The data were acquired in the BRL 
Transonic Range on a typical kinetic energy projectile. The plots present the 
angle of attack versus the angie of sideslip as the round moves through the 
200 m facility. With distance downrange, the yaw level decreases; therefore, 
the maximum yaw occurs near the weapon. Two blots are presented representing 
a case with low and high sabot discard perturbations, Figures 5a and 5b, 
respectively. The nature of this categorization can be seen if one considers 
the condition of the projectile at separation from the pun tube. 


Within the Lore, the projectile is constrained by clearances to yaw 
leveis on the order of 0.1 deg; although the angular velocity may be 
appreciable. If a statically stable projectile enters into free flight with 
such initial dynamics (i.e., near zero initial yaw, but finite yaw rate) the 
resulting angular motion should be nearly planar yaw. In fact, this type yaw 
15 well represented by the data shown in Figure ба. This plot would suggest 
that sabot discard interactions either were not significant or are consistent 
with norwal free flight aerodynamic loads. 


In contrast, a second round of the same type demonstrates the effect 
of sabot discard interaction (Figure 5b). Here, the yawing motion is not 
planar, but elliptical. If спе were to postulate initial launch dynamics for 
this round, they would be initial yaw of roughly 4 deg and a finite yaw 
velocity orthogonal to yaw angle. Obviously, the projectile did not have a 
4 deg yaw at release from the gun tube; nor is it reasonable to assume that 
disengagement from the tube forced the yaw to build to such a level. Rather, 
the sabot discard process produced perturbetions which significantly altered 
the launch dynamics of the round. 


The sabot discard perturbations consist of both mechanical and 
aerodynamic interactions. At the muzzle of the weapon, the sabot components 
and projectile are in direct mechanical contact due to the constraints imposed 
by the tube and the various bands or seals of the sabot assembly. After 
clearing the tube, elastic decompression, spin, and gasdynamic loads act to 
break these bands and to lift the sabot away from the projectile. Depending 
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upon the design of the sabot, the lift process may be rapid and clean ог may 
involve pivoting about a point of contact on the projectile. Alternatively, 
the sabot components may initially breal. contact only to reimpinge оп the 
projectile at a later stage of discard. Asymmetry in the contact will 
generate lateral momentum transfer between the sabot components and the 
projectile. This alters both the projectile trajectory and the symmetry of 
sabot discard. Geometric asymmetry in the sabot discard results in 
aerodynamic asymmetry in the mutuaily interacting flowfields associated with 
the sabot components and the projectile. 


The relative magnitudes of in bore and sabot discard perturbations 
have been estimated using data acquired from x-ray measurements of near muzzle 
projectile motion [3,4]. The results of a comparison of transverse angular 
impulse are presented in Table 1, below. The data indicate that the level of 
impulse due to sabot discard perturbations, is the same as that due to in bore 
disturbances. This conclusion is supported by data acquired by Biele [5] who 
measures the dynamics of a 120mm gun tube and the subsequent projectile 
trajectory. He finds that the gun muzzie motion accounts for roughly one-half 
of the measured trajectory jump angle. The remainder 13 associated with 
disengagement and sabot discard dynamics. 


Table 1. Comparison of Transverse Angular Impulse 


Round No. In Bore Angular Impulse Sabot Discard Angular Impulse 
(rad/s) (rad/s) 
1 6.37 5.50 
2 4.17 10.26 
3 3.03 ‚21 
4 3.58 3.28 
5 7.45 1.82 
6 1.71 6.30 
T 14.28 5.88 
8 1.29 9.04 
9 5.97 9.08 


Since sabot discard can significantly alter the projectile 
trajectory, it is of interest to consider the crigin of asymmetry in discard. 
Conceptually, there could be asymmetry associated with the failure of bands or 
retaining rings; however, this behavior is difficult to measure or model, 
Alternatively, the asymmetry could be directly related to the initial state of 
the sabot and projectile upon disengagement frum th ibe. During travel 
within the gun, the projectile moves along a curved patı  Jetermined Бу the 
tube bore profile and tube dynamic response. In addition, the projectile and 
sabot may be oscillating both relative to the bore and to each other. 
Finally, during release of constraints at shot exit, transverse loads may be 
imparted. 


Plostins [4] assumed that a major factor influencing the sabot 
discard interactions was the in-bore oscillation of tne projectile relative to 
the sabot. He examined the magnitude of measured sabot perturbation relative 
to the parameter © '/ + , For small values of this ratio, the projectile is 


iv-99 








SCHMIDT AND KOCHENDERFER 


assumed to be near the peak of its in-bore yaw. Conversely, large values of 
the ratio imply the projectile is being launched near the minimum of in bore 
yaw. The correlation indicates that sabot discard interactions are large for 
small values of fie / ฉั น ° This would correspond to ล condition where the 
projectile angular displacement is greatest and, presumably, where the 
differential strain in the sabot components is appreciable. The minimum 
discard interaction is measured for the case of large ^,' / ^. , which should 
be the case for low differential strain in the sabot components. 


The results point out the need for a design tradeoff. To reduce 
aerodynamic jump, the launch angular rate must be minimized; however, 
Plostins'! correlation suggests that as the rate is reduced the sabot discard 
interaction begins to build up. This implies that there could be an optimal 
launch condition for the minimization of the sum of the sabot discard 
interaction and aerodynamic jump. 


5. FREE FLIGHT 


In this section, consideration is given to the influence of the free 
flight aerodynamics upon accurecy. First, aerodynamic jump will be discussed. 
Second, the manner in which exterior ballistics is implemented in the fire 
control solution is addressed. 


Aerodynamic jump is the deviation of the trajectory associated with 
the yawing motion of the projectile. As yaw builds up, a lift is produced 
which results in lateral acceleration and displacement of the round. Murphy 
[6] integrates tne equations of motion to produce the following expression for 
the aerodynamic jump: C 


d L. [ , 
0 = — £' + HE 
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To examine the importance of the initial angular rate and angle of attack, 
consider two cases representing the launch conditions cf a fin-stabilized 
projectile: 


Case 1: É. = 5 deg, 5 = 0 rad/s 


Case 2: & = 0 deg, - = 14.5 rad/s 

To give yawing motion which is similar in magnitude (but not in 
phase) for each case, the first maximum of yaw will be roughly 5 deg. For 
Case 1, the projectile is assumed to be launched at this angle. For Case 2, 
the projectile is launched with an angular rate which will cause the yaw to 
build to a maximum of roughly 5 deg. Іп both cases, the initial velocity 
vector is assumed tc be along the desired lay angle. 


The resulting trajectory for each case was computed using a 
six-degree-of-freedom code (Figure 7). It is observed that the jump for 
Case 1 is negligible; however, for Case 2, there is significant trajectory 
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deflection. Since this case is more nearly that of typical gun launch, the 
computations indicate that the yawing motion of projectiles must be considered 
as an important source of deflection. This is especially true for yaw levels 
which are greater than one or two degrees. For comparison, the value of 
aerodynamic jump predicted by the above equation is shown on the plot. The 
two v^lues are in good agreement. 16 is interesting to note that trajectory 
def lc. ‚п angle is sensitive to projectile yaw for ranges less that 200 m. 
Beyond this distance, the jump approaches its asymptotic limit. 


The second free flight effect to be considered is associated with 
the implementation of gun lay. The fire control system, inclusive of the tank 
crew, must sense the ambient conditions and correct the pointing angle of the 
gun to provide the correct trajectory. If the ambient conditions are not 
properly sensed or input to the fire control computer, then the round will 
deviate from the desired point of impact. Since the flight characteristics 
are sensitive to the projectile design and launch conditions, two typical fin 
stabilized tank rounds аге considered: full bore, spike-nosed and pa 
sub-caliber, saboted long rod shapes. For a given level of uncertainty in “ 
ambient conditions, the shift in vertical target impact location is given as a M 
funetion of range in Table 2. Eo 


Table 2. Shift in Vertical Target Impact Location (in mils) due to 
Error in Fire Control Inputs where FB = fuil bore round and 
SC = sub-caliber round. 


De Range 10% error in 10% error in 10 ш/о แน 2216 5 ш/о Crosowind - 
Qe (m) density ambient temp. velocity error estimation error m 
FB SC FB SC FB SC FB SC 

0 0 0 0 0 0 0 0 0 > 

500 ‚018 ‚002 . 0011 ‚000 „036 ‚012 ‚290 ‚035 

1000 ‚090 ‚005 „020 .002 ‚083 ‚022 ‚632 O72 

1500 .245 .013 .059 . 004 .139 .035 1.003 .107 

2000 .560 ‚025 . 141 .008 ‚220 .047 1.484 ‚146 

2500 1.114 . 010 .286 .013 .321 .060 2.003 .182 
3000 2.184 . 060 .558 .021 478 ‚075 2.718 .226 

3500 4,002 ‚08 1.000 ‚030 „686 ‚087 3.490 ‚268 

4000 6.941 .116 1.617 O42 1.010 . 104 4,545 .310 


It is readily apparent that the sensitivity to errors іп fire 
control inputs is significantly lower for the sub-caliber round than it is for 
the full-bore round. In addition, the variation with range is highly 
non-linear. To gauge the importance of these shifts in impact location, it is 
necessary vo consider whether the estimation errors аге reasonable, For 
example, a 10% variation in air temperature is about 30 deg C from the 15 deg 
C standard. This could reasonably result from a seasonal excursion or from a 
day to night shange in a high desert environment. A 10% change in air density 
is at the extreme when variations are considered at a given altitude; 
however, if the weapon is moved from sea level to 1 km altitude, such a 

om density change is possible. А 10 m/s variation in muzzle velocity Пав been 
Г observed in tank guns even with a correction for propellant temperature and 
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has been attributed to uncertainty in tube wear and propellant aging 
characteristics. A 5 m/s error in tne estimation of crosswind is an upper 
bound on this property. 


Due the short engagement ranges of tank cannon, the influence of 
rotation of the earth on the trajectory is neglected. The magnitude of this 
effect on the azimuthal fall of shot is illustrated in Fig. 7. Even at 3 km 
range for the slower round, the correction is relatively small. 


The algorithm in the fire control system used to aim the weapon also 
influences the ability of the gunner to successfully engage targets. To 
compute gun lay angles, current systems use either а polynomial fitting to 
standard conditions plus corrections (unit effects) or solve two-dimensional 
equations of motion (with a closed form add-on equation for the cross plane). 
As long as computer memory space permits and the conditions are not widely 
different from standard, both solution methods yield comparable results. 
However, if there is a reasonably large divergence from standard such as a 15% 
density decrease coupled with a 10 m/s cross wind at an engagement range of 
2000 meters when firing a full caliber finner, the interaction of density and 
wind, if uncompensated, would introduce an error of about one mil in the cross 
plane. Fire control algorithms which utilize polyncmial fits for nonstandard 
conditions could suffer from the nonlinearity of the corrections but this 
shortcoming is minor and, generally, may be ignored. It should be pointed out 
that the fire control system attempts to compensate for the nonstandard 
effects addressed in Table 2; however, there is error in measuring each of 
these quantities which will influence the fall of shot. The perturbations 
listed, such as 10% in air density, can be properly accounted for if the 
correct value is input to the fire control. Unfortunately, the default values 
are often used in the solution of the aiming data because better data are 
unavailable or because the importance of these quantities is not recognized. 


6. SUMMARY 


The influence of perturbations to a projectile trajectory following 
separation from the gun tube is examined. Muzzle blast, Sabot discard, and 
free flight effects are considered. Muzzle blast has a minimal influence. 
Sabot discard can cause changes in the trajectory similar to those due to 
transverse loads within the gun tube. Generally, free flight loads are well 
understood; however, uncertainty in ambient conditions can generate errors in 
the lay of the weapon which can be Significant in some cases. 
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Figure 1. Schematic of Launch Environnent 
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Figure 3. Variation cf Transverse Force Through Muzzle Flow Field 
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Figure 5a. Measured Yawing Motion for "Low Disturbance" Launch 
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Figure 5b. Measured Yawing Motion for ‘Large Disturbance" Launch 
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Figure 6. Correlation of Sabot Discard Disturbance witn Launch State 
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Figure 7. Jump Angle as a Function of Range for Two Different Launch 
Conditions 
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Figure 8. Azimuthal Deflection Due to Rotation of the Earth for 3 km Range 
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ABSTRACT: 


Computer graphics is rapidly developing as ап engineering desipn and 
analysis tool, The advenc of interactive computing aad array processing оп 
minicomputers has made it possible to iaterface both experimental data and 
computer simulation with graphic displays to facilitate analysis. The use of 
color as a moaitor of stress level or thermal effects, for example, provides 
an enhanced capability for the recognition aod understanding of ballistic 
phenomena. 


This paper describes a computer code which provides three-dimeasional 
color output to display the dynamic motion of a gun tube, Using a description 
of the tube from a database and the centerline of the tube at rest, the static 
tube is displayed in both two and three dimensions. Data for the ceaterline, 
obtained from either a gun dynamics modeling program or experimental data, are 
used to calculate the gun tube display coordinates at a given time, 
Successive frames can be generated at appropriate time intervals to be used to 
make a film of the tube motion. 
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А THREE-DIMENSIONAL ОЛОК COMPUTER GRAPHICS PROGRAM FOR 
OF PROTECTILE GUN DYNAMICS 


DISPLAY 


KATHLEEN L. ZIMMERMAN 
U.S. Army Ballistic Research Laboratory 
Aberdeen Proving Ground, MD 21005-5066 


l. Introductioa 


Graphical represeatatioa of data is an old, widely used analysis tool. 
Computer algorithms to automatically generate 'X-Y' plots were developed 
quickly by the iadustry. The advent of the low-cost graphic terminal and 
development of miai- and micro- computers led to a tremeadous growth in the 
graphics field during the early seveaties. Graphics, suddealy available at a 
reasonable cost, captured the interest of many users from a host of disci- 
pliaes. Almost instantaneously, computer graphics were being used to display 
three-dimensional objects in 'compucer art', movies, and veadors' promotional 
material. Even the change in language from 'plottiag routines’ to ‘computer 
graphics’ reflects technological advances which allow us to make artistic use 
of colar to show density, stress, temperature, texture, etc. 


A survey of commercially available software was made before this project 
was undercakea la earaest. The leadiag coatender, PATRAN, has a fantastic post 
processor for display of three-dimeasional objects; the graphics modeling 
capability includes surface smoothing and highlighting which produce almost 
Photographic quality results. Unfortunately the company which sells it does 
not plan to market a version which will rua oa any of our three candidate 
computer systems: an НРІ000-Е minicomputer, an HP2836C microcomputer, or a 
РЕЗ252 minicomputer with a MEGATEK 7200 graphics engine. It took very little 
time to learn that buying software desigoed for our particular hardware 
combinations was aext to impossible. 


Since the ultimate goal of this project is to make a movie of the tube 
motion, it was decided to conceatrate development efforts using the Perkia- 
Elmer minicomputer, the MEGATEK, and a MATRIX camera. DISSPLA, version 9, has 
been installed oa the PE and is the curreat graphics package belng used. It 
has the advantage of familiarity білсе It is also installed oa the malasite 
CYBER 173 and CYBER 825 computers. The 3-D graphics available with DISSPLA 
canaot be used for this application since DISSPLA expects а siarle-valued 
function to describe in x aud y the surface. 
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Su The cbject of this paper is to present a simple, direct method for 

* showing gua tube motion as predicted by a gun modeling code, The program has 
beea written so that the user can generate the graphics picture usiag familiar 
two-dimensional plotting routiaes. Further, the program has been designed so 
that only oue subroutine needs to be re-written when the code is 'ported' to 
other computers. With very little effort, the graphics subroutine has been 
converted to run on the HP 1000-Е and the CYBER 825. 


2% Graphic Concepts 


The assumption has been made that gun modeling codes can produce 
centerline coordinates (x,y) for each node at a zivea time t and/or three- 
dimensional information which gives the centerline coordinates (x,y,z) and the 
rotation of point on the surface about each axis for each node at a given 
tine t, The radius of the tube at each node is obtained from a database of 
gun tube descriptions. Surface coordinates are calculated at 25 points around 
the circumfereace of the tube at each node and stored ia a three-dimensional 
matrix by node number, circufereace point number, and axis coordinate. This 
tube matrix is similar to a finite element mesh and can be used to display a 
‘stick figure' representation of the tube, 


In order to use the hardware polygoa fill or software routiaes of a 
graphics package like DISSPLA, the data needs to be presented in a form so 
that coasecutive vertices of а polygon аге defined. А plotting matrix is 
computed using the information in the gun tube matrix just calculated. The 
surface coordinates are arranged so that each group of five describes a 

е polyzon where the Fifth vertex is identical to tue first and is used to close 
the polygoa. The plotting matrix will become much larger than the tube matrix. 


ALL that remains now is to transform this three-dimensional data into 
information which caa be used for a two-dimensional display medium. At the 
present time, a 4 x 4 transformation matrix is calculated which performs an 
isometric projection of the plotting data aad then projects it onto the z=0 
plane. Every polygon is displayed aad filled, including polygoas which are 
not visible ia the picture. This takes a iot of time, but does aot alter the 
final view. Ia order to give the three-dimensional sense, the centerline of 
the tube as weil as the surface circumference points at each node are drawn. 
A picture is taken of the tube at its initial posicion. 


To calculate the plotting matrix for the next time step, it 1s necessary 
to perform the required translation of the ceaterline coordinates and the 
corresponding rotations at each node. This is done using the tube matrix and 
the plotting matrix is thea defined again. 


It is obvious that a hiddea line or hiddea surface routine must be added 
to this program in order to speed up the actual display of data on the 
graphics terminal. Since a gua tube is a series of convex volumes, the 
Roberts hidden line algorithm (reference 1) can be used to compute hidden 
lines. New eadpoiats are computed for partially visible lines and can replace 
the corresponding hidden vertices; thus a new polygon describing the visible 
portion of the original polygon is obtained. Totally invisible polygons, 
defined by the hidden vertices of their hidden edges, are aot plotted. 











A EZ AT VN ОЛЫ" gee. ME T у о, 


BOTAS TU cO cgo сат ЫР Mh ata ev, Mo MU. UM b b S ШЕ. NS ME CC S РЕ ACTUS 


ZI MMERMAN 


Visibility, that is to plot or not to plot, is indicated by putting a one or 
zero in the fourth columa of the plotting matrix. Even though the routine is 
CPU intensive, this should greatly reduce the plottiag time for each frame. 
It must be emphasized that this addition will save time; it will aot change 
the picture at each time step. 


3. Program Details 


The program reads input from units 5, 10, aud 11. Uait 5 is usually the 
default input unit in FORTRAN and associated with the terminal in aa interac- 
tive mode. Units 10 and 11 are attached to data files DTUBE and CUBE, Each 
input required is described below. 


The user is prompted for an enlargement factor (EF) aad the desired tube 
name (ВМАМЕ). Respoase is expected on unit 5, The y coordinate ot the 
centerline (aad eventually the rotations) are multiplied by EF so that the 
movement will be noticeable to the user. RNAM is used to pick the correct 
tube description data from the file DTUBE. 


The 
engineering drawings 


database of tuhe descriptions is created by the user from the 
and recorded ia the file DTUBE. This input file is 


attached te unit 10. Each tube is described in the following manner: 
Line 1: (А20,15) 
DNAME = Name of gun tube, i.e. 105mm ~ M68 
NSEC = Number of sections in iube 
Line 2: (12,1Х,А!,1Х,4 10.4) - repeated NSEC times 


[SEC = Section number 
ТУРЕ = Cross-section description 

R -> rectangular 

T -> trapezoidal 

D -> discontinuous 
ХВ = x coordinate at beginning of section 
YB = y coordinate at beginning of section 
XE = x coordinate at end of section 
YE = y coordinate at end of section 


Figure 1 shows the diagram 
entered tn the file DTUBE. 


of the 105mm gun tube aad Table 1 shows the data as 





Illustration of the l05mm - M68 Gun Tube. 


Figure l. 
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Table 1. Example of data as entered in file DTUBE. 


105mm - M68 6 

LR 0.0000 8.9000 22.5000 8.9000 
2 R 22.5000 8.9000 59,5000 8.9000 
3 T 59.0000 8.9000 103.2500 7 .0100 
á T 103.2500 7.0100 126.8500 6 .2500 
5T 126 .8500 6.2500 202.4375 6.0500 
6 D 202 4375 6.0500 210.5000 6.0500 


The file CTUBE, which is attached to unit 11, is created by the дип 
modeling program. Tt must contain the following iaformation: 


Line 1: (NODES (6X,I5) = Number of nodes chosea (does not change 
during the run) 


Line 2: TIME (6X,E17.10) = time at each step 


Line 3: (6E12.7) -> to be used for а 3-D modeling program 
С(Л,1) = x coordinate on centerline 
C(J ,2) y coordinate on centerline 
C(J,3) z coordinate oa centerline 
R(J,1) = rotation of point on surface about x axis 
R(.1,2) = rotation of point oa surface about y axis 
R(J,3) = rotation of poiat ол surface about z axis 
a where J = nodal point number < 20 and this line 
De is repeated for each nodal poiat. 


! 


H 


Data at subsequent time steps repeat Line 2 once and line 3 CNODES number 
of times. Line 1 is never repeated. Table 2 shows the data generated by the 
Boresi modeling program (reference 2) for the 105mm at time zero. This is a 


2-0 program aad the format of line 3 is (11X,2E18.10) since only the x and y 
coordíaates of the ceatevriine are provided. 


Table 2. Example of data as expected ia file CTUBE from the Boresi model. 
8 CNODES 


.0000000000Е 00 INTTIAL QONDITLONS 
1 .0000000000Е+02 .31579880285-04 


2 .2250000000Е%02 .2731650263Е-03 
3 .4090000000Е+02 -.1280901859Е-03 
4 .5950000000Е+02 -.1519039589Е-02 
5 .1032500000Е+03 -,8859127947Е-02 
6 .1268500000E*03 -.14970799451-0] 
7 .2024375000E*03 -.4300180835Е-0! 
8 .2105000000E*03 -.4623445425Е-01 


The program uses the gun tube description data to find the radius at each 
of the centerline nodal points. Tt has been assumed that the modeling program 
will include the ead points of each tube section in the nodal point set. The 

"1 program will create two points at the beginning of each discontinuous section. 
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Table 3 shows the plotting nodes which have been determined from the tube 
description dota and the gun model data. Note that the y coordinates have 
been multiplied by the enlargement factor ЕЕ=1000. 


Table 3. Plotting nodes determined from DTUBE and CTUBE input. 


Plott inp, 


Node x y z 
l 0.00000 0.03168 0.00000 
2 22 .50000 0.27317 0 „00000 
3 40 .89999 -0.12809 0 .00000 
4 59 £50000 -1.51904 0.00000 
5 103 .25000 -8.85913 0 .00000 
6 126 .85000 -14.97080 0 ,00000 
7 202 .53750 -43.00182 0.00000 
8 202 .53750 -43 .00182 0.00000 
9 210 .50000 -46.23455 0.00000 


Th  (v,z) coordinates оп the surface of the tube at each node x are 
calculated in subroutine CIRCLE. This subroutine uses a parametric equation 
of a circle so that the coordiaates are evenly spaced around the círcumfereace 
of the circle rather than evenly spaced on the y-axis. The results are stored 
ia the tube matrix G(KJ,KC,l) where KJ is the node number, KC is the number of 
the circumference point on the surface, and I tadicates the x, v, or z 
coordiaate. The first and Jasr poiais are identical in order to close ihe 
circle. 


The plottiag array ZP(N,3) is created from the tube matrix G(KJ,KC,3). 
As stated above, the surface ccordinates are arranged in groups of five to 
describe the polygoas which approximate the surface of the tube. The size of 
N is determined by 


5 x (# pts va circlo -1) x (£ plotting nodes -1). 


Assuniag that 19 nodes are specified by the gua modeling program and that 
there is just one discontiauous section, thea the number of plotting nodes is 
20. The program always calculates 25 points on the circle. Therefore, in 
this case, N = 5 x 24 x 19 = 2280 aud 6840 words are required to store the 
entire plottiüg array. Tt is obvious that the size of N grows quickly and 
that care must be exercised to keep the arrays within the memory iimits of the 
computer being used, The current limit for N is 2500. 


Subroutiae GRAPH3D is the driver routine for calculating the homogeneous 
coordinates which will be displayed on a two-dimensional medium. It expects 
data іл aa N x ^ plotting matrix where the first three columns define the x, 
y, aad z coordinates and the fourth columa indicates if the point is to be 
plotted. The subroutine call is 


SUBROUTINE GRAPH3D(ZP , IR, T9 , MH, INITL) 
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where ZP = input data matrix 
IR = number of rows in ZP 
T9 = 4x4 transformation matrix from subroutine TRANS 
ZH = homogeneous coordinate matrix 
INITL = O. 


Normally, options can be selected to describe the type and order of transfor- 
mations to be made to the data. For this application, am isometric projection 
of the data onto the z=0 plane has been preset. In general, the transformation 
matrices T(4,4) are calculated and concatenated in the order specified by the 
user. Тһе crtions supplied by GRAPH3D are used in subroutine TRANS to select 
successive 4 x 4 matrix multiplications until tne final transformation matrix 
1s calculated. These multiplications are done first in order to avoid succes- 
sive matrix multiplications with tne much larger plotting matrix. It is 
important that the options specify the transformations in the correct order. 
The last transformation matrix for a given set of options 


Л 


[T] = іту] х (1-1 x ee XUD] 
is used in GRAPH3D to calculate the transformed data: 


[ZP] x [T] [ZT] 


Nx4 4x4 Nx4* 


All transformations assume a right-handed coordinate system; the user is 
responsible for making any adjustments for left-handed coordinate systems. 


The homogeneous coordinates are determined by dividing the elements of 
each row by the fourth element of thet row: 


ZH(IROW,ICOL) = ZT(IROW,ICOL)/ZT( IROW, 4) 
where ICOL = 1 to 4. 


"us, all fourth column elements now equal one. Since the transformation 
usually includes a projection onto the z = 9 plane, the homogeneous coordinate 
matrix [ZH] contains the two-dimensional representation of the  chree- 
dimensional data in the form [x y O 1]. 


Detailed descriptions of the mathematics used for the translation, 
rotation, and projection matrices used in subroutine GRAPH3D and TRANS can be 
found in reference 3. This is an excellent text for mathematical techniques 
required to develop computer graphic algorithms. 


Now, since the fourth column information has been preserved in the array 
ZP(N,4), it can be inserted in the fourth column of the homogeneous array 


ZH. Finally, the data is ready to be displayed on whatever medium chosen. 
This data can be passed to a subroutine for actual plotting or written to an 
output file to be plotted later. Because of the amount of computation 


required, this prograu writes to file DGRPGP on unit 12. This decision has 
given the fiexibility to try the graphics using different graphics packages on 
other computers without worrying about user program space or memory size. 
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If the user wants to remove uidden lines, this must be done before the 
data is projected onto the 2-0 plane and after all other transformations are 
calculated. An option is available to do just that in GRAPH3D. Subroutine 
ROBERT is the implementation of a hidden line algorithm [reference 1] 
developed by L.G. Roberts at М.Т.Т. in 1963. This algorithm looks to see if 
an edge or line is hidden from view by any object in the scene. Each volume 
or object must be convex in itself; the total scene need not be. Each object 
is described in several ways: 


1. By a volume matrix composed of the plane equations 
for each plane -n the object, 

2. Ву a corner matrix which gives the x, y, and z 
Coordinates of each vertex, and 

3. By a plane/vertex matrix which contains the vertex 
numbers describing each plane. 


The general equation tor a plane 
А.х + By + Ciz + D = 0 
j ง 9 j 


provides the coetticients for each column of the volume matrix 


Aj Ajay te 
[v] = в, Bia oce 
Cj Diane 
Dj Оууу cee 


where j counts the planes. Using the plotting matrix, which is conveniently 
ordered by vertices defining each plane, this volume matrix is determined in 
subroutine VOLMAT. Three vertices on each plane are used to calculate the 
plane equation ccetficients A, B, and C; coefficient D is set to one. Tbe 
following is applied 


“о Yo 45 -1 - B 
m” = ` 
х, Y3 23 1 C 


to each ot the j plunes describing the object. Since an arbitrary 25 points 
were chosen tor the circumterence points, there are 24 planes around the tube 
and one at each end; therefore the volume matrix is 4 x 26. The volume matrix 
must be adjusted so that the dot product of the position vector of a point 
inside the volume cud the volume matrix vields a positive result. ‘The point s 
їз found by taking the average for each X, Y, and Z: 
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y = en * Y, ОЧЕ х Y;)/j 
z = (Z + % + oe. + 2)/j 
and its position vector s = [x y z 1]. If s [V] < 0, thea the corresponding 


columi of the volume matrix is multiplied by -1. 


The corner matrix is calculated in subroutine CORMAT. This matrix is 
simply a reordered subset of the vertices contaiaed in the plotting matrix. 


The plaae/vertex matrix is determiaed ia subroutiae PVMAT. This is 
another bookkeeping type routine which orders the vertex aumbers of each plane 
counterclockwise. The 'sides' of each object, which is the section of the gua 
tube between each node, are planes with 24 vertices. 


As previously stated, the transformations required are concatenated up to 
the potat of projection onto the z=0 plane. The inverse of the transformation 
matrix [T] is stored ia matrix [T9]. The volume matrix is pre-multiplied by 
the inverse transformation matrix to give the volume matrix in the transformed 
space: 


[T] l = [r9] 
[r9] [v] = [ут]. 


The corner matrix is transformed by post-multíplyiag it by the transformetion 
matrix: 


[CN] [7] = [cr]. 
Using the plane/vertex matrix, an edge matrix is computed which gives the 
two vertices describing the edge, the two planes which intersect at the едиз, 


aad а flag indicating visibility of the edge. Initially, the flag is set co 
oae to indicate visibility. 


The determiaatioa of hidden planes 15 made іл subroutine HIDE]. The 
transformed volume matrix is pre-multiplied by the eye point vector 


C = [00 1 0] 
to yield a vector 

G [VT] = PP. 
If each component of the resultiag vector. PP, is positive, then the planes 
are aot seen. The edges formed by the intersection of these planes are also 
invisible. The flag for that edge is then set to zero la the edge matrix for 


the voiume being tested. 


Testing for visible aad non-visible lines begins ia subroutine HIDE2. A 
quick test for complete visibility is done first. Given the two eadpoiats s 
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and ро of аа edge, the liae becweea them is visible is both endpotats are 
visible. The following vectors are calculated: 


d = рә ~ 8 
ру = з [VT] 
. = d T 
q; [УТ] 
. = УТ]. 
v; G [vT] 


1f, for any j, the two conditions 


" ‚ «0 
1 Р, P 


24 Pj 


where G = (0 0 1 0] 


+ q; < 6 and м} <) w 


are satisfied, thea the line is visible. If the edge is visible, its flag ia 
the edge matrix remains set to one and the next edge is processed. If the 
edge is not visible, the flag is set to two aad further testing is done. 


At this point, use is made of a parametric represencation of a line to 
determine its visibility. Given 


P(t) = P| + (PoP) )t where 0 < t x 1 (1) 
V = s + dt 
where 0 <t <l, a >l 


m 


, 
= Py and d = Р-Р 


U 


Qa ,t) = u=s + dt + аб 
Substituting, 


u [УТ] 
u [VT] 


i! 


s [VT] + td [VT] + a СУТ] 


Pj + tq; + aw (2) 


i 


Equation (2) yields a set of j linear equations to be solved for a aad t. 
Each pair of equatioas is solved ia HIDE3. 11 there is a solution for the 
given pair of equations subroutine TLIMIT is called to determine if t is a 


miaimum or a maximum, This subroutine determines the MAX(t ) and MINCt ) 


тіп тах 
for the eatire set of solutions to equation (2). The portion of the line 


which is visible is given for the region of t 
e 2 


< з а < < 
09 st < Liz Vad лах t < | 


aad the Invisible portion [or values of t 


<t $t 


tmin max ° -` 
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The endpoiats are calculated usia t = tnin and t = tax and substituting the 

values da equatton (1). The values of tnin and E dx Are used to determine the 
ง € 

number of visible seymeais from veach ойдо, Invisiblo segments are 


discacded. Visible segments from each edge are stored in a temporary array 
aad are tested адалы! any rematatayz volumes before proceediag to tne next 
ease. Those ѕершеліѕ which are found to be totally visible are stored іл а 
iiag matrix іл subtoutiae H1DE2. 


After ail the edges are checked agaiast all the volumes, except its own, 
the lino segments from che Line matrix aad aay edges which were found to be 
totally visible after the autek check are put {афо an array for plotting. 
Projection onto the 22 plane fr dou just prior to writiag the plotting 
array, 


4.  Conclusioas 


Although calculations have aot beea done with as small а time step as 
will be necessary to produce enough frames for a movie, the method has been 
demonstrated. Тһе hidden Line algorithm has oot been applied to this specific 
problem, bat it has been tested for a scene composed of three rectangular 
paralleipipeds. The routine will have to be modified to handle the 24-sided 
planes of each tube section. 


Oaly the Boresi gun model has been uses for input so far. Thus there has 
been ay test of the votacions; айаіс, the subroutine involved has been used to 
calculate rotations af other obiects ала its results have beon confirmed, 


Simple eashaacemeazs are planned like usin; a contrasting colored stripe 
on the surface of the tube to show twist and a marker of some sort to show the 
positioa of the projectile. Eventually, the motion of the projectile will 
also be displayed. In this zase, the iovisible section of the tube wili be 
drawa with the visible portioa ot the projectile. Highlighting, shading, aad 
surface smoothing are desirable but cot р1алаг! бог the near future. 
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ABSTRACT: 


A study of the projectile-gun tube interactions using a two point force 
technique was presented at the Third US Army Symposium on Gun Dynamics. It 
showed that the large computation-experiment.l discrepancy of one to two 
orders of magnitude for vibrations in the M68 gun tube could be eliminated. 


Tnis pzper presents the work related to the two point force formulztion of 
the general case of the transverse vibration of an elastic gun tube with 
uniform rifling. It outlines the formulations of projectile motion and tube 
vibration, presents correlations of computations and tests, aad discusses the 
theory of formulation. The success of the formulations is in the use of 
elastic deformations to replace the stereomechanical impact for projectile 
motion and the use of a two point force technique to replace the point mass 
approximation for tube vibration. The computations show the important effect 
of projectile c.g. eccentricity and tube-projectile clearance. A detailed 
сотр" tation of the muzzle motion of the M68 gun tube is obtained using a 
typical pressure-time curve. From the computed results it is seen that the 
gun vibrations fluctuate depending on the variations of initial conditions. 
For given realistic initial condition, the computed results are comparable to 
the test data and there exists no computation-test data discrepancy. The wide 
range in the variation of the computed results can be explained bv recalling 
the principles of resonance of vibration. The evaiuation of the computed 
results gives an insight into gun tube vi^rations. 
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SZU HSIUNC CHU, PH.D. 
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LARGE CALIRER WEAPON SYSTEMS LABORATORY 
DOVER, NEW JERSEY 07801 


INTRODUCTION 


The existing state in gun dynamics js that the computation of the 
transverse gun tube vibration generates results that are one to two orders of 
magnitude less than the corresponding experimental data (1, 2)*. However, the 
torsional vibrations of the gun tute have good computation-test agreement (2, 
3). In order to improve the computations, a two point force technique was 
proposed at the Second Meeting of the Gun Dynamics Steering Committee (4). 
This technique considers a projectile as a rigid body except at the bourrelet 
and the rotating band where elastic deformations may occur. The transverse 
projectile-gun tube interactions computed from a six degrees of freedom 
formulation of such projectile motion are then used as the two exciting forces 
with opposite sign at the bourrelet and the rotating band locations of the gun 
tube. Thus, in eddition to the effect of the conventional point mass forces, 
the effect of a monent induced by the base pressure resultant and the 
projectile c.g. eccentricity is aiso included. The M483 projectile 
computations (4) show that tie results are dependert a the projectile c.g. 
eccentricity and agree with the available test data. This technicue together 
with a simple example was later presented at the Third U.S. Army Symposium on 
Gun Dynamics (5). Reference (6) shows the effect of a moving couple 
introduced by the projectile eccentricity and (7) presents a si: degrees of 
freedom formulation. Both use similar concepts, ard indicate tha. the 
technique has many advantages. 


This paper documents in more detail the two point force technique and the 
rclated gun dynamics work performed in the Applied Sciences Division of the 
LCWSL, ARDC. its aim is to present some observations and insights into the 
gun dynamics field. The modification of current formulations using the point 
mass approximation may be easily done by substituting the original forcing 
function with the two transverse interactions at the bourrelet and the 
rotating band, or by an equivalent force and moment. 


This paper wiil emphasize the principles of problem formulation, 
descriptions of forces and deformations, aud correlations of computations and 
tests. The equations of motion wili be mentioned briefly, since the 
derivation is rather a general mathematical manipulation once forces and 
masses are defined. 











*Numbers in brackets ( ) in the text denotes the reference listed at the end 
of the paper. 
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The presentation begins with the formulation of the projectile motion. 
After the motion and forces are determined, the opposite of the transverse 
forces acting on the projectile are used to formulate the transverse vibration 
of an elastic gun tube with uniform ritling. Methods and examples of solution 
are presented. A detailed result of the computation of the muzzle motion of a 
M68 gun tube is obtained using a typical pressure-time curve. Various cases 
of different initial conditions are considered. Some correlations between the 
computations and che field tests are mentioned. The discussions of the 
formulations and the computed results indicate the advantages of this 
formulation and may give some insight into the complicated gun dynamics 
problem. 


FORMULATION OF PROJECTILE MOTION 


The forces acting on a projectile during launch are the base pressure, 
gravity, projectile-tube interaction, air resistance and the inertia force. 
Knowing these foices are essential for designing projectiles and fuzes, 
defining initial conditions of exterior ballistics ard generating tube 
vibrations. 


In 1971, a research program was initiated to determine the forces to be 
used in fuze and projectile design. After surveying related literatures, an 


extensive formulation was performed and documented (8). This formulation has 
the following features: 


1. Six degrees of freedom to replace the ccaventional three degrees of 
Treedom 


2. Physical reasoning is used to define the applied forces instead of 
mathematically assuming the three components of a force; 


3. Introduction of elastic deformation instead of traditional 
stereomechanical impact at the bourrelet contact; 


4. Emphasis of the projectile c.g. eccentricity from its geometrical] 
axis! 


5. Consideration of the effect of non-unifoim band engraving of the 
rotatiug band; 


6. Inciusion of the influence of tube curvature generated by its own 
weight or temperature difference; 


7. Initial position of the projectile; 


8. Using the base pressure directiy as input to the system of motion 
equations; 


9. Including the effect of transverse tube vibration; 


10. Considering rifling effects; aud 


11. Throw off at muzzle. 
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The program has undergone successive improvements and documented in (4, 5, 
9-15). Additional features as follows are considered: 


1. The air resistance or avrodynamic forces at the front portion of the 
projectile; 


2. The friction forces at the bourrelet and the rotating band; 
3. The effect of tube wear and copper deposit іп a worn tube; 


4. The center of the engraved rotating band is not confined to move along 
the tube axis, .nd lateral motion is permitted; and 

5. The effect of the resisting moment of the rotating band, which оррсвез 
the transverse rotation of the projectile. 


The six degre2»s of freedom formulation (8, 11, 15) is based on the basic 
principles of rigid bcdy dynamics, that is, Newton's second law of motion for 
trauslational motions and its complement in th? law of moment of momentum for 
rotational motions. For economic purposes, the products of inertia of the 
projectile is ignored. Consequently the following Newton's and Euler's 
equations of motion are used, 


та, = Fy Ч) 
' Sy (2) 
ша, = F, (3) 
bu - Q oh. + Qh, = My (4) 
hy - Shy + Bh, = М; (5) 
h. - 2h, + Qn, = Mj (6) 


where 
ш = mass of projectile, 


ах, Ayr ay = accelerations in the x, y, 2 directions, 

Pgs Fyr Е, = forces in the xy y, 2 directions, 

hi; hor h, = angular momenta in the 1, 2, 3 directions, 
2%! “2, น ว = anpular velocities iu the 1, 2, 3 directions, 
мү, м. . М, = soments in the 1,2, 2 directions, and 


dot over a quantity devotes it. time derivative. 


in addition, spim ing is governed by the rifling eguation, 
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Wtany = R(j + фсов8) (7) 
where oe 
W = displacement of rotating band in z direction, 
ү = twisting angle of rifling, 
R = bore radius, and 


(565707 Euler's angies. 


To properly define the motion and forces, several sets cf fixed and moving 
right-handed rectangular coordinate systems and Euler's angles are employed. 
These coordinate systems are described in Tables 1 and 2. The fixed X, Y, 2 
coordinate system is the basic reference system and it is fixed to the 
stationary portion of the gun tube cr the ground. The moving X', Y', Z' 
system is parallel to X, Y, Z system but its origin moves with the driving 
band. The moving Xj, Үр, 2) system translates with the rotating Land but 
rotates so that the Z,~axis is always tangent to the yun tube axis, The 
moving coordinate system, 1, 2, 3 and 1', 2', 3' are parailel to each other 
but with different origins. The body-fixed i", 2", 3" system is fixed in the 
projectile and is equivalent to the 1', 2', 3' system rotated an angle фу 
about its 3'-axis. 


The Euler's angles, у, $ and 9 , are defined as follows: 


1. wẹ is the angle of precession, which is the angle between the nutation 
axis (intersection line of the driving band plane and the X-Y plane of system = 
X, Y, Z ) and the X-axis; “- 
2 + is the angle of spin, which is the angle between the body- fixed 1"- 


axis in the rocating band plane and the nutatisn axis, and 
, 


3. 0 is the angle of nutation which is Formed by the Z-axis and the spin 
axis. 


Superscript ' and subseript 1 are uses for Euler's angles with respect to 
the reference coordinate systems X', Y', 2" and Ху, Үр » 4 respectively. 


The forces and moments a:ting on the projectile are basically visualized 
or the assumption that the projectile has elastic defomuatior at the bourrelet 
and tie rotating band, and othezwise rigid everywhere. All forces aud moments 
considered are shown in Fig. 1%. The subscripts used in this figure have the 
following meanings: 


А = aerodynamic force or air resistance, 
а = inertia force due to motion of gun systea, 
DRT = bourrelet normal and frictior force, 





*Fipures are shown at the end of the paper 
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. TABLE 1. 








SYSTEM 


NAME 
X,Y.2 
(fixed) 


X',Y',7' 
(moving) 


Y, Y 
he 12 


(moving) 


1 


SYSTEM 
MAME „с. 


1,2,3 
(moving) 


1 2',3! 
Y 


l', 
(moving) 


| оа 
(body- 
fixed) 


X',Y',Z' system plane axis 
TABLE 2. RIGHT-HANDED CARTESIAN COORDINATE SYSTEMS-SHELL AXES 
FIRST SECOND THIRD s 

_ _ ORIGIN ` AS 71 ฝ ล 7 ล AXIS — AXIS 

At shell C.G Parallel to nuta- Normal to shell Parallel to E 

tion axis 1' axis shell axis 

At moving point Nutation axis in Parallel to 2- Parallel to 3- 

where gun tube axis center plane of axis axis ^ 

intersects rotating rotating band 

band plane. Same 

às origin of X',Y', 

H ' 

Fixed on shell at Intersects 3- Normal to shell Parallel to 

point coinciding axis at negative axis shell axis 

with origin of side осе with 

1',2',3' systen 3J'axis) 


ORIGIN 


Fixed at point where 
gun tube axis inter- 
sects rotating band 
plane at t-0 


Always at moving 
point where gun 
tube axis inter- 
sects rotating 
band plane 


Same as origin of 


Mi FIRST 
MIS _ 





Horizontal, point- 
ing to right when 
viewed facing 
muzzle 


Parallel to Х- 
axis 


Same as X-axis 


ЖЕН; ooo 


SECOND 


Іп vertical 
plane 


Parallel to Y- 
axis 


In vertical 














THIRD ` 


RIGHT-HANDED CARTESIAN COORDINATE SYSTEMS -GUN TUBE AXES. 


Жыл ____ 


CREE 


Tangent to gun 
tube axis 


Parallel 
axis 


io 2- 


Tangent to tube 
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g = gravity force, 

p = base pressure, 55 
КВ = rotating band normal and friction force, 

RF = rifling normal and friction force, and 


x,y,z = components in the x,y,z directions 


The forces at the bourrelet is worthy of special attention. An elastic 
deformation or spring force is introduced instead cf the traditional use of 
rigid body contact or stereomechanical impact force (see Discussion 
section). This technique avoids the jerky computed results. A general 
configuration of bourrelet deformation is determined by the contact condition 
for an unevenly engraved rotating band case. The associated bourrelet force 
is the force resultant acting on the contact area. The simple expression for 
the normal bourrelet force is: 


кб (8) 


с 
il 


normal bourrelet contact force, 


= 
v 
" 


> 
с 
ย 


spring constant, and 
we = deflection at the contact. pru 


The rotating band forces may be obtained by integration around the 
rotating band. To avoid not including the effect of the uniformly distributed 
forces, the formulation computes the uniformly distributed force and an 
unbalanced force due to ihe lateral motion separately. The unbaianced normal 
force has an expression similar to the normal bourrelet force. Furthermore, 
for a rotating band of very narrow width, the resisting moment due to uneven 
force distribution in the tube axis direction may be ignored. In case this 
шөшелі is not neglected, it is expressed by 


M = Cs (9) 


= 
i 


resisting moment, 


1 


yaw anglc, and 


coefficient of resisting moment. 


The air resistance or aerodynamic force at the front portion of the 
projectile is not formulated in the same manner as the base pressure. Rather, 
it is considered as forces acting on an inclined object moving in the air 
confined inside the gun tube. Consequently, equations similar to the 
aerodynamic force formula are used. 


The rifting force is determined from the rifling condition (8). When a - 
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projectile is launched in a smooth pun tube, the situation is quite different 
and the formulation is much simpler (14). 


The equations of motion for both drooped and vibrating gun tube cases 
present additional computation problems over and above the original equations 
programmed for the fixed straight gun tube. These problems were solved by 
referring the motion to the ground fixed coordinate systems and using a number 
of coordinate conversion equations. A simpler set of equations have resulted 
from considering the relative motion with respect to a moving coordinate 
system which translates and rotates as it moves along the gun tube axis. This 
technique computes the projectile c.g. accelaration according to the following 


general acceleration expression and makes corresponding modifications in the 
formulation. 


The general vector expression for the accleration, a, of a point referred 
to a moving coordinate system is from the theory of dynamics of a rigid body, 


a= В + oxtuxp) + фхр + D + 2%>5 (10) 
c r r 


where 


өл 
h 


є = acceleration vector of origin of the moving coordinate system, 


ш rotation vector of the moving coordinate system, 


ล อ บ P e displacement, velocity, acceleration vector of the point 
relative to the moving coordinate system. 


Manipulating the cross product operation and separating the acceleration 
and associated components in the Xi, Уу, Z4 directions, these components of 
acceleration are 


T M . | : "E 2 
a = Bu + XQ - ARA + 27. - MN. (ш to.) 
з — & + 2. to 
bi Y (ox ү, SR 1 "z] X] d 
a = R y. + PEL. 27 ๓ + 2X4 У. fou А, v? ) 
Yi C l l х1 21 1 X1 
n SG FR + ú (12 
+ 21 (iy Mz, x, i x "Y, GC ) 
аша GR 4 HB - 2X, + 211. -Z (5 + „2 ) 
21 CZ] 1 1 Yi 1x4 17х Yi 
Ox. ہس‎ seu ) + V (ш, e, oues ) (13) 
l 21 Za Yi 1 Yi 21 хү 


In the case of a drcoped gun tube which is considered stationary with 
curvature in the vertical Yj7 Z, plane, the motion of the An, Үү, 2 
coordinate system is coufined in this plane and hence its acceleration 
- components are only Boy and VM . and the angular velocity and acceleration 
| components are only uy and Ex]. Consequently, the acceleration 
compoaen s of the mass center of the projectile reduce to: 
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a = 1 
xi i (14) E 
= 3 ¥ 22 aY wl - gu (15) SS 
зүү Rey) + У} - 2196 1” ey 1 ху 
5 T ; 2 я 
32 : ล t төх, SR 
1 


The moving coordinate translation and rotation terms in these equations 
are detcimined by differentiating the equation of the tube axis with respect 
to time and the travel of the projectile. The curvature of the gun tube axis 
may be experimentally determined or computed. A simpler method is as follows: 


The non-uniform gun tube is divided into many segments or stations, 
counted from the fixed end. The deflection and slope of the gun at 
station n are determined from the equations for a general cantilever beam, 











namely (25), 2 3 
oy 5 Май п + Solo 
1 ม eet ж ฮิ ก -1 ไ ใก + 2EI 3EI (17) 
^n n 
4222 
4 , " Matn Sn*n 
= + —— 
n n-l ELp 2EI, (18) 
where 
Y = deflection, 
9 = slope, š 
М = moment, = = 
< 
2 = length of beam element, 
I = section area moment of inertia, 
S = shear force, and 
n = subscript to denote station. 
In the case of a transversly vibrating gun tube the complete Equations 
(11), (12) and (13) are used. These equations show that the transverse and 
the torsional vibrations are coupled, since they contain terms with 921 
and 9:24, which are due to the rotation of gun tube about its axis or 
torsion. When torsion is not considered, these terms are omitted, and the 
acceleration components become: 
Š 52 В Y 27.1 - „2 + š ) + m 
day М, pn d^ 6 7 түзүлүү ae 7 (19) 
. 2 . 
= р үз. oae Эў - Y4u - Zio + Xj wy W 
ay Roy) + Yi сроку 1 Xj l X, 1 X1 Yi (20) 
ล = Rezy + Z. = 2X ш + 2¥ دا‎ - 2 (๒ 9 + wë ) - хо + ND (21) 
21 1 1 Yq 1 хү l хі 1 1 Y1 1 x) 


These equations contain coupling terms “x, and шү and their time 
derivatives. This means that th. transverse vibrationà in the vertical and 
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horizontal planes wil) not be independent and the usual separate computation 
technique is indeed only an approximation. The moving coordinate system 
translation and retation terms are determined by the gun tube vibration which 
is described in the next section. 


At the same time Picatinny Arsenal started this work in 1971, similar 
investigations were also performed at the Sandia Laboratories, Albuquerque, NM 
(16, 17, 18), later at Naval Surface Weapons Center, Dahlgran Laboratory, VA 
(20, 21), and recently at S&D Dynamics, Inc., Huntington, NY (7). 


FORMULATION OF TUBE VIBRATION 





In most transverse tube vibration studies, the projectile is considered to 
be a point mass since it is small compared to the gun tube, and the gun tube 
is approximated as a cantilever beam since it is a long tube with the breech 
end supported. In reality the transverse vibration of the gun tube is not 
confined to a plane. However, to simplifiy the formulation, it is usually 
considered that the tube vibration has components in two perpendicular planes 
and may be analyzed separately. In general, Euler's or Timoshenko's beam 
theory is used in the formulation of the equation of vibration. The general 
equation derived from Euler's theory for the X-Y (vertical) plane vibration is 


(EIY')" + my = s (X, - ХЕ (22) 
P P 
where 
m = mass of the gun tube per unit length, 
E = Young's modulus, 
I = area moment of inertia of the gun tube element section, 
X = location of the gun tube element, 


X, = location of the projectile c.g., 
Y = normal or Y-axis displacement of the gun tube element, 


Е = projectile-gun tube interaction force exerted on the gun tube by 
the projectile, and 


ê (Xp ~ X) = Dirac delta function. 


The gravity force of the beam is usually i,nored in the formulation. A 
similar equation is used for the X-Z (horizontal) plane vibration. 


The interaction force F, may be derived from the point mass approximation 
(22). An important point is that this formulation ignored tne effect of a 
moment which is the product of the base pressure resultant and the projectile 
с.р. eccentricity (see Discussion section). This moment is not small for high 
base pressure cases. Consequently, there is no satisfactory agreement between 
the computed results and the experimental data (1, 2, 3, 23). Nevertheless, 
the computation of tocsional tube vibration did show good agreement with the 
test data (3). This difference of agreement is perplexing since all theories 
and techniques of dynamics are well developed and rvailable. 
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To improve the focmulation, the technique or concept of a two point force 
approach is introduced (&, 5, 15). The equation of vibration becomes 


5 . 23) 
ПУ") "+ mY = - ~ 2 - XN ~- (X +h - X)N ( 
(ELY") l E S ( х ) by 
where 
ม = axial distar:e of the projectile c.g. to the rotating band 
ceuter section, 
h = axial distance of the projectile c.g. to the bourrelet, 
Noy = Y component of normal force acting at the rotating band, and 
му = Y component of normal force acting at the bourrelet. 
The forces -N and -N are the components of the two projectile-gun tube 


oy 
interactions acting at the rotating band and the bourrelet location 


respectively. They are computed from the projectile motion but with opposite 

sign. They may also be represented by an equivalent force and a moment at the 
projectile :.g., since the equivalent of two parallel forces is a force and a 

moment. Нсчеуег, the two point force representation is preferred, since they 

represent completely the force actions until the rotating band leaves the 

muzz le. 


The acceleration and rotation of the gun tube are obtained by 
differentiating the displacement and the slope respectively. These quantities 
are used in the equations of motion of the projectile as mentioned before. 


METHOD OF SOLUTION 


Traditionally the computations of the projectile and tube vibration are 
performed separately. Usually, the equations of motion of the projectile are 
solved fir.t to obtain the displacement, velocity and acceleration data 
without considering the effect of tube vibration. The computed projectile 
data are then used as known values to solve the equation of tube vibration. 

Yn reality, there is a coupling cffect between the projectile and the tube. 
This may be seen from the appearance of the coupling terms such as interaction 
forces in the vibration equation, and tube accelerations and rotation terms in 
the equations of motion of the projectile. Therefore, for accurate results, 
all these equations should be solved simultaneously. 


The computation time of solving simultaueous equations increases with the 
number of equations involved. To utilize the existing projectile analysis 
program and save some computation time, an alternate numerical integration 
technique is used. This approach solves the equations of motion of the 
projectile and the tube vibration alternately at each time step of numerical 
solution. The only data required from the interior ballistics computation or 
test is the base or chamber pressure. At the first time step the projectile 
equations are solved with the gun tube at rest or no vibration. At the second 
time step, the result of the projectile solution at the previous time step is 
used as the known input data to solve the vibration equations of the tube. 
The results of this solution is then used to solve the projectile equation. 
The same procedure is repeated for each time step until the projectil- 

rotating band is out of muzzle. This technique has the advantage of i th 
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simplifying the computation and computing the effect of projectile-tube 
interaction since the time step of integration ia taken small. 


Any of the standard numerical integrati on methods may be used in the 
solution. A simpler technique is the Euler's method and the constant 
acceleration technique (24). The basic equations used are 


d ER д 
Ynil Ya y At (24) 
2 
= + <ç + ыу = Р 
Yol Ze Tas Y. (25) 
where 


ЗЕ = time step size, and 
n, ntl = subscript to denote the time step number. 


Modal analysis is used to solve the equation of motion of the gun tube. 
The пабитз1 frequencies aud normal functions are solved using the Myklestad's 
technique. These techniques are well doeumented in many text books and 
technical reports, such as Reference 25, and not repeated here. These 
techniques transform the tube vibration equation into the following equations: 


n 
ү = $ Y,(X)q, (t) 


iri i (26) 
Ж 2 _ _i : ‘о 
di + =ч, =- gr Noy¥io" NpyYiy) (27) 
n 2 E Y 
mee Ap (28) 
Hj = ز1ز"‎ 
Where 
Y; = ith normal function or mode, 
q; = ith normal coordinate or moda: response, 
w; = ith natural frequency, 


m. = jth lumped mass of the gun tube, 
M. = ith generalized mass, 


Yip = ith normal function or mode at the rotating band, 
the bourrelet contact point Location, and 


n = total number of lumped mass of the gun tube. 


Computer programs have been generated based on the above mentioned 
technique to solve the following examples of computations. 


EXAMPLES OF COMPUTATIONS 


The input data for the computations are the geometrical dimensions, 
physical properties, initial conditions of the projectile and the gun tube, 
and the firing pressure. They are obtained from design data, test results or 
interior ballistic computations of propellant charges. The equations of 
motion are solved by standard numerical integration methods. 
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The solution of the projectile motien can calculate all comoonents of the 
moticn and forces. Some of chese teras avet projectile travel, velocity and : 
acceleraticn; c.g. and bourrelet center uisplacement (polar diagram), T 
bourrelet contact point and deflection; absolute i;teral velocity of c.g., 
normal accelerations of c.g. sud roce specified axia; points; yaw angle, 
velocity and acceleration; cioss spin rate; force coupouents at bourrejet and 
rotating band; tota! forces at c.p; etc. These computations could not be 
vbtained with previous formulations. The computed values fron tnis 
formulation are used 25 design reference values fur projectiles and iures. 
Fxampies cf nroje.tiles thet have been computed ere the XM673, 41712, M39242, 
X4829, М107, M453, M549 and so forth. Some typica! curves of the computed 
values are shown on Fig. 2. 


In the computation of gun tube vibration only fixed breech end cases are 
considered at present. The input data include the outside diameters, area 
moment of inertia at different elem mt stations along the tuhe in aecitica to 
those required for projectile analysis. Тһе computes results represent the 
tube configuration at any time, that is displacement, velocity and 
acceleration of differeat stat con peints, Special attention is paid to the 
muzzle motion and its displacement, slope, velocity and acceleration are 
computed. Some of the computed results of the 120 ะ ห ท snd M68 105mm guns асе 
shown in Fig. 3-5 and Table 5. 


CORFELATIOR OF FIELD TESTS 
Most of the cest data used for cases mentioned in this paper are from 
design drawings, similar items, sad simplificacions since no actual test data 
are available. Even in cases wücn test data are available, many parawecers or & 
constants required in the formulation are not recorded airing the rest. 
Therefore, exact comparison between the computations and the tests are not 
feasible. However, some correlations may be mentioned to show the advantages 
of tais formulation- 


During 1974, a thorough flight evaluatioa was conducted on the 195ша M483 

pou ple fired in tne MLSS gun to determine the cause of short rounds 
previously eacountered in cold weather test at Nicolet, Canada. la order to 
determine the importance of interior ballistic factors cu this progran; the 
bsllot ue motion of the projectile was analyzed using the method described in 
Lhis paper- Іптегіог ballistic performance from various M185 tubes was 
calculated, including a new tube and tubes with a significant amount ot 

wear. It is noted that tube no. 22530 (11) had a restriction in its center 
portion (180" to 75" from muzzle), due to the deposit of copper, which tended 
to reduce the in-bore clearance between projectile and tube Due to this 
croudition, the frequency of impacts at the bourrelet Lended Lo increase 
comared to the cther tubes of better wear conditions and the magnitude of the 
contact force at the bourrelet was significantly higher. A comparison of the 
calculated balloting behavicr with test firings is indicated in the {following 
Table 4. It is noted that the coppered tube (Tube No. 22520) produced a 
higher cross spin rate and higher bourvelet force than a tube with 
siguificantly higher wear in both the origin of rifling aud the muzzle, and 
thac rhe calculated performance corresponded with the higher first maximum yaw 
and higher incidence of short rounds noted in the field firings. After tiring is 
an additional 800 rounds, howevec, the star paging of the tube indicated that 
the reduced diameter section wa removed and perforuance appeared to improve. 
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TABLE 4. EFFECT ОЕ GUA TURE 
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Tube 

No. Type 

22535 Neu 

22539 Ceppeced 

22530 BOO rds 
latec 

782.11 Condanned 


An additional 


= А Corp 
Peak Bourrelec 
Impact Porce 


{lbs.) 
1100 
8900 


5700 


$100 


Instance 


CONDITION - M185 GUN 























prelictions aud the field experience was found. 


the field tests at Yuma Proving Ground 
balance, and the position and degree of unbalance recorde?. 
for initial values cf ths precession an,le, V , а 
180°, 270°; and 270°, 180° respectively, the center to gravity of an 
ition in the gun tube. 


90°, 0°; 


eccentric projectile is in the 12 o'clock pos 





utatfíor о Жук т; 2 Field Experience 
Cross Spin Rate ist Мак. Short Rd. 
at Muzzle Yav, Avg. Rate 
(rad/sec) (deg.) 
«11 1.6 1/20 
1.65 4.7 4/20 
3:19: Not Taken 3/60 
1.01 2.0 0/1C 





of correspondence betweea the analyticai 


The 1482 projectile used in 


were all tested for projectile mass 


It is noted that 
nd spin ange, , of 09, 909; 


Average 


calculated values for initisl conditions with the heavy shell side lecated in 


the 
the 
sole а. Fam - 
Wea D 


ТАМА 5. 


- — 


Venter of 
Gravity 
Orientation 
(o'clock) 


A 


6 
9 


12 


3, 6, 9 and 12 o'tleck positions are piven in Table 5. 
initial orientation calculated to yieid a higber value of cross spin rate 


Computed Peak 


Bourrelet 
Force 
(lbs) 


4500 
4700 
3500 


6200 


№685 PROJECTILE 


It 1s noted that 


CO APO POE а higher value of first maximum yaw in the field firings. 


EFFECT OF INITIAL PROJECTILZ CENTER OF GRAVITY ORTH NYATIGN 


Computed 

Cross Spina Average 156 
Tee Мах Yaw 
(rad/scc) (deg.) 


1.08 


.96 








The 1551 projectiles М107, MARIAN, and M49, fired in the M18 howitzer, 
yy Əspeclaliy iu the Lourcelet engravings found in field 


perform differen} 


tests. 


(actual 
in oz) were analyz 


-omputed results a 


variations of the yaw angle and velocity. 
the normal accelerations av two axial points (6.35 сш and 25.4 cm 


spin rate, 


to the projectile nosey, the lateral terces астар at 


ed. Both new 


and «orn yun tube 


An analysis was рессотаей to ideacify and compare thei. iu-bore 
lateral motion cheracteristics (12). There proj.citles with en assumed с.р. 


data not available} unbalance of l4lmu-N (20 1и-07} and 353mmn-N (50 


s were considerel. The 


га presented in plots and tablas tor the peak values and 
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. bourrelet (Fig. 2C) and the c.g., and the paths of the c.g. and bourrelet 
| center. These results show thac the M549 projectile has the fargest balloting 1754 
i affect of the three projectiles, and agree with the field cest results that 
the M549 rounds have more severe ептгауїлд at the bourrelet. 
All analyses of projectiles with large c.g. eccentricity fired in a cow | 
rifled gun tube, show that the bourrelet aad the rotating band forces increase ^ 
toward the muzzle (Fig. 28). ‘The tube wear at the muzzle end region are due a 
to the friction between the projeccile and tube. Consequently, with the М 
rotating band and bourrelet forces increased, the tube wrar vil! be increased | 
accordingly. This corresponis wall to the general observation of more tube 
wear near the muzzle. 
A computation in 1977 of the M344Al projectile in the 106mm M206 cannon 
showed the bourrelet and the rotating band :orces has а maximum value as shown ы 
in Fig. 2A. Later, а gun tube failed at the location where the computed force 
is maximun- к 
` 3 
Another interesting correlation шау be mentioned here. As show by à 
computation results, the bourrelet force (Figs. 2B and 2C? is inzreased with x 
c.g. eccentricity, and travel along the rifling. Thevefore there are manv m 
pulses of bourrelet force acting along the rifling during Laach, Lf the X 
force pulse is large enough there vili be many dents or local bending aiong S 
the rifling. This may explaic the test fat that Local bending was observed к 
in the vicinity of the projectile (рр 14 of (27)), and these deflections аге E. 
| ео ызаа pear | 
Le "RE | | E 
Tube vibration is of primary concern 1n gun tube dynamics. The large i 
discrepancy between the computed results and tre experimental measurements Es 
puzz led investigators for a long time. Tue computations based оз the point M 
mass approximation are much less than the test data (1, 23). However, the З , 
computations from the two point force technique are well ccmparable to the м 
test data depending on the parameters of c.g. eccentricity, projectile-tuhe £ 
clearance, spring constants and so forth. For example, the BRL test data (26) ^ 
of mean vertical muzzle deflection of the MIOSAL howitzer at shot exit is 75 
2.029mm and the computed result with several c.f. eccentricity values a А 
shown in Table 6. The average salue cf tbe four computations is 1.972run. " 
TABLE ©. COMPUTED MUZZLE VERIICAL DISP! ENIS OF 15548, NIO)AL ү, 
HOWITZER WITH май ROUND ТЕТ 
“базе | € ¿pu FAC ` Cav, Imb. T Uy SSES 80 у 
No. mn wm Ñ Defi. mn 
А (in) (in ะ a2 ` (in) 
l 154 97 авт 2 ื 0 ๐ | | 
( oun? 0) (.04) Ы 
2 IER 84.23 1.54 Ж 
(.001) (19) ( .05) E 
3 ‚202 91.80 2.03 Ba 
(.0.:3) (13) (.n8) ; 
| 4 .231 105.92 1.32 W 
(.009) (15) (14) ۳ 
NOTE: Projectile Launched with zone 3 charge, c.g. at 5 эссе; сөгітісеа a 
veri А 
5 
UT 
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The mos: ฉ่ อ ร ว ชา ไร น์ and compuced „азе Ха Yil2crature is the transverse gun 
vibration of the 468 105та van with а М262А2 piesecctile. test datz and 
zomat tions cre weil reported (23). A detailed computation on the Timed 
breech end case using the metho? mentioned іс 4828 paper sie performed to 
compare the results. The peak rest cata for the vertical esd the z:rizontal 
muzzle displicemants аге 2.25 anl .25ma resucctively ‘Test No. ident 06 of 
(23)). The input data used for analyses ате fron uesige d rwings, similar 
items, and simalifications. Some ai tho computed results are shown ir Table 3 
ard also shown jn Figures 4 to 5. From the gable ard figures it 75 seen that 
depending on the bourrelet contact conditioas, 2.0» eccentricity, sv. ing 
constants acd initial projectile positions, the computed peal mite 
displacemer.t values are either larger or smaller than th: test dats. Whereas 
the computations of point mass approximation always give amaller ง น เน ค ร . This 
shows that the method mentioned in Lhis paper has the capability to compute 
the comparable values, The problex is how to obcain prop2r inp. data af 
computation. 
DISCUSS IONS 

А major difficulty in solving gun dynamics p-coiczs is Tov to define che 
torce relations among the various components of the wen syscem avi how they 
behave during firing. 


The point mass approximation his the advantage of simplirying the 
formulation and the subsequent solution of the resulting equations. However, 
in considering an object as a point mass, sowe chacacteristics uf the object 
may be inadvertentiy iguored. In gur dynamics, the projec.ile is small 
compared то the gun tube and the yaw notion is severely vest-icted by the 
tube-bourrelet clearance, апі Белге traoditionoliy ths point mass approximation 
is used. This leads to two discdventages in the solutinn. 

The first adverse effect 1. in che formulation of the uation of the 
projectile. Tiadicionally, а project sie is first considered ^s a point mass 
to compule the forward ассеіесаії са. To compensate for spin effect, an 
equivalent (heavier) mas. is used. Later, the centritugal force is added to 
account for the c.g. eccentitcity effect.  Furtheruore, čo compute toc yaw 
effect, it is considered as ว spinning body with а yas motion. The point mass 
approximation excludes body deformation. This Irado to Еге solid body contact 
force or the stereouzecbanical impact formulation io deseribe the bourrelet and 
rotating forces. The solid Lody couract force formulation computes а less 
varied or fixed yaw angle while the stareemectanical impact approach computes 
резку motlor The stereomecharical iwpact formoiation is incapable of 
describing the transieut Өмеез or deformations proaucad, and is limited to a 
specificatton of initial and termival velocity states ot tae objects and the 
applied liazeas or angular impulss (28). Furthermore, the pun is tired in a 
rather Short time duration, and thug not including the transient deformation 
leads to ¿naccurace computation. Satista.tury computation: sre therefore 
wunebtainible. The тот. ої mentioned in tuis piper uses the six degrees of 
freedom foraulat:on aad considers the prejectile as а rigid body except at the 
bourrelet auc tne vocatitg bani where deformetion exists, and consequently 
introduces the elastic contact or spring forces ได ง ขิ อ ร อ of the 
atereomechanical โน หง ๆ ท ย Torres.,  Yho successtul cempatations mentioned in 
previous ร อ บ ชู โก ท ย and the fart taat the s.ereemechanical impact formulations 
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are subsequently replaced by similar elastic forces by others in rhe field nu 
Ех Show that this technique of formulation is practical and useful. Vowever, the ki, 
ks small time step size of numerical integration is essential to the ง ў 
computation. Unless some check measure is employed in the computer program to А 
make sure the time step size is small enough, the corputation will be AV 
incorrect. ЖУ, 
Next, the poin mass approximation has a great iniluence ох the transverse a 

tube vibration formulation. Traditionally the projectile is always considered 

as a point mass in the problem formulation since, as mentioned above, the ns 
projectile is apparently very small compared to the gun tube and there is 52 

serious restrictions of yaw motion. The approximation does have the advantage 3 
of simplifying the formulation by reducing the number of equations of à 3 


motion. However, the tube-projectile interactions or the transverse forces 
exerted on the gun tube by the projectile are only one force resultant. 
Conssquentiy, the moment effect which is the product of the base pressure Š 


resultant and the projeccile c.g. eccentricity, is inadvertently ignored. น > 
Unfortunately, this momcat is so large that it should not be ignored for large e 
guns or in cases when high base pressure is used. An approximate computation К 
of this neglected moment of several projectiles 3, shown in Table 7. » S 
% x 


TARLE 7. MOMENTS DUE TO PROJECTILE C.G. ECCENTRICITY AND BASE PRESSURE 














^ Projectile | Time, ms Base Pressure, psi Moment, in-lb 
€) ( MPa) Cm-kg) 
! "XEM ES ЕРТЕН an «fle, Bdge NEEDED ا‎ E 
M106 4.9 37300 18590 
(257) (214) 5 
E 
м107 4.9 35 350 10240 NS 
(244) (118) х TA 
M39242 3.9 49500 4480 m 
(342) Gn к 
M456 2.2 54140 7230 ele 
(373) (85) De 
XM5459 5.0 38250 11076 š 
(264) (127) y 
Sis. Gee cmd ร ร ERU = ง ล a, Se 7 } 
D e ๑ . . D Ы 
NOTE: AL] moments computed with Ol in c.g. eccentricity and firing SN 
pressure from Heppner, Leo D., "Methodology Investigation on ^s 
Setback and Spin for Artillery Mortar, Recoilless Rifle, and 
Tank Ammunition’, TECOM Proj. No. 9-CO-011-075, Repor. No. 
АРС-МТ-4503, Materiel Testing Directorate, АРС, MD, Sept. 1974. 
ea The tube vibration depends on the forcing function Wi.“ lesa exciting 
force the vibration decreases accordingly. This is why all computations using m 
the point mass approximation gave results of much less value than the ps 
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experimental data (1, 2, 23). 





The coucept of noint mass approximation 15 
wibration investigations. This is seen by the fact thar though the rotating 
baud and bourrelev forces have been computed using the six degrees of freedom 
formulazioa, yet these Forces are not used as the two interaction forces to 
excite the tube vibration. Inscead, a single force which is their resultant 
is used. 


Tables 3 aad 6, an’ Figs. 4 ami 5 show that the computation results of gun 
vibrations fluctuate depending on the variations of initial conditions. For 
certain initial conditions, the computed result may be larger, less than or 
match the test data, and there exists io computatiou-test data discrepancy. 
The wide талде of variation of the computed results is not difficult to 
explain if oae recalls the principle of resonance of vibration. The balloting 
motion, spinning, rifling reaction and different natural frequencies of 
vibration пау cause many variations. This insight into gan vibrations may 
help to pin-point important conditions for good weapon system design. The 
successful computations using the method described in this paper and tae 
apparent importance of the initial conditions indicate that this method is a 
proper way to perSue more accurate computations. The problem now is how to 
o^ tain more accurate data to be used for comprtation.s and perform поте tests 
to get ectual data to refine the detalls of the foimulatioa. 


Many coerficients or constaats used in the method cre assumed to be 
obtained, before computation, from experimental measurements, зім: 1аг items or 
may be cemputed using general elastic:iy equations oc finite element 
programs. Either average value or ta ulated values ma; be used. The method ue 
of computing coefficients at each time step is not used and thase constants 
are saved once they have been computed previousty. This reduces computing 
time. 


The resisting monent at che rotation band which sota against the 
transverse turning of the projectile may be outaiued by either test or 
computations. Since the yaw angle is rather small, this moaert is ignored in 
rhe computation by using а zeo coefficient of тезтекаиг moment. This 
onmittance шау have rather srail effect when the width of the rotatinp band is 
small compared to Lhe bourrelet-rotating band distance. 


In the formulation of the aerodynamic forc^s or air resistance, which acts 
ou the front portion of the projectile, a techaique similar to that of 
aerodynamics is used instead of computing a uniformly distributed sir pressure 
on a flat surface. The air is compressed thraugh the vueven tube-bourrelet 
clearance and acts oa the cylindrica’ portion ef the projectile, in addition 
to acting oa the nose portion. ‘he projectile is not always parallel to the 
tube axis, cr there is a yaw motion. Thus the cir pressure з net 
necessarily uniformly applied on the [ront postien of the projectile. Iu 
theory the total force should be ovtained from the integralion of air pressure 
on ali sucfaces of this feont portion. This procedure ig similar to computing 


y 


cerodyuanic forces. The problem now is how to determine experimentally the 
a па data available, 
and some assumed values based or test data of solid objects moving iu air are = 


related coefficients and parameters. AL present, there i 


used, or tae ferc:s ате ignored by using aero coefficients. 
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From all computations, it is seen that the projectile c.g. eccentricity 
SM and bourrelet contact are essential to the projectile motion and the tube 
SE vibration. The moment induced by the base pressure resultant or setback force 
at the projectile c.g., depends on the location of the projectile c.2. with 
respect to the projectile itself and the tube axis. The spring constants, 
tube-bourrelet clearance, engraving of rotating band, etc., further affect 
this state. To reduce the projectile balloting motion and tube vibrations it 
is essential to decrease c.g. eccentricity, and the tube-bourrelet clearance. 


CONCLUS ION 


Successful computations in gun dynamics depends on the proper formulations 
of forces acting at the projectile, and the description of the projectile-gun 
tube irteraction forces which excite the tube vibration. The intrcduction of ; 
an elastic deformation force at the bourrelet and the rotating hand instead of 
the traditional stereomechanical impact make the compuiation practical and 
accurate. Projectile in-bore motions, such as thet of M483, XM673, XM712, 
M107, erc., have been computed. There are many correlations between the 
computatlons and the field tests. 


The introduction of the two point force technique to replace the point 
mass approximation includes the important effect of moment which is the 
product of the base pressurc resultant and thc projectile c.g. eccentricity. 
Thís increases the exciting force of tube vibration aud thus the computed 
results are no longer small as computed with the traditional point mass 
approximation. Thus, this technique can eliminate the computation-test 
discrepancy which has long puzzled scie tists and engineers in the field. ‘The 

Up computations of the MIUJA gun with the 2483 projectile and the M68 gun with 
the ™392A2 projec cile both show that the computed results can exceed, he less 
than or match the test data, depending oo the variations of parameters, such 
as spring constants, projectilc-guu tube clearance, etc., and especially the 
projectile c.g. ecceniricity. 


From all computations, it is seen that the magnitude of the projectile 
c.g. eccentricity and tube-projectile clearame have a large effect oa the 
projectile balloting motion and the tube vibration. The тошеп induced by the 
base pressure resultant or setback force at the projectile c.g. depends on the 
locatioa of the projectile c.g. with respect to the projectile itself and the 
tube axis. The spring coustaats, tube-cbourrelet cloavance, and the engraving 
of rotating band further affect the condition. То reduce the projectile 


balloting and tube vibrations, it is essential to dezrcase the effect of these ài 
parameters, especially the c.g. eccentricity and tube-projectile clearance. 2 
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B ABSTRACT: 


Traditionally, the terms lightweight and long renge, when applied 
j to artillery weapons, have been mutually exclusive goals. To obtain a 
sipnitdcant increase іп range, the impulse applied to the weapon structure 
must also increase. This increased impulse negatively atfects the stability 
ol the piece. Compounding the stability problem is the decreased weapon 
mass due to the lightweight system requirement. Significant weight reduction can 
be aceomplished through the use of lightweight materials (i.e. composites), 


à a reduction in design safety таспаға, or a combination of both. The load 
appiied to the weapon structure becomes critical as safety factors and 

t State-ot-thoe-arct. Lightweight materials are utilized. 

A Traditional recoil design involves throttling hydraulic oil through 
ауасына 1107(( ірің orrides 1 ol reeoti stroke and is 

4 Le mechanically а Гіхес system. In essence, the eritice profile is pre- 


presrazned amd designed around ideal parameters which rarely exist. The 
1 А 


result is a пол-еріішілес cnergy dissipation system which cannot be toler- 
ated dn a Livhitwei;sht Weapon, 


this paper presents the ceneept of a microprocessor centre] led 


А worvovalve. whieh, threush closed loop reedback control principles will 
ensure consistent operation., Iwo Control alsṣoritbms are presented and their 
1 арен te recoil mechanism design are discussed. Finally, the concept 
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reveil and "eett" recoil, auamented with microprocessor 
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LIGHT ARTILLERY RECOIL MECHANISMS 


STEPHEN G. FLOROFF 
NORMAN T LIONETTi 
U.S.AKMY ARMAMENT RESEARCH AND DEVELOPMENT CENTER 
LARGE CALIBER WEAPONS SYSTEMS LABORATORY 
WEAPONS DIVISION 
DOVER, N.J. 07801-5001 


1 introduction 


Traditionally, the terms lightweight aud long range, when 
applied to artillery weapons, have been mutually exclusive . In 
order to obtain long range capability, the impulse applied to 
the projectile must be large. Due to the realities of 
conservation of momentum, this same impulse is also imposed upon 
the supporting structure. This applied impulse has a negative 
aftect upon stability.  Compounding the stability problem is a 
decrease in weapon mass resulting from lightweight design 
requirements. Reduced tiring loads are mandatory if truly 
lightweight artillery with acceptable range characteristics is to 
be realized. Consequently novel approaches to recoil energy 
Managem E are necessary. 


Th. numerous design considerations essential to attaining a 
lightweizht, long range weapon are outlined in the following 
sections of this paper.In sectíon 2, performance characteristics 
for a lightweight 155mm howitzer are discussed and methods of 
obtaining these characteristics are described. In section 3, 
basic recoil mechanism design is reviewed and shortcomings of 
present systems are presented. In section 4, the concept of 
electronic feedback recoil control is introduced. Finally, in 
section 5, two approaches to lightweight recoil mechanism design, 
augmented with electronic control, are described. 


2. Design Approach for a Lightweight Howitzer 


The Army Research and Development Center, Dover, NJ, is 
developing a full scale lightgweight, 155mm towed howitzer 


teasibility demonstrator. The weapon must have the same 
performance characterisitcs as the M198 towed howitzer but weigh 
9,000 pounds. Seiected performance criteria of the M198 are 


listed below. 
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M198 Towed Howitzer 





Total weight - 16,000 1b 
Maximum breech force - 1,500,000 1b 
Maximum applied impulse- 11,400 lb-sec 
(with muzzle brake) 

Weight of recoiling parts- 7,000 1b 
Maximum trunnion force- 80,000 1b 


Thus, overall weight reduction of approximately 44 percent is the 
desired goal. The use of state-of-the-art composite materials 
including glass/epoxy, graphite/epoxy, Kevlar and metal matrix 
are being investigated as alternatives to more traditional 
homogeneous metal construction. Further weight savings can be 
realized by reducing structural safety factors, Through the use 
of finite element analysis, critical areas of the weapon can be 
determined and maximum stress levels can be established. This 
analysis presupposes a specific input force; this force being 
generated through the recoil mecharism. If structural safety 
factors are reduced it is imperative that the applied load be 
consistent. Unfortunately traditional recoil mechanism design 
cannot be relied upon to provide consistent applied force. 


3. Recoil Mechanism Design 





Conventional large caliber artillery recoil mechanisms are 
comprised of three basic components; a recoil brake, a 
counterre2oi] mechanism, and a counterrecoil buffer. The recoil 
brake provides controlled resistance to weapon recoil by 
throttling hydraulic fluid through a variable orifice. The 
couaterrecoil mechanism, or "recuperator" returns the recoiling 
parts to the initial tiring position by storing and releasing а 
portíon of the recoil energy. The counterrecoil buffer reduces 
counterrecoil velocity of the moving parts to zero through a 
hydraulic fluid throttling process similar to the recoil brake. 


Ideally the recoil brake should throttle hydraulic oil such 
tnat a rectangular retarding force verses recoil distance is 
obta.ned. See figure l. Since the area under this curve 
represents the work necessary to stop the recoiling mass , a 
rectangular curve wili yield the lowest retarding force tor a 
given recoil length. Theoretically, the only limit to this 
reduction in retarding force is the physical constraints 
associated with length of recoil. These constraints are 
attributabie primarily to weapon configuration. 
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Figure 1 -"Ideal" Force Versus Stroke Relationship 


The M198 towed howitzer, as an example, can accomodate the 
following recoil lengths 


Propelling Charge Nominal Recoil Length (in) 


Designation 







โว ว ว จุ อ อ 


25 
M4 45 
M119 62 
(reduced to 50 for high elevation) oo 
“2 
M203 70 


(reduced tc 50 for high elevation) 


Reduced recoil length at high firing elevations is necessary 
to prevent the recoiling parts from striking tl» ground. Thus 
the weapon supporting structure must be designed to withstand 
this “worst case" retarding force. In the instance of the M198 
towed howitzer, the worst case is the retarding force associated 
with the impulse from the M203 charge and a 50- inch (nominal) 
recoil length. 


In order to generate a constant retarding ferce, a variable 
throttling orifice is required. This is necessary due to the very 
high input force which must be attenuated. The procedure for 
determining a preliminary orifice "profile" is well established, 
(refs. 1,2,and 3) and is not repeated here. However, it must be 
stated that for any preliminary recoil brake orifice profile 
design, two parameters are essential: 


9 Total applied impulse (1Һ-вес) 
$ Total leugth available for recoil (in) 


This information is just the first of numerous assumptions ‘Tz 
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made in the design process. Ultimately due to practical 
realities of machining operations; hydraulic fluid 
compressibility, and turbulent fluid flow, many iterations are 
necessary to obtain a workable system- Computer models must be 
noditied with rluid discharge coetticients to obtain a suitable 
match to live fire data. 


As discussed previously. maximum force reduction is possible 
only when maximum recoil stroke is utilized. While the oritice 
profile is designed to do just that, many system variables fend 
to upset this ideal torce versus stroke relationship. These 
include 


° Variations in maximum impulse due to production 
tolerances in propellant manufacture, 


9 Propellant temperature variations due to varying 
climatic conditions which change the maxinum impulse profile. 2 
9 Manutacturing tolerances іп the tirottling 


oritice(s). 


° Temperature induced hydraulic fluid viscosity 
changes. 


In essence, the recoil brake throttling orifice is a pre- 
programmed device designed around ideal parameters “hich rarely 
exist. This may result in a uon-optimized force versus stroke 
relationship during firing. If weapon wcight and stability are 
not concerns, this approach to recoil design is satisfactory. | 
The lightweight artillery challenge however, requires an 
etiective and consistent energy dissipation system, 


4. Microprocessor Control of Weapon Recoil Energy, 

In order to ensure consistent recoil operat:on regardless of J 
the system variables mentioned above, the concept of closed-loop Ей 
feedback controi has been investigated., The concept of throttling 
recoil oil through a constantly variable orifice controlled by a 
Ғосаһаск systom is not new. It was described in a 1977 Rock 
Island Arsenal Technical Report (ref.4) but was never pursued. 

The reason for this, i. the author's opinion, was twofold. 


9 Microprocessor technology was not thor sutficiently ` 
developed 
° The high degree ot "tinesse" alluded to was 


unnecessary 


The advantages of electronic control include the tollowing: 
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a. Variable recoil length - Many artillery recoil 
mechanisms have provision to shorten the recoil stroke for high 
elevation firings. As discussed in section 3, this is necessary 
to prevent the recoiling parts from contacting the ground during 
recoil. The mechanism necessary to accomplish this reduction in 
stroke adds considerable mechanical complexity to the recoil 
system. This can be alleviated by providing the microprocessor 
with a weapon elevation ínput. The microprocessor could then 
perform the trigonometric calculations necessary to determiue 
available recoil length.  Recoil resistance could be adjusted 
accordingly. 


b. Counterrecoil control - The system used to control 
recoil could be used to control counterrecoil as well. 
Throttling of hydraulic oil can be programmed to ensure 
consistent buffing action and return-to-battery. This wouid 
eliminate the need for separate counterrecoil passages, thus 
further reducing mechanical complexity. 


c. Elimination of control orifices - Traditional 
artillery recoil mechanisms employ control rods or grooves to 
throttle oii. These orifices are precision machined and 
therefore costly to manufacture. Furthermore the final design 
is often done by tedious iteration with prototype hardware. This 
is both expensive and time consuming. Precision-made control 
oritices can be eliminated through incorporation of a constantly 
variable throttling orifice. 


d. “Tailoring” of the recoil force versus stroke 
profile- As discussed in section 2, one approach to weight 
savings to to apply a consistent load such that structural satery 
factors can be reduced. The weight savings gained through 
incorporation of reduced safety factors is of little value, 
however, ií weapon stabiiity is not maintained. Practically 
speaking, the weapon must not jump" or "hop" when fired. 
Unfortunately, a lightweight artillery piece is, by nature, more 
prone to instability due to the reduction of overall mass. 
Consider the free body diagram of forces acting on the weapon 
(fig.2.) 
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By summing moments about point А it сап be seen that the gun 
system will remain rotationally stable as long as the counter- 
c ockwise moments from the weights of the recoiling parts (Wr), 
and carriage assembly (We) exceed the clockwise moment from the 
force that the ¢cecoiling mass exerts on the carriage. The 
condition Гог siabilii» is 


RECOTLFORCELJ<W(L-*)+W.L, = © А 


Note that as recoil progresses to the right, the counter- 

clockwise moment trom the r.-coiling weight decreases due to the n 
moving center oi gravity of the recoiting parts; therefore the 
condition lor stability is most critical at the end of the recoil 
stroke. While this condition does not override the primary goal 
of maintaining a consistent recoil force as low as possible, it 
demonstrates that the recoil force should not increase during the 
final portion of the stroke. It would be entirely poss.ble for 
the microprocessor to tailor the force versus stroke profile to 
improve stability of a lightweight weapon. 


The constantly variable fluid throttling orifice thus 
described is envisioned as a microprocessor-activatod servovalve 
operating as a fluid bypass. A schematic of this system is shown 
in figure 3. 
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7 — RECOIL BRAKE 


FIGURE 3. - MICROPROCESSOR-CONTROLLED BYPASS VALVE 


Assuming various factors such as recoil force, position, 
velocity and acceleration arc available to the microprocessor in 
rea time during the recoil stroke, optimum recoil energy 
dissipation can be programmed. Two control algorithms have been 
developed (ref 5) to accomplish this and aru described below 


Level ! Control - Maintain a Preset Recoil Force. The recoil 
force or the force transmitted to the weapon supporting ез 
structure сап be monitored and maintained at а preselected 

value. This value could be mathematically or empirically 

determined such that the total available recoil stroke is always 

utilized for a specific impulse input. If actual recoil force is 

less than the preselected value, the servovalve would be 

commanded to shitt to the closed position If recoil (orco is 

greater than the preselected value the valve would shiit to the 

full open position. 


Th controi algorithm described above is very simple, vet it 
has the potential for true optimization of each recoil stroke 


through the concept of closed loop teedback control. 


evel 2 Control- Compute Recoil Force During Firing A 1 ore 
е 1 approacn to effective recoil energy dissipation is 
to dynamically determine required recoil force during recoil , 

This has the advantage of not requiring advance knowledge of what 


impulse is to be expected. 


This technique equates the mechanical energy of the recoiling 
nass to the amount of werk necessary to stop the recoiling 
parts. The work energy reiationship is tirst established. For a 
simple, one-dimensional recoil mechanism, motion is governed by 
Newton's second law 








^i 
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gud. * . .. 
Ta Ё Ж. М X (2) 


where: F= sum of all forces acting on the recoil 
mechanism in the direction of motion 
(1b) 

Mass ot recoiling parts (slugs) 
Acceleration of recoiling parts 
(it/sec ). 


ж 
ТІМ! 


Figure 4 illustrates a treo body diagram of torces acting on a 
typical artillery recoil mechanism. 


C 
Sx ce DIRECTION OF 


p Figure 4. - Forces Acting on au Artillery Recoil Mechanism 





where: B(t)=Propellant gas force as a function 

ot time (1b) 

K =Recoil resistance (1b) 

Е -Mechanicai friction due to recoil 
mechanism movement (lb ) 

W Weight ot r coiling parts (lb ) 

Ө =Angle ot elevation of weapon 
(degrees) 


Arbitrarily assuming a positive direction as shown in figure 3, 
equation 2 can be re-written as 


M.x = btt )« Ve / sing- -K- E (3) 


Recoil неа m t with the พ Д3 substitution: 


oN 
== i 
S иаи — Re)twsin@ -К- © 
M gv. = В) А | 
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44.» 
Integration үа respect to distance yields "s 


БЕСІ 


Ву selec "m. m. initial T LUE อ อ ค อ อ อ อน at some 
arbitrary distance,x, (where the velocity 15 V), a tinal 
condition at maximum recoil stroke.xnax (where recoil os ae is 
Zero), and assuming recoil torce (k) is to be maintatned at 
constant value, equation 6 becomes 


mex 
- I^ ES - aide W sine rax) = Крета х)- fd 


Mn this Buation tor recoil resistance (k) vis Ts in 


/ та xX Aa Х 


= س‎ Е УИ жах жү у?Ө (8) 
XMmax~X, 
Equation 8 becomes the ee equation" used to compute 
recoil resistance during the recoil process. While the 
microprocessor SE probably update the computation quickly 
enough, the mechanical transducers. A-D converters and servovalve 
action may be "too slow" to effect the proper change. Although 
this has not yet been established through сезі, ii ceriainly i 
would be desirable to simplify equation & to reduce computation x 
time, This could be accomplished using the followimg logic- 


A comparison of the magnitude of forces used in equation 8 as 
a function of recoil streke for a typical artillery recoil 
mechanism is depicted in figure 5, It is evident that the 
friction and weight terms are not large contributers to the work- 
energy relationship. Realizing that friction tends to reduce the 
recotl resistance force while the weight contribution tends to 
increase it, tt would be prudent to include th. weight 
computation and to simplify equation 8 „omit friction. This wil! 
ensure the stability requirement is met since recoil force will 
not inerease as recoil action progresses. 


Further simplitication of equation 8 is possible if the 
breech force component is ignored. This may appear radical since 
the breech force is tho largest force applied to the weapon; 
however, it acts for only a traction of the total recoil 
stroke.Furthermore, the computed value would only be correct 
atter Blt) b. comes negligible. Ignoring breech force would 
require the servovalve to control сессії force to a preset value 
until the majority of it (the breech force) has been applied. 
This value could be the maximum recotl force the weapon is 
designed to handle. Incorporating these decisions; equation 8 can 
be rewritten as 
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The obvious question is at what point during recoil should 
control transition from preset rodpull to the calculated rodpull 
computed in equation 9. This could occur when the acceleration 
ot the v. coiling parts changes sign trom positive to negative. 

li the weapou employs a muzzle brake, transition would occur when 
the muzzle brake activates. lt a muzzle brake is not used, 
acceleration will change sign when the breech torce is reduced to 
à level below the combined recoil resistive force and trictional 
lorce 


There are both physical and computational problems with 
equation 9 during the final portion of the recoil process which 
need to be addressed. Due to the inevitable leakage of piston 
Seals and variations in frictional and recuperator forces the 
recoiling mass would never stop at precisely the specified xmax. 
This, in turn, causes К to become undefined as xmax is 
approached. ln order to prevent this,a controlled closure ot the 
servovalve should be provided at a specified distance from xmax. 


This approach to recoil mechanism control has been simulated 
(ret£.5). The results асе eucouraging, and prove the viability of 
this approach to recoil control. 


Implementation of the microprocessor recoil approach results 
in some concerns. They include: 


a. The substitution of electronic complexity for 
mechanical complexity. It must be acknowledged however, that 
electronics have proven their dependability in the acrospace and 
automotive fields. therefore,This may not be a deterrent. 


b. Input ot energy is necessary to "manage" recoil 
energy dissipation.( Electrical energy is required to power the 
nicroprocessor and a high pressure fluid source is required to 
shift the servovalve.) 


с. Micropreecessor tailure could be criteal it а 
redundant backup is not provided, 


These comments are not meant to discourage the approach, but 


merely to point out unique problems not previously encountered In 
recoil mechanism design. 


Lowered recoil force is maudatory tor a lightweight artillery 
piece. There is no magic involved in achieving this; lowered 
force levels can be achieved only by increasing eftective recoil 
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is always utilized. 

Two approaches to lightweight recoi! mechanisms are desertbed 
in the tollowing paragraphs. These concepts are not new but with 
the addition ot microprocessor control, become viabie 
alternatives for a lightwelght artillery application. 


Super Long Recoil - Using a conventional cycle; i.e. 1ісе-гесоіі- 
return to battery, but incorporating 120 inches for recoil 
travel, the torce imparted to the weapon can be reduced by 50 
percent.This reduction iu force is due to more than doubling the 
Length ot recoil. An artists sketch of one possible approach is 


shown in figure 6. 


Utilization ot this approach does have drawbacks. these 
inelude: 


а, Drastic changes when compared to а conventional 
weapon carriage. For example, a significant increase in veapon 
height. This would, in turn, negatively а іесс the vehicie's 
silhouette, making it more "spottable". 


b. Loading the weapon would be more diiticult because 
the breech is 10 ft. off the ground. This nrohiem could be 
aileviated by programming the microprocessor to bring the 
recoiling parts to a load position,toiílowed by re-positioning for 
firing. 


5010 Recoil- Ferhaps the most intriyuiay approach to 
lightweight artillery recoil systems is the concept of soft 
recoil. The basic idea ot soft recoil is embodied by its 
fundamentally different sequence of operations when compared te 
the conventional recoil cycle. This dittercnee is graphically 
illustrated in figure 7. 


In conventional recoil mechanisms, the firing momentum is 
directly transferred to recoiling parts that are at rest prior to 
tiring. This momentum is then absorbed by the supporting 
structure, and the recoiling parts are brought to rest. In soit 
recoil, approxiamately halt the moneatum of the projectile is 
applied te the recoiling parts by accelerating them іп the firing 
direction prior to firing the projectile. Firing momentum thus 


transferred to the toward moving recoiling parts stops and 
propels chem to the rear. In ettect, recoil length is doubled 
and imparted foward momentum is absorbed. The soft recoil cycle 
can reduce forces applied to the supporting structure by more 
than hall that imposed by a conventional system.(ref. 6 ) 
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Length- li this can be accomplíssed, the microploceszor can sup- 
ervise the recoil process to ensure that the total recoil stroke 
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Due to cail-sate devices incorporated in past soft recoil 
mechanisms,the existing mechanisms tend to be heavy and complex 
simply due to the practical realities of building such a 
weapon. Specifically, two conventional "recoil mechanisms" are 
required, one at each end of travel. This is due to the potential 
for two overload conditions which may occur as follows: 


a. The weapon must be programmed to tire at a foward 
recoiling parts velocity commensurate with the propelling charge 
energy teleased during combustion. The possibility exists to 
fire the highest zone (maximum impulse) at lowest foward 
velocity. This results in a rear overload condition. This 
condition can also occur due to propellant "cook off", in a hot 
tube prior to the initiation of foward velocity. in either case 
the imparted energy must be dissipated usually through а 
traditional oil throttling process. 


b, After foward velocity is initiated, ignition delay 
or misfire (failure to fire) will cause a forward overload 
condition. This energy must also be dissipated again, usually 
through an oil throttling process. 

In essence, the soft recoil approach requires a high degree 
of controi during the recoil process. Timing is critical. In 


addition, stability problems are inherent in the design if it a 
doves not function properly. K 2 


It is suggested that microprocessor technology can be 
incorporated to "manage" the entire soft recoil cycle. This 
management could consist of the following: 


a, Prior to firing,ensure that the proper zone, i.e. 
requir:d foward velocity, is input into the weapon control. This 
could consist of a keyboard excerise, or manual verification of a 
firing command. This could eliminate the rear overload 
condition. 


b. Electronically monitor foward velocity and initiate 
firing. The microprocessor could monitor propellant temperature, 
changing frictional torces and store past firing performance so 
that firing is initiated at just the right instant based on the 
action of the particular weapon. This will ensure optimum cnergy 
dissipation. 


c, Assuming a misiire or ignition dalay occurs, a 
controlled shutdown is initiated by throttling oil through a 
servovalve. The recoil shutdown torce profile could be designed 
such that weapon stability will be maintained. This could Бе 
accomplished with “tipping” sensors to maintain overall 
stability. 


ee 
^ 
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d. On~board diagnostics will ensure the weapon wil! not 
be destroyed due to failure of a replacable fail-soft 
component. Assuminmg an overload conditíon occurred, such that a 
sacrificial energy absorber was used, for instance,in Lieu of, or 
in addition to, servovalve throttling, firing would be prohibited 
until a replacement was made. 


These "management" options are only suggestions of what is 
possible with microprocessors controlling soft recoil. It is 
essential to review state-of-the-art soft recoil to fully 
appreciate the problems associated with this concept. The U.S. 
Army has constructed па type classified a i05mm soft recoil 
howitzer, designated M204. An excellent synopsis concerning the 
development of this weapon and can be found in referenco 6. 
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Abstract 


In order to arrive at a decision in regard to the development of a flick 
гаттег for the Howitzer Improvement Program the compatibility of the flick 
ramming environment with large caliber fuzes was investigated. The fuzes 
investigated were the M557, М759, M577, ХМ762, M732 and M572. The only 
possible area of concern uncovered in this investigation was the М1 delay 
plunger (used in the M557, M572 and M739 fuzes) and Its successor the impact 
delay module (IDM) which is used in the M739A1 fuze. 


The centrifugal pin in the М1 delay plunger has been undercut to prevent 
functioning in the absence of a spin environment. A range of setforward 
decelerations, plunger spring prefoads and frictional forces were simulated. 
In al! cases the undercut сп the centrifugal pins caused the plunger assembly 
to lock up before it could function, 
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Robert X. Brennan 
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Picatinny Arsenal 
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1. Introduction 


To insure the interchangeabiiity of current and future US & NATO artillery 
fuzes ammunition used in the flick ramming/autoloading environment, PM-CAWS at 
Picatinny Arsenal initiated en experimental and theoretical investigation of 
the compatibility of the M557, M739, M577, ХМ762, M732 and M572 fuzes with the 
flick ramming environment, 


In May 1983 the deceleration G-level required to seat the 155mm round in 
the M199 gun tube of the M198 Towed Howitzer equipped with a flick rammer was 
measured with triaxial! accelerometers. Figure 1 shows the accelerometers 
mounted in the fuze well of the 155mm M107 projectile. 


Ten ramming tests were conducted through a range of elevation angles 0 to 


70 degrees In steps of 10 degrees. These tests ylelded a total! of 80 


deceleration-time traces. The traces were digitized by the Hybrid Division of 
ARDC. A subset was used to provide a quantitative time dependent forcing 
function for the computer simulation. The highest magnitude of pulse observed 
on these traces was 1800 G's, 


As part of the study the fuze mechanisms were discussed with the 
respective project engineers. Those components about which there was some 
doubt as to whether they could withstand the flick ramming environment wert 
subjected to analysis and/or testing. 
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The first step in determining the response of fuze mechanisms was to 
measure the deceleration u's required to seat the 155mm round In the M199 gun 
tube of the M198 towed Howitzer equipped with a "Flick Rammer". In the fuze 
weil of the 155mm round the deceleration G's were monitored by a triaxial 
accelerometer. The mounting of the accelerometer Is shown In figure 1. Ten 
flick ramming tests were conducted at quadrant elevation angles of 0°, 10°, 
209, 309, 409, 509, 60°, and 709. !17 These tests yielded a total of 80 
deceleration-time traces. The results of these tests are shown in figures 2 
and 3, Figure 2 shows the maximum measured values of the deceleration at 
the fuze well as a function of the quadrant elevation angie of the gun tube. 
Figure 5 gives the range of ramming speeds generated by the flick rammer. E 


These 80 anaiog traces were digitized and processed by the Hybrid 
Division of ARDC using a program called "Hydra™*’. The maximum axial БЕ 
deceleration amplitude was about 1800 G's. А typical trace for quadrant 5 
elevation angle 30° is shown in figure 4. The x and y axis are transverse to 
the round and they represent balloting forces The z component is directed 
along the longitudinal axis of the round. The HYDRA program can expand the 
time scale and the plot is shown with the expanded time scale. Since the 
round deceierates when being seated by the flick rammer the sign asscciated EM 
with the Gis should be negative. Accordingly, in the analysis only the AE 
negative portion of the z component was used as input to the analysis. The 
positive portion of the trace was assumed to be noise. These pulses have a 
typical time duration of Z milliseconds. These durations agree well with 
the measurements conduct gd by GE for the Bundesamt fur Wehrtechnik and 
Beschaffung (BWB) 1981 which resulted in pulse durations of 1.8 
milliseconds. In а series of tests conducted at Мерреп in 1983 (4) by the 
BWB, 40 rounds with the M739 fuze set in the delay mode were tlick rammed at 
33 feet/second and all the rounds subsequently functioned properly cn 
impact. Twenty fuzes were тоџпіє і on the 115А1 projectile and fired at high 
and low temperatures. Another 20 were fired with tne M549A1 RAP round at 
high and low temperatures. The ramming Gts were not measured, but they 
could have been as high as 4,000 C's. 


As a comparison, measurements "95 also made in the M185 gun tube 


equipped with a hydraulic ғаттағ( . This trace was also 
processed by the HYDRA program. The pulse duration is again 2 
milliseconds, but its peak is only 218 G's. 


Euzes Investigated 


The fu: ss which we cheched for compatibility with the fiick rammer 
were the M557), 73906), echt, хМ762 9, and M732 игез. 

Discussions were held with the project engineers on these fuzes, and those 
Sub-assemblies which could possibly fail as а resul! of the flick rammer set- 
forward ferce were identified. The following sections of this report 
summarize the analysis and laboratory testing of these parts. 


M357 PD Element 


Б 


Figure 5 shows a sketch of the firing pin In the M557 PD element. Under 
set forward deceleration it is possible that the firing pin could be forced 
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through the washer and the crimp on the nase. Using a spring tester the load 
on the firing pin head was increased until the crimp sheared. The shear ing 
force measured was 69 Ibs and since the firing pin weight is 0.021 oz it would 
require 52000 G's of setforward deceleration to drive the fi-ing pin through 
the crimp. Since the highest flick ramming G's never exceeded 1855 G's, this 
failure mode can be disregarded. 


M739 Fuze Eiring Pin 


A possible mode of failure in the M/39 fuze Is for the firing pin head to 
come off the firing pin. According to the drawing specifications for part 
#9294606, it would take 80 Ibs to pull the head off the pin. In the actual 
asse ibly as shown in figure 6 more force would be required since the firing 
pin head is а!зо held by the firing pin tube. Since the pin weighs 0.03 oz, 
using 80 105 as a conservative estimate it would require 42000 G's setforward 
deceleration for this mode of failure to occur. Therefore, this type of 
failure can also be ruled out. 


M739 Fuze, Firing Pin Support 


Another potential mode of failure in the M739 fuze is tor the firing pin 
suppoct to collapse. According to the drawing part number 9258614, it takes 
about 160 Ibs to collapse the firing pin holder. Since the detonator housing 
weighs 0.075 cunces, this mode of failure would require a set forward 
deceleration of 21,000 G's. 


M577 Euze Time Assembly 


tube and the timer hairspring. A skeich of this mechanism is shown in figure 
7. According to drawing part number 9236712, the timer hairspring tube must 
be able to withstand a load of 800 ibs. The weight of the tube was calculated 
to be 0.054 ounces. To break this press fit would require a setfor ward 
deceleration of 235000 Gis. According to the same drawing, it would require 9 
165 to pull the + mer hairspring free of its collar. The weight of the 
hairspring was calculated to be 0.002 ounces. It would, therefore, require 
65000 G's to pul! the hairspring free. 


4732 Fuze 


Three M732 inert fuzes were mounted cn inert 155mm projectiles and rammed. 
All three {игез locked undamaged, hadn't turned and felt tight in a hand grip. 
The three fuzes were disasembled ang nothing unusualy was observed, Recent 
air gun tests at Picatinny Arsenal indicate that the time setting torque 
could be ร อ ท เอ น ร เห degraded when the fuzes were subjected to setforward 


decelerations between 1000 G's and 2000 G's for a duration of 5 milliseconds. 
ХМ762 Fuze 
This fuze is in the engineering development phase. Discussions with the 


project engineer did not point jo any part of the mechanism as being possibly 
sub ject to failure. 
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М1 Delay Plunger 


The M1 delay plunger which is built into the inventory of M557, M572 and 
M739 fuzes provides a functioning time delay. A sketch of this device is 
shown in figure 8. Since the delay element is designed to function due to the 
setforward deceleration generated when the projectile impacts and penetrates a 
wall, It is possible that if could also function when the projectile is seated 
by flick ramming. The major difference is that the projectile has a high spin 
when it impacts a target whereas there is no spin when it is flick rammed. As 
а result the setforward deceleration would have to be sufficient to retract 
the centrifugal pins. Since the centrifugal pins are skewed at an angle of 15 
degrees in the direction of flight it is possible that there would be a 
component of the deceleration G's sufficiently large to cause the pins to 
retract. As the centrifugal pins retract, the spin lock detent (not shown in 
the figure) will play no role since the plunger immediately rides forward over 
the plunger support. Once the plunger goes forward the delay element would be 
criven into the firing pin. The initiation of the delay element would in turn 
damage the S&A causing the round Tc be a "dud". Qualitatively there is the 
possibility of a malfunction of the M1 Delay Plunger and the question remains 
to be answered if quantitatively there is enough energy for this to occur. 


Centrifugal Plunger Pin and Spring 
In order to proceed quantitatively the differential equations for the 


motion of the centrifugal pins as well as for the plunger are required. The 
differential equation for the centrifual pins Is given by: 


1) мХ = -АХ - mg Gep – к mg Са (6) ces, — ma б. (Е) eoo Ө, + m4 Ga (é).oxe- Ө, 


x(0) = ๐ 
x(0) = O 


m is fhe mass of the contrifugal piunger pin, 0.116 X 10-4 Ib 


г 


sec2/inch 
x is the displacement of the centrigual plunger pin, inches 
k is the contrifugal plunger pin spring rate, 0.19 Ib/inch 
g is tne gravitational acceleration, 386.4 inches/sec 
бер centrifugal plunger pin spring preload, 17 G's (M!L-P-10480) 
6,(%) is the longitudinal component of setforward deceleration 


Ax coefficient of friction 


9, angle between the centrifugal plunger pin and the perpendicular 
tc tho spin axis, !5 degrees (see figure 8) . 


С-(1) transverse or balloting component of deceleration 
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For the purposes of analysis equaiion 1 can be written in the form 


2) Kaw X — 9 Gip- мк G, (Ó) ce» @, + д ( Ga (€. Ө, - G. (4) сео ө.) 


where МӘ 
3) Wee JE- 28.4 hera 
Plunger Assembly 


The differential equation for the motion of the plunger Is 
4) mp y = -kp y = трд G; + тра Ge (t)— My mp 4 Gr fé) 

Шу mass of the plunger body 2.74 X 107^ 1b-sec/ Inch 

РР plunger restraining rate 3.9 Ib/inch 

С, plunger restraining spring preload 13 G's 
For the purpose of integration equation 4 can be written as 
5) ў =- ову- و‎ Gs + 4 Ga (©) му و‎ Gr (€) 


where Wp = 20 herz 


Results of the Computer Simulation 


(г; 


The differential equations 1 and 4 were solved simultaneously employing а 
special purgqse computer program for solving ordinary differential equations 
cal led ACSL The results for test 29 rammed at QE of 30° are shown in 
figure 9, In figure 9 three curves are superimposed. These are the axial 
setforward G's (650 G's max.) and the displacement of the cen?rifugal plunger 
pins together with the plunger displacement. To cause initiation of the delay 
element the centrifugal pins must displace 0.095 inches before the pin and 
plunger assembly are rammed against ihe plunger support. іп figure 9 the 
coefficients of friction of both the centrifugal pins, My , and the plunger 
body, My , were taken equal to 0.1. Тһе relaticnship of the plunger assembly 
with respect to the plunger support is shown in figure 10. The centrifugal 
pins must release before the plunger moves 0.059 inches. In figure 9 the pins 
displace only 0.044 inches before being rammed against the plunger support. 


If the G trace were to be multiplied by 1.5 it would peak at 1000 C's. 
Maintaining the coefficients of friction wx and Му equal to 0.1 the results of 
а simulation are shown in figure 11. Once again the pins are rammed against the 
plunger support before they free the plunger. 


The simulation was rerun using the 650 G pulse and the coefficient of 
plunger fricticn was set equal to 0.5. Again the МТ delay plunger does not 
function. The simulation was repeated with the same coefficients of friction 
and а 1000 G pulse. Although the pins displace somewhat further there is 
no functioning of the M1 delay. 
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However, if the preload оп the plunger spring 15 increased from 15 G's to 
100 G's the МТ delay will almost function. The main effect of the high 
preload is to delay motion of the plunger until the pins release. 


A higher preload on the plunger spring (150 65) allows the pins to clear 
the plunger support permitting the plunger to move forward, However, as shown 
in figure 12 the plunger does not close the 0.2 inch gap to the firing pin as 
the plunger spring is too stiff. 


If the conditions of the previous simulation are repeated except that the 
peak G's are increased to 2000 G's the М1 plunger assembly will function, 
This is shown in figure 15, Тһе peak kinetic energy is over 200 inch ounces 
which would initiate the М1 delay element. On the other hand a 2000 9 pulse 
and a 13 G preload on the plunger causes the pins to jam the plunger. 


In conciusion, the race between the plunger and the centrigual pins is won 
by the plunger unless the plunger spring stiffness is increased by a factor of 
10. 


Effect of Plunger Support Interference 


The question remains to be answered as to whether the М1 delay assembly can 
function when the retraction of the centrifugal pins is interfered with by the 
plunger support. The head of the centrifugal pin is undercut by 13° which 
produces an additional normal force on the pins, Taking this additional normal 
tigure 14. For convenience the angle of the pin and the angle of the undercut 
were both taken as 6, = 159. From this figure the equation of equ!librium can 
be read directly. 


6a) М ceo 8, - ma б» Cov 86 - mq Gr am 0, = о 


6b) mx’ = mg Gz 4in8, - Rx - тд бір -mg Gr Сюд, = их New p,- Nain Ө, 


In these two equations terms of the order of Ay aim O have been 
neglected. Combining these equations yieids: 


7) X” = es X “= 9 Gip- 3 Gr ase D- инд (Gp cre O04 Grom ө.) 


Since ali the terms in this equation are negative the pins cannot move. 
Only the bal Ictting force G (+) can change sign. Since the two pins are 
diametrically opposed, if the ballotting drives one pin out, it drives the other 
pin into engagement with the plunger support. 
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Figure 3 








чв amm aw ља w s. 


Дл R.J K ET A ПУ SL лас TH рыга oe es es ar V T แร า พ ті; 


= V ALT 


"AT. 


25311480 60/51/70 


050990: 7590 5990 S 








т-а ma pet V QR сла CR + 





664° 49- 436104 XYEN 
22854256 135 Nd 20015 ХУЗа 


ааа‏ م س 








oo BEATE piepers ฟา ๆ ว 


Suave) Go 0 


ПОО. NOIIUSJIJOGU ข่ ว ม น บ ปี 


PN NN ا‎ 
"IDN Z-1H 
62 "DN ASIL 


==: e ee = 
IO 00б носон JNO Yung 


er 





U XE. 


wu s 


CSGNO2:3) SWI 
659.0 


5.9.0 


. 23-0 








V-68 





- M557 FUZE PD ELEMENT 


DETONATOR, M24 











ООР | == 
BALLS fF АА аран 
@ | FIRING PIN 


M55 PRIMER амы 





Figure 5 











M739 FUZE FIRING PIN 


FLASH TUBE 


M99 STAB CETONATOR 


RAIN CAP 
=Z Ж 





CROSS BARS 


FIRING PIN 5 


—Т 


Ус! FIRING PIN 


9% 
| 2536 
hee .CETONATCR 


[FIRING PIN ASSEMBLY 
: | 92525.8 406% 





A 
| 1 Му | 
! ER 2 ด 
LOA esas ละ ค ซะ ร 
Е у E eg „Сте 6 
} их ST 
1 EMEN EE & 24 
น a+ 
X 
Стыр BEN 282548504 — z ' ถ ว ก จ ชู «๕ 
NOTES уфа » Mor De P qs THEA (SING HIN 
--!036 ๑ 032----= ฯ “ое E 
Figure 6 
DA 
MM y 
V-70 








SETTING KEY 





/ 


ENT 


CRUSH ELEM 





EY 


SETTING Ki 
FLANGE 


` SQ ELEMENT 


TIMER HAIRSPRING TukE 


TIMER HAIRS PRING 





-71 


u 
` 











M1 DELAY PLUNGER 1557 AND M739 FUZES 
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Figure 8 
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Figure 9 
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The importance of computer simulations of performance of 
complex mechanical systems has long been recognized, The Army Materiel 
Systems Analysis Activity (AWšAA) has deveioped an armored vehicle 
mode] for studies or automotive performante, hit probability estimations, 
and relatea parameters. Since дип aynamics is known to play an important 
role in weapon accuracy, it was deemed apzronriate to replace the rigid 
qun tube portion of the AMSAA model with a relatively simple qun 
dynamics model supplied by the Ballistic Research Laboratory (BRL). 


In this paper we describe the intejrated model, including the 
lim: tations and simplifying assumptions concerning the qun dynamics 
portion, and vresent comparisons between computer simulation output of 
muz7le motion with actual rezerded muzzle motion from a moving vehicle 
Ove, various types of terrain. 


2. QVERVIEN OF THE INTEGRATED WEAPON ARMORED VEHICLE MODEL (IMAVM) 


The qun dynamics medel used in the simulation has as its basis 
a modeling technique deveicped by Boresi (1). This technique involves 
the uss 5* а finite element method to simulate the dynamic response 
of the qun tube to nounat motion ana baltistic loading. For the purposes 
of the си: rent computer modei, numerous significant changes have been 
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made by one of the authors to the original modeling technique. In 
particular, the recoil, the interior ballistics, and the mount motion 
sections of the model were completely reworked to reflect mcre realistic 
conditions, 


The original interior ballistics portion of the model used the 
so-calied standard equations for pressure and projectile acceleration 
and travel from Corner (2). These equations suffer from two main diffi- 
culties. Namely, the phasing between pressure and acceleration is not 
accurate, and an additional inflection point which results in incorrect 
curvature near the critical muzzle exit point. To correct this problem, 
the simulation was modified to use actual pressure, acceleration, velocity and 
travel curves for particular projectiles. 


Figures 1 through 4 depict the modified and augment Corner 
interior ballistic responses for the M392 projectile used in the simulation. 
These data curves were obtained hy averaging, using ammunition lot acceptance 
data. The fact that the interior ballistics were modifiea necessitated 
reworking the recoil section of the code, since the interior ballistic 
forces drive the recoil mechanism. In this case, the same scaling between 
projectile acceleration and recoil acceleratiun was used as in the original 
version. The breech is modeled es а solid mass with axis offset from the bore 
axis. 


To accommodate vehicle motion input to the trunnions, the model was 
modified to accept time series data, instead of using a Fourier expansion 
as forced base motion. In this way, actual data recorded from vehicle 
road and cross-country tests could be input to the gun dynamics portion of 
the model for comparison with response to simulated vehicle motion. This 
also allows actual firing data from moving vehicles to be compared with the 
simulation. It should be noted that in the current model, only gun tube 
motion in the vertical plane is simulated. It is anticipated that 
future developments will include a full six-decree-of-freedou qua dynamics 
model, for more rea'istic simulation. 


The gun dynamics model has been run independently of the vehicle 
model for verification of its output in terms of response to ballistic 
leading. On the basis of available experimental data for comparison, the 
моче] gives an adequate representation of the actual firing dynamics of tlie 
M68 105mm tank cannon. Comparisons of test data with computer runs of 
the full simulation snow gocd correlation, especially of muzzle vibration 
frequencies. 


The motivation for inteqrating the Boresi gur dynamics model 
into a detailed engineering model/simulation uf an urmored vehicle stems 
from an analytical need to: 
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a. Characterize/Quantify this error source and it's impact 
on weapon system delivery accuracy, especially for dynamic situations 
such as firing on the move, and 


b. Study possible fire control modifications/al gorithms 
designed to correct for muzzle flexure induced errors resulting from 


base motion disturbances and interior ballistic efforts. Besides enhanc- 


ing the delivery accuracy capability of a weapon study, this study has 
the potential of increasing taraet servicing rates by allowing firing 
of the weapon under extreme base motion disturbances resulting from 
high terrain severity conditions. Figure 5 depicts the functional 
block diagram of the "Integrated Weapon Armored Vehicle Model" (IWAVM) 


that has been jointly developed by the US Army Materiel Systems Analysis 


Activity and US Army Ballistic Research Laboratory. 


The digital computer simulation of an armqred vehicle weapon 
system that is being used in the study contains representations of 


the vehicle and target motion, suspension characteristics, weapon/turret 


servo drives, fire control, sight/reticle servo drives, human qunner 
model, human driver performance characteristics and terrain profiles. 
Briefly, the basic subroutine modules are: 

T a inertias, cen 


ST, in which total 5 
га run are initia 


a. SUBROU IR и 
nbal nd scenario conditions 


NE F 
gravity, mass unb се a 


+. 
ler 
liz 


b. SUBROUTINE MOTION, which generates the motion of the road 
wheel over either an APS Bump Course or Micro Terrain Profile and the 


anquiar motion of the hull on the suspension. Іп turin it calls SUBROUTINE 


BUMP, in which the terrain or APG Stabilization Bump Course is modeled, 


and SUBROUTINE SUSPEN, in which tne torsion bar, volute springs, suspension 


stops, and dampers are described. 


c. SUBROUTINE DRIVE, which describes both the elevetion and 


of 
ed. 


azimuth stabilization drives to include the analog compensation networks, 


hydraulic servo valves, load dynamics, gear box characteristics and sensors. 


d. SUBROUTINE GUNNER, in which the human dynamical responses, 
decision making, and action and visual thresholds are modeled. 


e. SUBROUTINE SIGHT, which simulates the gun-director type 


fire controi in the elevation sight axis and the driven-reticle model in 


the azimuth sight axis. 


f. SUBROUTINE RETCON, which simulates the digital control laws 


programmed in the ballistic computer for commanding the azimuth reticle 
servo. 


q. SUBROUTINE OFFSET, which emulates the Ml elevation and azimuth 


algorithms for computing the ballistic offsets for the kinetic-energy 
and heat rounds. 
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h. SUBROUTINE TRAJEC, which computes required weapon to target 
offsets and takes into account target motion as well as additional 
velocity components imparted to the projectile. 


i. SUBROUTINE GUNDYN, which simulates the independent weapon 
tube flexure responses based on the Boresi model formulations. The 
model has been modified to accept vertical and gun tube translational 
motion generated at the trunnions by the armored vehicle simulation. 


3. SUMMARY 


To date the Boresi gun dynamics model has been successfully 
integrated with the Ml engineering armored vehicle model and exercised 
for various verificaiton/validation scenarios. Figures 6 and 7 depict 
actual and simulated power spectral analysis responses for the M68 
105mm weapon. The scenario conditions were gravel road (0.5 inch RMS - 
Waterways Experimental Station Terrain Severity Criterion) at 15 mph. 
Both simulated and actual responses generally exhibit the same resonance 
frequency range and amplitude levels. 


Further verification and validation work as weil as additional 
refinements to the mode? are being pursued prior to application of the 
IWAVM simulation. Future work will focus on extending the model to 
include the extended 105mm as well as the 120mm gun tubes. 
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|. INTRODUCTION 
Since the seventeenth century many people have noticed that 
gum projectiles do not always land on their expected taraet. 
Attempts to improve qun accuracy are concerned with the 
scurces and mechanisms which cause a given projectile impact 
distribution at time target. 
% The flight of the projectile tc the target can be considered 
in three phases: 
(a) Interior ballistics - projectile motion within the 
barrel. 
(b) Intermediate ballistics - projectile launch, propellant 
gas blast and, for sub-calibre rounds, sabot separation. 
(a) Extericr ballistics - dewnrange flight and target 
impact. 


In this paper we will consider the motion of the gun barrel 
and its mounting up to projectile launch. The computer program 
ГАМА 34 was used. This program was written and deveioped by 
The Royal Military College of Science and a paper describing 
an earlier version was presented to the 2rd US Army Symposium 
лп Gun Dynamics (1). 


~. CESCRIFTION OF COMPUTER MODEL 
The mathematical model represents a flexible gun barre: 
mounted in a rigid body cradle. The cradie is able ts rotate 
about its trunnions, its rotational movement being rear.ted t, 
an elevating gear stiffness coupled to a slipping clutch «see 


Figure 1). The barrel can have breech and muzzle masses 
AT, 
Mh attacmed to it. 
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At the barrel/cradie interface are two identical cylindrical 
elastic bearings of finite length. The hearing positions ard 
parameters can be varied. 







CARREG. 


อ อ ล 64 ๐ 3 LTRUNNIONS. 


6 เล ห ศ ร ์ เอ น ่ 
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Figure 1 Schematic of the Breech, Barrel and Cradle 





The theoretical basis cf the model is the Euler-Bernoulli 


theory of thin Seams. The model is subiect to both flexural 
and axial vibrations. The resulting differential equations 
with their boundary conditions are numerically solved using an De 
implicit finite difference algorithm. The program is now very ๐ 


user friendly and data files are easily compiled, edited and 
recorded. 


The algorithm used for the solution of the differential 
equations is efficient and robust. The method cf solutisn was 
described in the earlier paper. 


Z.FARAMETERS INVESTIGATED 
In а pre. Lous paper presented to the ТТСЕ ETAS Second Worl shop 
โก September 1984 CU + the predicted effect of local 
stiffening and mass redistribution of the gun barrel was 
examined. Jne parameter not examined in the f TAG paper was 
that of cradle bearing stiffness. 


t should be emphasized that the gun system parameters used in 

116 investigation are for illustrative purposes anly and dc 
поё represent a current tank gun nor any proposed future tank 
jun. Тһе effects аце to changes in three parameters will be 
examined: 


(a) Gun barrel stiffness 
(b? Cradle bearing spacing = 
(=) Cradle bearing stiffness Ф 
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4.СНОТСЕ OF QUTFUT VAFIABLES 


m ` کے‎ . . ——— 


If we intend to produce a ‘better’ gun we must decide what 2 
measured criteria we will accept aS producing an improvement. ы 
In this paper the vibration and yaw angle of the projectile К 
within the. gun barrel will not be considered. although some Y 
work in this area does look promising. The assessment will 


simply compare gun muzzle transverse velocity, muzzle slope 
and gun cradle rotation. 


It could be claimed that a gun system giving а consistenkiy 


high ‘gun jump’ figure (see Figure 2), would give a constant 
bias which could be compensated for within the signting 
system. It should be mentioned, however, that s qun which i 


sensitive to small production Variations will probably lead to 
variable accuracy whereas а gun which demonstrates a lower 
sensitivety to changes in these parameters would probably he 
termed an ‘improvement’. 


Мен 
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Figure 2 Definition of “Gun Jump” 


Two gun barrel configurations were examined and are shown in 
cvtline in Figure 2. 


The first configuration was used to examine tne effect or 
inserting a more compliant barrel section between the breech 
and the forward cradle bearing whilst retaining the same 
recoiling mags. The barrels, of course. no longer retain the 
same decree of balance. 


ex In almost every case Configuration No.l is the worst cate 


in d 

ық terms of the maximum disturbance amplitude, although clase te р 
shot exit (just over Утв), there is often little difference. 4 d 

Comparison ot barrel transverse velocity at the muzzle 15 s 


V-95 











eae Ab EM 


SRL CML UTI алша Адуч.а RINNE ESSE a ภา ร ดา ¿ ыу ыз WT N SE ee ee eee 


TAM с w. шул w mw re ee ee ee Ава mo m узшу oe ee -.- 


FERNY 2, KING 3 e 


shown in # 4; mete the difference іп the rates cf chance 
at shot exit. At “ms the cradle rotation is about dcuble that 
ст the Second configuration. 
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Figure 4 


n more detailed e:xanination of the results is com “inad T 
Conftiguration Ho, 2, 


= 


2.2 Forward cradle bearing position 


The forward cradle bearing position was increased in o. Em 
increments from іт to 2m. The muzzle transv rse velocity dzee 
not exhibit а marked change before shot exit (see Figure 5). 


The muzzle slope figures also show the master differences to 
eccur after shot exit (see Figure 6). The cradle rotation "INE 


reaches a maximum value earlier as the cradle bearing spacing 
15 increas d (from 9.5ms down to ms) and the ma: ru: um 
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amplitude 15 reduced from 9.8mr to 0. бте. At shot exit, 
howe.er, the cradle rotation increases with the cradle bearingc 
spacing since larger moments are being passed ints the 
cradle. no ๆ 
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Figure 6 Muzzle Slope 
Graph numbering Genotes cradle bearing spacing in metres 
5.2 Cradle bearing stitfness 
А ล number of cradle bearing stiftnesses were examined but only 
$3 three are shown. The units cf Dearing stiffness are Newtons 
" 


ner metre deflection per metre lerigth ct bearing. The figure 
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D 
2f 1.0: 107 Nm/m 15 often used and these results are ccompared 
with 0.95107 Mm/m and 5.9410’ Mm/m. 
Examination of the muzzle transverse velocity, muzzle slope 
and cradle rotation graphs indicated that the effect of the 
softer bearings was small, whereas the effect of the stiffer 
bearing was quite noticeable; in particular note the lower 
rates of change with the stiffer bearings. It should he noted 
that the stiffer bearings were much stiffer from the standard 
figure than the degree of softening. The results are shown in 
Figures 7, 8 and 9. 
1.0 
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Figure 7 Muzzle Transverse Velocity Б Ж 
Graph numbering denoted cradie bearing stiffness x1077N/m 
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Figure З Массе Slope Е ам 
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Figure 9 Cradle Rotation d 
Graph numbering denotes cradle bearing stiffness xit ° Nm/m 


The results are only true for the combinations of barrel and 
bearing stiffnesses together with the qun balance and bearing. 
locations. 


The computer model used for this work had linear elastic 
bearings with no backlash. The model indicates a substantially 
changed response with different stiffness cradle bearings. The 
proaram is in the process of incorporating a cradle bearing 
model which includes backlash. A paper by Bulman indicates the 
possible effects of bearing backlash 53). 


Since the cradle is considered as а rigid body the qustead 
bearing stiffness includes the effective cradle stiffness. 


The sutput variables chosen are not the only variables that 
can be examined and өсте projectiles are known to be sensiti.s 
to certain launch conditions. As Powell demonstrates, the 
projectile dynamics within the gun barrel can be important 
(4). However, based on our simple criteria the preferred 
system would be Configuration No.2 with the stiffer cradle 
bearings. 


It 15 considered that the computer model has improved sur 
appreciation of the dynamic interactions within the weapon 
system and that properly supported Бу triais results it will 
assist us in producing a more effective weapon system. 
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